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1.  INTRODUCTION 

The  goal  of  this  project  is  to  characterize  the  mechanisms  leading  to  hypermutated  prostate 
cancer  and  to  integrate  tumor  hypermutation  status  with  clinical  decision  making  and  therapy  to 
improve  the  care  of  men  with  advanced  prostate  cancer.  Using  Next-Gen  sequencing 
approaches  my  colleagues  at  the  University  of  Washington  recently  identified  a  hypermutated 
phenotype/genotype  in  10-20%  of  advanced  prostate  cancers.  This  phenotype  was 
subsequently  observed  in  primary  prostate  cancer.  Prostate  cancer  hypermutation  is  a 
promising  target  for  precision  therapy,  but  the  mechanisms  leading  to  hypermutation,  optimal 
methods  to  measure  hypermutation  status  in  the  clinic,  and  clinical  implications  for  prostate 
cancer  patients  are  not  yet  understood.  Our  hypothesis  is  that  hypermutated  advanced  prostate 
cancer  is  caused  by  defects  in  genes  regulating  DNA  repair  pathways,  which  can  be  accurately 
identified  using  existing  clinical  diagnostics,  and  that  hypermutation  status  can  predict 
responses  to  therapy. 


2.  KEYWORDS 

Prostate  cancer,  hypermutation,  hyper-mutation,  microsatellite  instability,  MSI,  MLH1,  MSH2, 
MSH6,  PMS2,  metastasis,  precision  medicine 
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3.  ACCOMPLISHMENTS 

Accomplishments  in  the  first  year  for  research-specific  tasks  are  reported  according  to  major 
goals  of  the  project  in  the  approved  SOW,  and  organized  by  specific  aim. 

3.1  What  were  the  major  goals  of  the  project? 


Specific  Aim  1:  Identify  mechanisms  that  drive  the 
hypermutated  phenotype  in  advanced  prostate  cancer. 

Months 

Completed  in 
Year  1? 

Major  Task:  Sequence  DNA  repair  pathway  genes  in 
advanced  prostate  cancer  tumor  samples 

1-12 

Yes 

Subtask  1:  Examine  hypermutated  and  non-hypermutated 
UW  prostate  cancer  rapid  autopsy  samples  using  BROCA 
and  UW-OncoPlex  assays 

1-6 

Yes 

Subtask  2:  Assess  for  functional  loss  of  DNA  repair 
pathway  gene  expression  by  IHC,  and  MSI  PCR 

3-12 

Yes 

Milestone(s)  Achieved:  identification  of  specific  mutated 
DNA  repair  pathway  genes  in  hypermutated  prostate 
cancer 

12 

Yes 

Specific  Aim  2:  Determine  unique  vulnerabilities  of 
hypermutated  prostate  cancer  to  therapy  in  xenograft 
models. 

Months 

Completed  in 
Year  1? 

Major  Task:  Assess  differential  responses  to 
chemotherapy  and  targeted  therapy  in  LuCaP  tumor  cell 
lines  xenografted  in  mice 

12-36 

No 

Subtask  1:  Use  xenograft  LuCaP  hypermutated  prostate 
cancer  cells  lines  58,  73,  and  147  and  3  non-hypermutated 
control  cell  lines.  Assess  xenograft  tumor  responses  to 
chemotherapy  and  targeted  therapy  agents 

12-36 

Partially 

Milestone(s)  Achieved:  Identification  of  differential  efficacy 
of  targeted  therapies  in  hypermutated  prostate  cancer 

24-36 

No 
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Specific  Aim  3:  Develop  and  validate  a  clinical  diagnostic 
approach  to  determine  hypermutation  status  in  advanced 
prostate  cancer. 

Months 

Completed  in 
Year  1? 

Major  Task:  Establish  a  clinical  assay(s)  to  detect  tumor 
hypermutation 

1-24 

Partially 

Subtask  1:  Develop  bioinformatics  methods  to  accurately  detect 
hypermutation  and  microsatellite  instability  using  the  UW- 
OncoPlex  assay 

1-12 

Yes 

Subtask  2:  Establish  the  performance  characteristics  of  MSI- 
PCR  and  IHC-based  approaches  to  detect  hypermutation 
compared  to  the  UW-OncoPlex  genomic  sequencing 

12-24 

No 

Milestone(s)  Achieved:  Clinically  validated  approach  to  detect 
the  hypermutated  subtype  of  advanced  prostate  cancer 
established 

24 

No 

Milestone(s)  Achieved:  Original  manuscript  on  bioinformatics 
method  on  detected  MSI  by  next-generation  sequencing 

12-24 

Yes 

Specific  Aim  4:  Implement  diagnostic  testing  for 
hypermutation  status  in  the  UW-OncoPlex  program  for 
precision  cancer  medicine. 

Months 

Completed  in 
Year  1? 

Major  Task:  Clinical  trial  of  UW-OncoPlex  testing  in  advanced 
prostate  cancer  that  includes  assessment  of  hypermutation 
status 

24-36 

No 

Subtask  1:  Establish  a  clinical  trial  that  includes  hypermutation 
testing  by  UW-OncoPlex  with  or  without  additional  MSI- 
PCR/MSI-IHC  tests  depending  on  results  of  Aim  3 

24-36 

No 

Subtask  2:  Report  hypermutation  status  results  to  medical 
oncologists  in  prostate  cancer  precision  tumor  board  meetings 
and  document  treatment  decisions  and  short-term  outcomes. 

24-36 

No 

Milestone(s)  Achieved:  Hypermutation  status  is  used  in  clinical 
decision  making  for  men  with  advanced  prostate  cancer  with 
feedback  on  outcomes 

36 

No 

Milestone(s)  Achieved:  Manuscript  describing  the  clinical  role  of 
tumor  hypermutation  status  as  a  predictive  biomarker  for 
advanced  prostate  cancer 

36 

No 
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3.12  What  was  accomplished  under  these  goals? 

Specific  Aim  1:  Identify  mechanisms  that  drive  the  hypermutated  phenotype  in  advanced 
prostate  cancer 

Work  on  Specific  Aim  1  was  largely  completed  in  Year  1  and  is  summarized  below.  We 
published  a  manuscript  in  Nature  Communications  based  on  the  work  accomplished  in  Aim  1 
(Pritchard  et  al.  Nat  Commun.  2014  5:4988,  see  Appendix  1). 

Specific  Aim  1,  Subtask  1:  Examine  hypermutated  and  non-hypermutated  UW prostate  cancer 
rapid  autopsy  samples  using  BROCA  and  UW-OncoPlex  assays 

We  hypothesized  that  mutations  in  key  DNA  repair  pathway  genes  lead  to  the  hypermutated 
subtype  of  advanced  prostate  cancer,  most  likely  mutations  in  DNA  mismatch  repair  genes.  To 
test  this  hypothesis  we  performed  targeted  deep  sequencing  of  DNA  repair  genes  in 
hypermutated  and  non-hypermutated  advanced  prostate  cancer  samples  from  two  sources: 
LuCaP  xenograft  lines  and  tumors  from  the  UWMC  rapid  autopsy  program.  Both  tumor  sources 
consisted  primarily  of  castration  resistant  prostate  cancer  (CRPC).  Using  exome  sequencing 
we  identified  3  hypermutated  patient-derived  xenograft  (PDX)  lines  (LuCaP  58,  LuCaP  73,  and 
LuCaP  147)  and  5  of  50  rapid  autopsy  patients  with  hypermutated  tumors  (05-165,  03-130,  06- 
134,  00-010,  05-123).  There  was  partial  overlap  between  the  PDX  and  the  autopsy  cases 
because  some  LuCaP  lines  had  been  derived  from  the  autopsy  patients.  There  were  a  total  of  7 
out  of  60  unique  patients  who  had  hypermutated  tumors,  for  an  overall  prevalence  of  11.6%  in 
our  cohort. 


We  performed  the  BROCA  targeted  DNA  capture  and  massively  parallel  sequencing  assay  that 
assesses  single  nucleotide  variants  (SNVs),  small  insertions  and  deletions  (indels),  copy 
number  variants  (CNVs),  and  structural  variants  (SVs)  in  DNA  repair  genes  simultaneously. 
Importantly,  the  BROCA  assay  includes  capture  of  complete  genes  including  introns  and 
flanking  sequences,  which  is  in  contract  to  exome  sequencing  which  captures  exons  only.  This 
detail  proved  to  be  crucial  to  our  success  in  this  research  aim.  We  sequenced  samples  to  an 
average  of  ~800x  depth,  multiplexing  24  samples  per  lane  on  a  HiSeq2500.  The  BROCA  assay 
uses  the  Agilent  SureSelect  enrichment  system  to  capture  the  coding  exons  and  flanking  splice 
sites  of  genes  listed  in  Table  1. 


Table  1:  BROCA 


genes  (assay  version  6) 


DNA 

Repair 

Pathways 

ATM 

CHEK1 

PALB2 

RBBP8 

ATR 

CHEK2 

PMS2 

TP53 

BAP1 

F AMU  5 A 

PRSS1 

TP53BP1 

BARD1 

MLH1 

PTEN 

XRCC2 

BRCA1 

MRE11A 

RAD50 

BRCA2 

MSH2 

RAD51B 

BRCC3 

MSH6 

RAD  51 C 

BRIP1 

NBN 

RAD51D 

Additional 

AKT1 

APC 

BMPR1A 

CDH1 

CDK4 

CDKN2A 

CTNNA1 

GALNT12 

Cancer- 

GEN1 

GREM1 

HOXB13 

MEN1 

MUTYH 

PIK3CA 

POLD1 

POLE 

Related 

PPM1D 

RET 

SDHB 

SDHC 

SDHD 

SMAD4 

STK11 

VHL 
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To  assess  within-patient  tumor  mutation  heterogeneity  we  tested  up  to  4  different  metastatic 
sites  in  a  subset  of  patients.  For  each  patient  we  also  tested  matched  normal  (non-tumor) 
tissue  to  determine  if  mutations  were  inherited  or  somatic. 

All  three  PDX  hypermutated  tumors  had  complex  structural  rearrangements  in  MSH2,  MSH6  or 
both  genes  (Table  2),  while  only  1  of  20  non-hypermutated  xenografts  had  mutations  in  these 
genes  (LuCaP  145,  derived  from  a  patient  with  neuroendocrine  prostate  cancer,  Supplementary 
Fig.  4).  A  second  loss-of-function  mutation  in  MSH2  or  MSH6  was  detected  in  the  three 
hypermutated  PDX  tumors,  but  not  in  LuCaP  145,  supporting  a  requirement  for  bi-allelic  gene 
inactivation  underlying  the  hypermutated  genome. 

We  performed  the  BROCA  targeted  DNA  capture  and  massively  parallel  sequencing  assays 
that  assesses  single  nucleotide  variants  (SNVs),  small  insertions  and  deletions  (indels),  copy 
number  variants  (CNVs),  and  structural  variants  (SVs)  in  DNA  repair  genes  simultaneously. 

All  three  PDX  hypermutated  tumors  had  complex  structural  rearrangements  in  MSH2,  MSH6  or 
both  genes  (Table  2),  while  only  1  of  20  non-hypermutated  xenografts  had  mutations  in  these 
genes  (LuCaP  145,  derived  from  a  patient  with  neuroendocrine  prostate  cancer).  A  second 
loss-of-function  mutation  in  MSH2  or  MSH6  was  detected  in  the  three  hypermutated  PDX 
tumors,  but  not  in  LuCaP  145,  supporting  a  requirement  for  bi-allelic  gene  inactivation 
underlying  the  hypermutated  genome. 


Table  2:  Mismatch  Repair  (MMR)  Gene  Mutations  Detected  in  All  Hypermutated  Prostate 
Cancers 


Tumor 

Source 

Hyper¬ 

mutated 

MSI 

MMR  Gene  Mutation(s)* 

LuCaP  58 

PDX 

Yes 

Yes 

1 )  MSH6  del  exon  8  to  3'UTR 

2)  MSH6  frameshift  (c.3799  3800del) 

LuCaP  73 

PDX 

Yes 

Yes 

1)  MSH2  and  MSH6  copy  loss  (del  3Mb) 

2)  MSH2-FBX011  inversion 

PDX 

1)  MSH2-KCNK12  74kb  inversion 

LuCaP  147** 

Yes 

Yes 

2)  MSH2-KCNK12  40kb  inversion 

3)  MSH2-C2orf61  inversion 

1)  MSH2-KCNK12  74kb  inversion 

05-165** 

Autopsy 

Yes 

Yes 

2)  MSH2-KCNK12  40kb  inversion 

3)  MSH2-C2orf61  inversion 

03-130 

Autopsy 

Yes 

Yes 

1)  MSH2  translocation  splits  the  gene  t(2;  1 8) 

2)  MSH2  copy  loss 

3)  MSH6  frameshift  (c.2690del) 

4)  MSH6  copy  loss 

06-134 

Autopsy 

Yes 

Yes 

MLH1  homozygous  copy  loss 

00-010 

Autopsy 

Yes 

Yes 

MSH2  frameshift  (c.2364  2365insTACA) 

05-123 

Autopsy 

Yes 

Yes 

1 )  MSH2  frameshift  (c.  1 1 24  11 25insG) 

2)  MSH2  frameshift  (c.1082del) 

3)  MLH1  frameshift  (c.1310del),  lymph  node  only 
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‘Mosaic  MSH6  frameshift  mutations  observed  in  a  poly  G  tract  in  exon  5  (c.3261dup/del)  and 
poly  A  tract  in  exon  7  (c.3573del)  were  detected  in  several  hypermutated  samples  and  are  not 
included  in  the  table  because  they  are  presumed  to  be  due  to  MSI. 

**LuCaP  147  is  derived  from  patient  05-165 


We  detected  mutations  with  predicted  loss-of-f unction  in  MSH2,  MSH6,  or  both  genes  in  4  of  5 
rapid  autopsy  patients  with  hypermutated  tumors.  Mutations  included  complex  structural 
rearrangements,  copy  losses,  and  frameshift  mutations  (Table  2,  Figure  1).  Two  hypermutated 
patients  had  mutations  in  the  MMR  gene  MLH1.  In  all  patients  hypermutation  status  and  MMR 
mutations  were  concordant  at  different  metastatic  sites  tested  in  the  same  patient.  MMR 
mutations  were  also  concordant  between  primary  tumor  and  metastasis  except  for  a  single 
MLH1  frameshift  mutation  in  patient  05-123  not  found  in  the  primary  (Table  2).  No  MMR 
mutations  were  detected  in  patient-matched  non-tumor  tumor  tissue,  indicating  that  that  none  of 
the  MMR  mutations  were  inherited  in  the  germline. 

A)  MSH2  Structural  Rearrangement  in  Hypermutated  Autopsy  Tumor  05-165  and  LuCaP  147 


Inversion  1  (343kb)  Inversion  2  (74kb)  Inversion  3  (40kb) 


C2orf61  intron  3  MSH2  exon  14  MSH2  exon  14  KCNK12  intron  1  MSH2  intron  15  KCNK1 2  intro n  1 


B)  MSH2  Structural  Rearrangement  in  Hypermutated 
Autopsy  T  umor  03-1 30 


Ht 


chr18 


i«i  MSH2  1 HHH — I — I  II I  mil 


chrl  8:45941 596 


chr2:47673506 


'  :■  »  -  :  m1- 


G  A  G  I  I  C  I  II  G  C  A  G  G  A  G 


chrl 8  q21.1  (intergenic)  MSH2  intron  8 


Figure  1:  Examples  of  complex 
MSH2  structural  rearrangement 
detected  in  hypermutated  prostate 
tumors.  A)  In  autopsy  sample  OS- 
165  and  patient-derived  xenograft 
LuCaP  147  is  a  representative 
complex  MSH2  rearrangement 
(LuCaP  147  was  derived  from 
autopsy  patient  05-165).  B)  MSH2 
structural  rearrangement  in 
hypermutated  autopsy  tumor  03-130. 
Breakpoints  were  confirmed  by 
Sanger  sequencing.  Genomic 
coordinates  are  build  hg19.  A  total 
of  4  or  7  hypermutated  cases  had 
complex  rearrangements  in  MSH2 
and  MSH6  or  both  genes. 
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To  cross-validate  mutation  calling  for  our  UW-OncoPlex  targeted  sequencing  platform  that  is  the 
focus  of  clinical  sequencing  work  for  precision  medicine  we  tested  one  hypermutated  rapid 
autopsy  prostate  cancer  case  (00-010)  and  two  non-hypermutated  autopsy  cases  (00-029  and 
00-090)  using  UW-OncoPlex.  We  also  tested  one  LuCaP  line  (LuCaP  23.1).  Among  the  genes 
that  overlap  the  two  panels  there  was  100%  concordance  of  somatic  coding  mutation  calls  that 
were  present  at  >5%  variant  allele  fraction  between  the  two  platforms. 


Specific  Aim  1,  Subtask  2:  Assess  for  functional  loss  of  DN A  repair  pathway  gene  expression 
by  IHC,  and  MSI  PCR 


MSH2  and  MSH6  are  mismatch  DNA  repair  genes  that  act  together  as  a  heterodimer,  and  bi- 
allelic  inactivating  mutations  of  either  gene  are  predicted  to  result  in  microsatellite  instability 
(MSI).  PCR  of  microsatellite  loci  revealed  MSI  in  all  hypermutated  tumors,  from  both  PDX  and 
autopsy  patients  (Figure  2,  Table  2).  IHC  for  DNA  mismatch  repair  proteins  in  hypermutated 
tumors  demonstrated  complete  loss  of  MSH2  and/or  MSH6  in  a  pattern  consistent  with  the 
inactivating  mutations  detected  by  sequencing  (Figure  3).  Non-hypermutated  tumors  were 
microsatellite  stable  and  had  intact  MSH2  and  MSH6  protein. 
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Figure  2:  Hypermutated  tumors 
are  MSI-High.  Hypermutated 
tumors  exhibited  microsatellite 
instability  by  PCR.  Shown  is 
representative  data  for  LuCaP  58 
which  is  positive  for  MSI  in  3/5 
mononucleotide  marker  systems 
(BAT25,  M0N027,  NR26,  arrows). 
All  hypermutated  tumors  were 
MSI-PCR  positive  in  at  least  2/5 
loci. 
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Figure  3:  Hypermutated  have  loss  of  MSH2  and  MSH6  protein  by  IHC.  Similar  results  were 
observed  in  hypermutated  tumors  from  rapid  autopsy  patients  (see  Appendix  1). 


The  findings  support  the  conclusion  that  the  hypermutated  subtype  of  prostate  cancer  is  chiefly 
due  to  loss-of-function  mutations  in  MSH2  and  MSH6  that  result  in  MSI.  Most  interestingly,  4  of 
7  hypermutated  cases  had  complex  structural  rearrangements  in  MSH2  and  MSH6  that  were 
not  detected  by  exome  sequencing  in  the  same  samples,  and  would  also  not  be  expected  to  be 
detected  by  traditional  exon-based  Sanger  sequencing  methods.  Previous  studies  have 
reported  MMR  protein  loss  and  MSI  in  both  primary  and  advanced  prostate  cancers,  but  very 
few  MMR  mutations  have  been  identified.  We  speculate  that  technical  limitations  have  led  to  an 
underestimation  of  MMR  gene  mutations  in  prostate  cancer. 


Specific  Aim  2:  Determine  unique  vulnerabilities  of  hypermutated  prostate  cancer  to 
therapy  in  xenograft  models. 

Aim  2,  Subtask  1:  Use  xenograft  LuCaP  hypermutated  prostate  cancer  cells  lines  58,  73,  and 
147  and  3  non-hypermutated  control  cell  lines.  Assess  xenograft  tumor  responses  to 
chemotherapy  and  targeted  therapy  agents 

Work  on  aim  2  will  largely  be  done  in  Year  2  and  3.  In  Year  1  IACUC  approval  was  obtained  to 
carry  out  animal  studies  and  preliminary  experiments  begun.  Note  that  no  funding  for  animal 
studies  is  provided  by  this  award.  The  goal  of  this  aim  is  to  carry  out  a  pilot  study  using  patient- 
derived  xenograft  (PDX)  preclinical  models  as  ‘tumor  avatars’  to  test  anti-cancer  therapies  in 
comparison  to  responses  in  non-hypermutated  LuCaP  xenograft  lines.  We  are  assessing 
responses  to  currently  used  and  approved  chemotherapeutics  including  docetaxel,  carboplatin 
and  5-fluorouricil,  to  determine  if  the  hypermutated  subtype  is  more  or  less  susceptible  to  drugs 
that  can  be  immediately  used  in  clinical  practice.  We  have  identified  3  hypermutated  PDX  lines 
(LuCaP  58,  73,  and  147)  that  we  are  using  to  assess  selective  responses  of  these  therapies,  in 
collaboration  with  Drs.  Colm  Morrissey,  Robert  Vessella,  and  Eva  Corey  at  the  University  of 
Washington  GU  Cancer  Research  Laboratory. 
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For  hypermutated  PDX  tumor  LuCaP  147  we  obtained  4  different  metastatic  sites  from  the 
patient  from  whom  the  xenograft  line  was  derived  and  found  that  the  same  complex  MSH2 
structural  rearrangements  were  present  in  all  metastatic  sites  in  the  pre-xenografted  tumors, 
demonstrating  that  the  MMR  gene  structural  rearrangements  are  not  an  artifact  of  xenografting. 

We  begun  by  testing  responses  to  docetaxel  because  it  is  the  most  commonly  used  genotoxic 
chemotherapy  in  advanced  prostate  cancer  patients.  In  collaboration  with  Drs.  Corey,  Vessella, 
and  Morrissey  we  evaluated  the  efficacy  of  docetaxel  in  hypermutated  LuCaP  xenografts 
(LuCaP  58,  73  and  147)  and  non-hypermutated  LuCaP  xenografts  (LuCaP  105CR,  96CR,  77, 
96  and  35CR).  LuCaP  tumors  were  subcutaneously  implanted  into  SCID  male  mice.  When 
tumor  exceeded  100mm3,  animals  were  randomized  and  enrolled  into  following  groups:  1) 
Docetaxel  treatment  at  10  mg/kg,  2)  Docetaxel  treatment  at  20  mg/kg,  3)  Vehicle  controls/no 
treatment  animals,  and  4)  Castration  animals.  Tumor  volumes  and  body  weight  were  measured 
once  weekly. 

Treatment  responses  for  the  two  different  dosages  (10  mg/kg  and  20  mg/kg)  varied  across  the 
hypermutated  and  non-hypermutated  tumor  models.  For  hypermutated  models,  10  mg/kg 
treatment  responses  ranged  from  major  tumor  growth  inhibition  in  LuCaP  58  to  mostly 
unimpeded  tumor  progression  in  LuCaP  147.  At  20  mg/kg,  LuCaP  73  exhibited  maximal 
responsiveness  as  opposed  to  LuCaP  58  and  147  (Figure  4).  We  also  assessed  single  animal 
response  to  docetaxel  (20  mg/kg)  within  the  hypermutated  LuCaP  PDX  and  found  heterogeneity 
in  responses.  Similar  to  hypermutated  LuCaP  PDX,  the  non-hypermutated  LuCaP  PDX 
exhibited  a  broad  range  of  susceptibility  to  docetaxel  (see  data-rich  PowerPoint).  Comparing 
the  responses  of  hypermutated  and  non-hypermutated  LuCaP  PDX,  we  did  not  find  a 
significant  difference  in  susceptibility  to  docetaxel,  or  a  differential  survival  benefit.  Our  results 
suggest  a  range  of  docetaxel  responsiveness  among  LuCaP  PDX  lines  that  is  not  strongly 
predicted  by  hypermutation  status. 

We  have  begun  similar  experiments  in  the  same  hypermutated  LuCaP  lines  using  carboplatin 
and  plan  to  test  5-FU  and  a  PARP  inhibitor.  We  anticipate  results  for  additional  drugs  in  Year  2 
and  3. 


LuCaP  73 


Control 

Castration 

Docetaxel  (10  mg/kg) 
Docetaxel  (20  mg/kg) 


Figure  4.  Responses 
to  docetaxel  in 
hypermutated  prostate 
cancer  PDX  LuCaP  73. 
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Specific  Aim  3:  Develop  and  validate  a  clinical  diagnostic  approach  to  determine 
hypermutation  status  in  advanced  prostate  cancer. 

Aim  3,  Subtask  1:  Develop  bioinformatics  methods  to  accurately  detect  hypermutation  and 
microsatellite  instability  using  the  UW-OncoPlex  assay 

Work  on  Aim  3,  has  been  a  focus  of  work  in  Year  1  and  will  continue  in  Year  2.  We  developed  a 
novel  method  for  inferring  MSI  and  hypermutation  from  next-generation  sequencing  data  that 
we  call  “mSINGS”.  We  recently  published  a  manuscript  on  this  method  for  which  Dr.  Pritchard 
was  the  senior  and  corresponding  author  (Salipante  et  al.  Clin  Chem.  2014  60:1192-9).  A 
graphical  depiction  of  how  the  mSINGS  method  works  is  given  in  Figure  5. 


Representative  MSI-PCR  Electropherograms  Representative  NGS  Virtual  Electropherograms 
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Figure  5:  Detection  of  microsatellite  instability  by  MSI-PCR  and  next-generation  DNA 
sequencing  using  “mSINGS”.  Representative  capillary  electrophoresis  results  from  MSI-PCR 
(top  panels)  and  “virtual  electropherograms”  of  next-generation  DNA  sequencing  data  (bottom 
panels),  where  the  length  (x-axis)  and  relative  abundance  (Y-axis)  of  variant  repeats  are  plotted. 
Loci  in  top  and  bottom  panels  are  not  equivalent,  and  are  from  different  genomic  locations. 

We  have  adapted  the  mSINGS  method  to  both  the  BROCA  and  UW-OncoPlex  genomic  deep 
sequencing  platforms  to  accurately  detect  both  phenotypic  MSI  and  hypermutation  status,  even 
when  matched  non-tumor  tissue  is  not  available  (Figure  6).  UW-OncoPlex  is  a  clinically- 
validated  diagnostic  platform  for  precision  cancer  medicine  developed  by  Dr.  Pritchard  that  has 
been  used  to  test  over  1,000  cancer  patients  to  date  (for  details  on  the  assay  see 
http://tests.labmed.washinqton.edu/UW-OncoPlex,  or  Google:  “UW-OncoPlex”).  We  have 
identified  65  mononucleotide  microsatellite  loci  that  are  captured  in  the  current  UW-OncoPlex 
assay  version  (version  4).  We  established  parameters  for  each  locus  to  be  called  unstable 
based  on  the  SD  of  peak  distribution  and  defined  MSI-High  as  having  at  20%  unstable  loci. 
Using  the  mSINGS  informatics  approach,  we  correctly  identified  all  known  MSI-High  cancer 
samples  (7/7)  including  one  hypermutated  prostate  cancer  sample  which  had  24/65  (37%)  loci 
unstable  (autopsy  sample  00-010,  liver  metastasis).  MSI-High  samples  had  35  +/-  12  unstable 
loci  (n=7,  mean  +/-  SD),  while  known  microsatellite  stable  samples  had  only  2  +/-  1 .5  unstable 
loci  (n=10). 
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Figure  6:  Detection  of  MSI  in  prostate  cancer  samples  using  mSINGS  applied  to  BROCA 
and  UW-OncoPlex  targeted  gene  sequencing  panels.  (Left  panel)  The  fraction  of  unstable 
microsatellite  loci  are  shown  for  BROCA  (left)  and  UW-OncoPlex  (right)  targeted  sequencing. 
Results  are  stratified  by  hypermutation  or  MSI  status.  The  threshold  used  for  interpreting  MSI 
status  is  indicated  by  a  dashed  line,  set  at  a  fraction  of  0.2  (20%  unstable  loci).  This  threshold 
perfectly  separated  hypermutated  (MSI  positive)  and  not  hypermutated  (MSI  negative)  tumors. 
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Figure  7.  Hypermutated  CRPC  cases  from 
SU2C  international  dream  team  have 
phenotypic  microsatellite  instability  (MSI) 
detected  by  mSINGS.  We  applied  an 
approach  to  measure  microsatellite  instability 
directly  from  next-generation  sequencing  data 
(mSINGS)  to  four  hypermutated  cases,  defined 
as  >300  nonsynonymous  mutations  in  exome 
sequencing.  A  threshold  fraction  of  0.2  (20%) 
unstable  loci  is  the  cutoff  for  microsatellite 
instability  using  the  mSINGS  method  (dashed 
line).  All  four  hypermutated  cases  were  MSI 
positive  and  had  somatic  mutations  in  mismatch 
repair  genes  ( MLH1 ,  MSH2 ,  or  MSH6). 

Selected  cases  with  less  than  300 
nonsynonymous  mutations  were  MSI  negative 
(bottom). 


In  Year  1  we  developed  a  simple  and  rapid  clinical  diagnostic  assay  based  on  the  mSINGS 
method  that  we  call  “MSIplus”.  This  method  using  amplicon  sequencing  of  17  microsatellite  loci, 
following  by  NGS.  We  validated  this  method  on  81  tumor  specimens  with  known  MSI  status, 
including  prostate  cancer  samples.  This  work  has  led  to  a  manuscript  which  is  currently  in 
press  in  the  Journal  of  Molecular  Diagnostics  (see  Appendix  4). 
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Aim  3,  Subtask  2:  Establish  the  performance  characteristics  of  MSI-PCR  and  IHC-based 
approaches  to  detect  hypermutation  compared  to  the  UW-OncoPlex  genomic  sequencing 

Work  on  this  subtask  has  not  yet  begun.  In  Year  2  we  plan  to  evaluate  MSI  detection  by 
MSIplus  and  UW-OncoPlex  in  hypermutated  and  non-hypermutated  prostate  cancer  samples 
that  have  been  previously  characterized  by  traditional  MSI-PCR.  Although  we  have  already 
shown  that  MSI  can  be  detected  in  hypermutated  prostate  cancer  samples  using  MSIplus,  we 
anticipate  that  a  subset  of  microsatellite  loci  will  be  optimal  for  detection  of  MSI  in  prostate 
cancer.  In  collaboration  with  Dr.  Stephen  Salipante  we  will  identify  a  set  of  optimal  loci  for  MSI 
detection  in  prostate  cancer  and  develop  bioinformatics  to  adapt  the  MSIplus  assay  to  evaluate 
these  loci  for  prostate  cancer  samples.  With  help  from  mentor  Dr.  Larry  True,  we  also  plan  to 
evaluate  IHC  patterns  in  hypermutated  tumors  to  determine  if  IHC  may  be  a  reliable  screening 
assay  to  identify  hypermutation. 


Specific  Aim  4:  Implement  diagnostic  testing  for  hypermutation  status  in  the  UW- 
OncoPlex  program  for  precision  cancer  medicine. 

Aim  4,  Subtask  1:  Establish  a  clinical  trial  that  includes  hypermutation  testing  by  UW-OncoPlex 
with  or  without  additional  MSI-PCR/MSI-IHC  tests  depending  on  results  of  Aim  3 

Work  on  Aim  4  is  in  early  stages  and  will  be  a  focus  of  work  in  Year  2  and  3.  No  funds  from  this 
award  have  yet  been  used  for  this  aim.  In  year  1  we  obtained  IRB  human  subjects  approval  for 
this  work.  We  have  begun  to  offer  clinical  UW-OncoPlex  testing  (Table  3)  for  prostate  cancer 
patients.  Clinical  reports  are  provided  and  results  discussed  directly  with  treating  oncologists 
and  urologists  at  a  monthly  precision  tumor  board  led  by  Dr.  Pritchard.  We  are  evaluating 
treatment  decision  making  in  prostate  cancer  patients  who  have  undergone  UW-OncoPlex 
testing.  We  plan  to  use  this  established  assay  and  tumor  board  framework  to  formally  test  the 
role  of  hypermutation  status  as  a  precision  biomarker  in  a  clinical  trial  in  Year  2  and  3  of  this 
research.  Recent  work  suggests  that  hypermutation  and  MSI  due  to  DNA  mismatch  repair 
deficiency  is  a  predictive  biomarker  for  anti-programmed  cell  death  1  (PD-1 )  immunotherapy  in 
several  cancers  (Le  et  al.  N  Engl  J  Med.  2015  372:2509-20).  In  collaboration  with  Dr.  Michael 
Schweizer,  we  are  working  on  establishing  a  protocol  to  screen  men  with  prostate  cancer  for 
MSI  using  MSIplus  as  a  qualifying  test  for  enrollment  on  anti-PD-1  therapeutic  trials.  This  work 
is  in  planning  stages. 


Aim  4,  Subtask  2:  Report  hypermutation  status  results  to  medical  oncologists  in  prostate 
cancer  precision  tumor  board  meetings  and  document  treatment  decisions  and  short-term 
outcomes. 

Work  on  this  subtask  will  largely  be  done  in  year  2  and  3.  We  have  already  identified  one 
patient  with  a  hypermutated  prostate  cancer  at  precision  tumor  board.  This  patient  had  MSI 
detected  by  our  mSINGS  method  using  UW-OncoPlex  and  an  underlying  tumor  MSH6  mutation 
with  associated  loss  of  heterozygosity.  The  discussion  at  precision  tumor  board  suggested  that 
the  patient  may  be  eligible  for  an  anti-PD-1  checkpoint  inhibitor  immunotherapy  trial.  Recent 
work  on  MSI-high  cancers  indicates  that  hypermutation  may  predict  response  to  anti-PDI 
therapy.  Unfortunately,  this  patient’s  disease  was  very  aggressive  and  the  patient  passed  away 
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before  trial  enrollment  could  be  accomplished.  In  year  2  and  3  we  anticipate  identifying  at  least 
3  to  4  additional  hypermutated  patients  through  a  100  patient  pilot  study. 


Table  3:  UW-OncoPlex™  genes  (assay  version  4) 


ABL1 

AKT1 

ALK 

AR 

ASXL1 

AURKA 

BAP1 

BCR 

BCL2L11 

Tier  1 : 

BRAF 

BRCA1 

BRCA2 

CCND1 

CCNE1 

CDK4 

CDK8 

CEBPA 

DDR2 

Currently 

DNMT3A 

EGFR 

EML4 

EPHB2 

ERBB2 

FGFR2 

FGFR4 

FLT3 

HI  FI  A 

IDH1 

IDH2 

JAK2 

KIF5B 

KIT 

KRAS 

MAP2K1 

MET 

MLL 

actionable 

MPL 

NKX2-1 

NPM1 

NRAS 

PDGFRA 

PIK3CA 

PML 

PTEN 

RARA 

NTRK1 

ROS1 

RET 

STK11 

TP53 

VHL 

PALB2 

CALR 

DNAJB1 

ABL2 

AKT2 

AKT3 

ATM 

AURKB 

BCOR 

CBL 

CBLB 

CDK6 

Tier  2: 

CHEK1 

CHEK2 

ERBB3 

ERBB4 

FBXW7 

FGFR1 

FLT1 

FLT4 

GATA2 

Actionable 

GNA11 

GNAQ 

GRM3 

HD  AC  4 

HRAS 

IGF1R 

IKZF1 

JAK3 

KDM6A 

in  the  near 
future 

KDR 

MAP2K2 

MAPK1 

MC1R 

MCL1 

MDM2 

MDM4 

MEN1 

MITF 

MLH1 

MRE11A 

MSH2 

MSH6 

MYC 

MYCN 

NF2 

NOTCH1 

PAX5 

PDGFRB 

PIK3R1 

PMS2 

RAF1 

RUNX1 

SMO 

SRSF2 

SUZ12 

TMPRSS2 

TSC1 

TSC2 

TET2 

TYR 

WT1 

POLE 

POLD1 

AXL 

ATRX 

BARD1 

BRIP1 

RAD51C 

RAD51D 

NBN 

FAM175A 

PHF6 

SHH 

GUI 

DAXX 

H3F3A 

SUFU 

TACSTD2 

APC 

BAK1 

BCL2 

CDH1 

CDKN2A 

CREBBP 

CRLF2 

CSF1R 

CTNNB1 

Tier  3: 

EPHA3 

EPHA5 

EPHB6 

ETV6 

EZH2 

FGFR3 

FOXA1 

GAB2 

GAT  A  3 

Frequently 

GATA1 

GNAS 

GRIN2A 

HNF1A 

IL7R 

JAK1 

MAP2K4 

MED12 

MUTYH 

MYCL1 

NF1 

NOTCH2 

PBRM1 

PRPF40B 

PTCH1 

PTPN11 

PTPRD 

RBI 

mutated 

RICTOR 

RPS14 

SF1 

SF3B1 

SMAD2 

SMAD3 

SMAD4 

SMARCA4 

SMARCB1 

SPOP 

SPRY4 

SRC 

TFG 

TGFBR2 

TRRAP 

U2AF1 

U2AF65 

ZRSR2 

ZBTB16 

CDKN1A 

CHD1 

MTAP 

MLH3 

PLK1 

PLK3 

PLK4 

PAK1 

TACC3 

DOCK7 

DEPDC5 

NPRL2 

NPRL3 

MIOS 

Germline 

ABCB1 

ABCC2 

ABCC4 

ABCG2 

C1orf144 

COMT 

CYP1B1 

CYP2C19 

CYP2C8 

CYP2D6 

CYP3A4 

CYP3A5 

DPYD 

EIF3A 

ERCC2 

ESR1 

ESR2 

FCGR1A 

pharmaco- 

FCGR2A 

FCGR3A 

GSTP1 

GUCY1A2 

ITPA 

LRP2 

MAN1B1 

MTHFR 

NQOI 

genomics 

NRP2 

UMPS 

SLC19A1 

SLC22A2 

SLC01B3 

SOD2 

SULT1A1 

TPMT 

TYMS 

UGT1A1 

http://tests.labmed.washington.edu/UW-OncoPlex 

This  is  a  clinically-available  comprehensive  gene  sequencing  platform  co-developed  and  offered 
clinically  by  Dr.  Pritchard’s  CLIA-certified  genetics  and  solid  tumors  laboratory. 
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3.3  What  opportunities  for  training  and  professional  development  has  the  project 
provided? 

Training-specific  tasks  from  the  approved  SOW  are  given  below.  Detail  related  to  training  goals 
in  the  first  year  is  provided  in  the  section  that  follows. 


Major  Task:  Training  and  educational  development  in 
prostate  cancer  research 

Months 

Completed  in 
Year  1? 

Subtask  1 :  Attend  the  Prostate  Cancer  Foundation  Annual 
Retreat  and  the  Association  of  Molecular  Pathology  Annual 
Conference 

1-36 

Yes,  see  National 
Conferences  and 
Committees 
section  below 

Subtask  2:  Present  research  at  the  monthly  mentor  group 
meetings,  and  at  least  once  per  year  at  Pacific  Northwest 
Prostate  Cancer  SPORE  research  conferences 

1-36 

Yes,  Seminars 
and  Interaction 
with  Mentors 
sections  below 

Subtask  3:  Lead  local  prostate  cancer  precision  tumor  board, 
including  review  of  genomic  sequencing  data  and  preparation 
of  presentations  that  integrate  prostate  cancer  patient  clinical 
histories  with  genomic  findings 

1-36 

Yes,  see  Prostate 
Precision  Tumor 
Board  section 
below 

Subtask  4:  Attend  face-to-face  meetings  as  part  of  the  Stand- 
Up-To-Cancer  (SU2C)  prostate  cancer  dream  team  2  to  3 
times/year.  Attend  monthly  conference  calls  for  the  SU2C 
prostate  cancer  dream  team  sequencing  and  analysis  group 

1-24 

Yes,  see  National 
Conferences  and 
Committees 
section  below 

Subtask  5:  Lead  local  “pipeline”  journal  club  focused  on 
application  of  new  genomic  technologies  in  the  clinic 

1-36 

Yes,  see  Journal 
Club  section 
below 

Subtask  6:  Train  senior  pathology  residents  and  fellows  in  the 
interpretation  and  clinical  reporting  of  prostate-cancer  focused 
precision  diagnostics 

1-36 

Yes,  see 
Teaching  section 
below 

Subtask  7:  Serve  as  PI  for  the  prostate  cancer  precision 
medicine  component  of  the  local  “ACT-SMART”  initiative  as 
part  of  the  institute  for  prostate  cancer  research 

1-36 

Yes 

Subtask  8:  Clinical  reporting  of  UW-OncoPlex  testing  applied 
to  advanced  prostate  cancer 

12-36 

Yes,  see  prostate 
cancer  precision 
tumor  board  and 
clinical  duties 
below 

Subtask  9:  Hands-on  training  in  prostate  cancer  with  Dr.  True  in 
pathology/immunohistochemistry  and  exposure  to  preclinical 

1-36 

Yes,  see  Hand’s- 
on  training 
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prostate  cancer  models  with  Drs.  Vessella  and  Morrissey. 

section  below 

Milestone(s)  Achieved:  Publication  of  original  research 

24-36 

Yes,  see 
appendices  for 
publications  in  the 
first  year 

Milestone(s)  Achieved:  Presentation  of  project  data  at  a 
national  meeting 

12-36 

Yes,  see  National 
Conferences  and 
Committees 
section  below 

National  Conferences  and  Committees:  In  year  1  I  attended  the  Prostate  Cancer  Foundation 
annual  retreat  where  I  presented  work  on  mechanisms  of  hypermutation  in  advanced  prostate 
cancer.  I  also  attended  the  Prostate  Cancer  Foundation  Coffey-Holden  Prostate  Cancer 
Academy  Meeting.  This  3-day  invite-only  meeting  was  a  think  tank  of  leading  prostate  cancer 
researchers  focused  on  the  question  of  oligo-metastatic  disease.  I  attended  the  Association  of 
Molecular  Pathology  Annual  Conference,  which  is  my  primary  professional  society.  This  year,  I 
was  invited  to  give  an  opening  plenary  session  at  the  2015  Association  of  Molecular  Pathology 
Annual  Conference  on  the  topic  of  bioinformatics  as  a  new  area  for  the  clinical  laboratory.  I  also 
attended  Stand  Up  to  Cancer  (SU2C)  Prostate  Cancer  International  Dream  Team  face-to-face 
meetings  and  participated  in  monthly  SU2C  sequencing  and  analysis  conference  calls.  As  part 
of  my  involvement  with  the  SU2C  Prostate  Cancer  International  Dream  Team  I  analyzed 
genomic  data,  prepared  figures  for  publication,  and  gave  a  formal  presentation  to  the  SU2C 
team  members. 

Seminars:  I  continue  to  attend  and  present  at  the  weekly  Pacific  Northwest  Prostate  Cancer 
SPORE  conference  series.  For  example,  at  the  most  recent  SPORE  talk  my  work  on  Aim  1  of 
this  research  was  presented  in  addition  to  my  collaborative  with  the  SU2C  international  dream 
team.  I  also  attend  weekly  Laboratory  Medicine  grand  rounds  and  Laboratory  Medicine 
research  rounds  which  include  a  wide  range  of  topics  related  to  clinical  diagnostic  medicine. 

Teaching:  In  year  1  of  funding  I  trained  4  senior  clinical  pathology  residents,  2  molecular 
genetic  pathology  fellows,  and  2  junior  molecular  pathology  faculty  (Dr.  Eric  Konnick  and  Dr. 
Tina  Lockwood)  in  the  clinical  interpretation  and  reporting  of  genomic  testing  for  prostate 
cancer.  I  have  continued  to  be  an  active  mentor  to  our  molecular  genetic  pathology  (MGP) 
fellowship  director,  MGP  fellows,  and  chief  resident  in  prostate  cancer-related  molecular 
diagnostics.  In  addition,  I  have  mentored  2  medical  genetics  senior  residents  in  molecular 
oncology  diagnostics,  including  one  resident  (Dr.  Mari  Tokita)  who  spent  9  months  of  dedicated 
time  in  my  research  laboratory  developing  circulating  tumor  DNA  diagnostics  methods.  As  part 
of  this  work  she  helped  to  assemble  ctDNA  samples  for  clinical  assay  validation  from  patients 
with  metastatic  prostate  cancer.  As  associated  director  of  the  genetics  and  solid  tumors 
laboratory  at  the  UW  I  have  also  mentored  6  junior  clinical  pathology  residents  in  4-week  basic 
genetics  training  rotations. 
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Journal  Club:  I  continue  to  lead  the  “pipeline”  monthly  journal  club  at  UW  which  is  attended  by 
about  20-30  faculty  and  senior  trainees  and  is  focused  on  genomic  technologies  applied  in  the 
clinic.  I  have  recruited  several  speakers  in  year  1  of  funding.  Examples  include  a  speaker 
(Mary  Goldman)  who  works  for  the  UCSC  genome  browser  and  outlined  new  tools  for  cancer 
data  visualization,  and  a  speaker  who  discussed  new  methods  for  gene  fusion  analysis  in 
cancer,  including  TMPRSS2-ERG  fusion  detection  in  prostate  cancer  (David  Wu). 

Prostate  Precision  Tumor  Board:  I  lead  the  local  prostate  precision  tumor  board  at  the 
University  of  Washington  where  we  review  genome  sequencing  data  from  advanced  prostate 
cancer  patients  and  make  treatment  recommendations  to  treating  oncologists.  This  activity  is  a 
substantial  effort  in  which  I  prepare  detailed  PowerPoint  presentations  for  each  patient  and  do  a 
thorough  literature  review  to  identify  potential  or  actual  therapeutic  implications  for  genomic 
findings.  This  monthly  tumor  board  is  attended  by  GU  medical  oncologists,  pathologists 
(including  co-mentor  Dr.  Larry  True),  urologists,  research  coordinators,  genetic  counselors,  and 
trainees  including  our  junior  faculty  and  molecular  genetic  pathology  fellows. 

Clinical  Duties:  I  am  a  primary  faculty  member  responsible  for  clinical  reporting  of  BROCA  and 
UW-OncoPlex  tests  that  are  available  through  my  CLIA-certified  laboratory.  I  have  already 
trained  2  other  junior  molecular  pathology  faculty  to  assist  me,  and  will  continue  to  mentor  and 
train  new  junior  faculty  to  assist  with  these  clinical  duties.  My  clinical  work  includes  close 
mentorship  from  Dr.  Mary-Claire  King,  who  personally  consults  on  all  BROCA  cases.  We  have 
weekly  “signout”  meetings  that  last  1  to  2  hours.  I  interpret  and  write  clinical  reports  for 
approximately  1 ,000  genomic  sequencing  cases  per  year.  As  part  of  my  clinical  work  I  interpret 
and  report  on  UW-OncoPlex  clinical  testing  for  all  prostate  cancer  patients  tested. 

Grant  Writing:  In  year  1,  in  collaboration  with  Dr.  Peter  Nelson  and  Dr.  Bruce  Montgomery,  I 
am  a  partnering  PI  for  a  DOD  IMPACT  proposal  focused  on  DNA  repair  defects  as  a  predictive 
biomarker  for  advanced  prostate  cancer.  I  am  also  preparing  an  R01  application  in 
collaboration  with  Heather  Hampel  at  The  Ohio  State  University  focused  on  genomic 
sequencing  as  a  tool  for  Lynch  syndrome  screening.  This  application  benefits  from  the  work  we 
have  done  to  develop  robust  methods  to  detect  hypermutation  and  MSI  in  cancer. 

Hand’s-on  training:  With  Dr.  True  I  have  received  personal  tutorials  in  prostate  cancer 
pathology,  including  one-on-one  formal  lectures  followed  by  teaching  slide  review.  In  addition, 
Dr.  True  provides  helpful  feedback  on  prostate  pathology-related  issues  at  the  monthly  prostate 
precision  tumor  board  that  I  lead.  I  have  also  received  ‘hand’s  on’  exposure  to  preclinical 
prostate  cancer  models  with  Drs.  Eva  Corey,  Robert  Vessella  and  Colm  Morrissey  in  the 
University  of  Washington  GU  Cancer  Research  Laboratory. 

Interaction  with  Mentors:  I  am  co-mentored  by  Dr.  Larry  True  and  Dr.  Mary-Claire  King  for 

this  training  award.  Dr.  True  is  an  internationally  prominent  genitourinary  pathologist  with  over 
30  years-experience  mentoring  junior  academic  researchers.  He  and  I  have  worked  together 
since  2000,  and  he  has  effectively  mentored  me  in  aspects  of  prostate  cancer  pathology,  and 
molecular  biomarker  development.  For  example,  I  served  on  the  NCI  Tissue-Based  biomarker 
subcommittee  of  the  Investigational  Drug  Steering  Committee  that  Dr.  True  chaired.  Dr.  King  is 
an  internationally  famous  geneticist  who  discovered  the  BRCA1  locus,  and  has  mentored 
dozens  of  highly  successful  faculty  in  the  area  of  cancer  genetics.  She  has  mentored  me  since 
1998  as  part  of  the  medical  scientist  training  program  and  has  been  an  even  closer  ongoing 
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mentor  to  me  over  the  past  3  years  in  the  clinical  implementation  of  the  BROCA  assay  that  was 
developed  in  her  laboratory  by  Dr.  Tom  Walsh.  I  meet  more  than  once  a  month  at  length  with 
Dr.  King  to  discuss  BROCA  test  results  and  discuss  research  directions,  and  I  meet  with  Dr. 
True  at  the  monthly  prostate  cancer  precision  tumor  board.  In  addition,  I  meet  with  Dr.  True  for 
formal  prostate  cancer  pathology  training. 

3.4  How  were  the  results  disseminated  to  communities  of  interest? 


Through  the  Institute  for  Prostate  Cancer  Research  (IPCR)  we  have  reached  out  to  patient 
advocates  in  the  region.  This  has  included  a  formal  presentation  on  prostate  cancer  precision 
medicine  that  I  gave  to  a  lay  audience  for  the  IPCR  on  the  topic  of  the  UW-OncoPlex  program 
for  prostate  cancer  precision  medicine.  In  year  2  we  plan  to  work  on  a  segment  to  air  on  UWTV 
that  highlights  the  prostate  cancer  precision  tumor  board  that  I  lead. 

3.5  What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 

Plans  in  Year  2 

Research-specific  tasks 

•  Evaluate  all  available  tumor  and  matched  normal  samples  from  the  UW  rapid  autopsy 
program  by  BROCA  to  identify  additional  hypermutated  cases. 

•  Continue  work  in  Aim  2  to  test  carboplatin  and  5-FU  therapies  in  hypermutated  PDX  tumors 
(no  funds  from  this  award  used  for  animal  studies). 

•  Implement  MSIplus  as  a  routine  clinical  assay  through  Dr.  Pritchard’s  CLIA-certified  lab. 
Evaluate  performance  characteristics  of  MSIplus  for  prostate  cancer  and  determine  the 
optimal  loci. 

•  Validate  mSINGS  analysis  by  UW-OncoPlex  for  clinical  use  and  reporting  in  prostate 
cancer. 

•  Test  additional  SU2C  international  dream  team  CRPC  cases  by  mSINGS  and  correlate  with 
MMR  mutation  status. 

•  Work  on  50  patient  pilot  project  using  clinical  UW-OncoPlex  for  men  with  prostate  cancer. 
Use  experience  gained  from  clinical  testing  to  inform  a  future  clinical  trial  for  men  with 
hypermutated  prostate  cancer. 

•  Evaluate  treatment  decision  making  in  prostate  cancer  patients  who  have  undergone  UW- 
OncoPlex  testing  through  review  with  treating  oncologists  at  the  precision  tumor  board.  The 
primary  measured  outcome  will  be  treatment  decisions  specifically  influenced  by 
hypermutation  status  results.  We  will  use  this  data  to  plan  future  clinical  trials  to  more 
formally  test  the  role  of  hypermutation  status  as  a  precision  biomarker. 

•  Investigate  protocols  to  refer  men  with  hypermutated  MSI  prostate  cancer  to  existing  clinical 
trials,  including  trials  of  anti  PD-1  immunotherapy. 

Training-specific  tasks 

•  Attend  PCF  annual  meeting  and  AMP  annual  meetings  as  well  as  SU2C  prostate 
international  dream  team  face-to-face  meetings. 

•  Present  data  at  SPORE  conference  and  at  national  conferences. 

•  Continue  resident  and  fellow  training  of  genomic  testing  in  prostate  cancer 

•  Continue  hand’s  on  training  meetings  with  Dr.  True  in  prostate  cancer  pathology 


20 


Colin  Pritchard  MD,  PhD  2014-2015  Year  1  Progress  Report  (AWARD:W81XWH-14-1-0448) 

Characterizing  the  Hypermutated  Subtype  of  Advanced  Prostate  Cancer  as  a  Predictive  Biomarker  for 
Precision  Medicine 

•  Continue  frequent  genomic  sequencing  signout  sessions  in  cancer  genetics  with  Dr.  Mary- 
Claire  King  (about  2  to  3  times  per  month). 

•  Continue  to  lead  monthly  Prostate  Precision  Tumor  Board 

•  Continue  to  lead  monthly  pipeline  journal  club 

•  Continue  grant  writing  activities 

4.  IMPACT 

4.1  What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 

The  work  in  year  1  of  this  research  has  led  to  the  discovery  of  the  mechanism  of  hypermutation 
in  advanced  prostate  cancer.  We  found  that  complex  somatic  MSH2  and  MSH6  mismatch  DNA 
repair  mutations  resulting  in  microsatellite  instability  are  the  chief  cause  of  hypermutation.  We 
also  found  that  hypermutation  is  more  common  in  advanced  prostate  cancer  than  previously 
expected,  with  7/60  (12%)  patients  identified  in  our  series.  Our  discovery  identifies  parallels  and 
differences  in  the  mechanisms  of  hypermutation  in  prostate  cancer  compared  with  other 
microsatellite  instability-associated  cancers.  Our  findings  have  important  implications  for 
prognosis  and  treatment.  If  hypermutation  can  be  targeted,  a  substantial  minority  of  patients 
with  advanced  prostate  cancer  may  benefit.  For  example,  cancers  with  mismatch  DNA  repair 
deficiency  have  recently  been  shown  to  be  responsive  to  anti-PD-1  immunotherapy.  Our 
research  has  also  facilitated  microsatellite  instability  and  immunohistochemistry-based  testing 
as  screening  tools  for  hypermutation  in  advanced  prostate  cancer. 

Our  work  in  year  1  has  also  led  to  the  development  of  highly  innovative  and  robust  methods  to 
detect  microsatellite  instability  (MSI)  that  is  associated  with  hypermutation.  We  recently 
developed  the  “mSINGS”  method  for  detection  of  MSI  directly  from  NGS.  This  has  facilitated 
MSI  analysis  of  exome  data  from  the  SU2C  prostate  cancer  international  dream  team 
consortium,  proving  that  all  hypermutated  prostate  cancer  cases  in  that  series  also  have  MSI 
associated  with  underlying  mismatch  DNA  repair  mutations. 

4.2  What  was  the  impact  on  other  disciplines? 

Our  work  builds  bridges  between  research  in  colorectal  and  endometrial  cancer  and  research  in 
prostate  cancer.  Hypermutation  and  MSI  are  well-studied  in  colorectal  and  endometrial  cancer. 
We  have  applied  the  mSINGS  method  we  developed  to  both  colorectal  and  endometrial  cancer, 
resulting  in  the  first  ever  tumor-based  DNA  sequencing  test  for  Lynch  syndrome,  ColoSeq 
Tumor.  We  collaborate  closely  with  colleagues  at  the  Ohio  State  University  on  Lynch  syndrome 
screening  research  (Heather  Hampel  and  Albert  de  la  Chapelle),  resulting  in  a  recent  NIH  R01 
grant  submission  that  harnesses  the  mSINGS  method. 

4.3  What  was  the  impact  on  technology  transfer? 

Nothing  to  Report 

4.4  What  was  the  impact  on  society  beyond  science  and  technology? 

Nothing  to  Report 
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5.  CHANGES/PROBLEMS 

5.1  Changes  in  approach  and  reasons  for  change 

Due  to  cost  constraints  we  were  required  to  slightly  modify  our  planned  experimental  design  for 
aim  2  to  evaluate  only  3  different  therapies  in  the  hypermutated  and  non-hypermutated  LuCaP 
xenograft  models.  The  three  therapies  chosen  are  docetaxel,  carboplatin,  and  5-FU.  Note  that 
this  change  does  not  affect  the  budget  of  this  award  because  no  vertebrate  animal  work  is 
funded  through  this  training  award. 

5.2  Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them 

Based  on  prior  experience  by  the  GU  cancer  research  laboratory  we  anticipate  potential 
problems  with  carboplatin  toxicity  in  future  LuCaP  xenograft  animal  studies.  This  is  particularly 
a  problem  for  the  LuCaP  58  hypermutated  line  in  which  frequent  ulceration  may  occur.  To 
address  this,  we  will  lower  the  dose  and  dosing  frequency  of  carboplatin  for  this  arm  of  the 
study. 

5.3  Changes  that  had  a  significant  impact  on  expenditures 

Nothing  to  report. 

5.4  Significant  changes  in  use  or  care  of  human  subjects,  vertebrate  animals,  biohazards, 
and/or  select  agents 

Nothing  to  report. 


6.  PRODUCTS 

6.1  Publications,  conference  papers,  and  presentations 
6.11  Journal  publications 

Pritchard  CC  (corresponding  author),  Morrissey  C,  Kumar  A,  Zhang  X,  Smith  C,  Coleman  I, 
Salipante  SJ,  Milbank  J,  Tait  JF,  Corey  E,  Vessella  RL,  Walsh  T,  Shendure  J,  Nelson  PS; 
Complex  MSH2  and  MSH6  Mutations  in  Hypermutated  Microsatellite  Unstable  Advanced 
Prostate  Cancer;  Nature  Communications',  5:  2014;  4988;  published;  acknowledgement  of 
federal  support  (yes). 

See  Appendix  1  and  2 

Robinson  D,  Van  Allen  EM,  Wu  Y,  Schultz  N.,  Lonigro  RJ,  Mosquera  J,  Montgomery  R,  Taplin 
ME,  Pritchard  CC  (co-second  author),  Attard  G,  Beltran  H,  Abida  WM,  Bradley  RK,  Vinson  J, 
Cao  X,  Vats  P,  Kunju  LP,  Hussain  M,  Feng  FY,  Tomlins  SA,  Cooney  KA,  Smith  DC,  Brennan  C, 
Siddiqui  J,  Mehra  R,  Scher  HI,  Chen  Y,  Rathkopf  DE,  Morris  MJ,  Solomon  SB,  Durack  JC, 
Reuter  VE,  Gopalan  A,  Gao  J,  Loda  M,  Lis  RT,  Bowden  M,  Balk  SP,  Gaviola  G,  Sougnez  C, 
Gupta  M,  Yu  EY,  Mostaghel  EA,  Cheng  HH,  Chew  FS,  True  LD,  Plymate  SR,  Dvinge  H, 
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Ferraldeschi  R,  Flohr  P,  Miranda  S,  Zafeiriou  Z,  Tunariu  N,  Mateo  J,  Demichelis  F,  Elemento  O, 
Robinson  BD,  Sboner  A,  Schiffman  MA,  Nanus  DM,  Tagawa  ST,  Sigaras  A,  Eng  KW,  Heath  E, 
Pienta  KJ,  Kantoff  P,  de  Bono  JS,  Rubin  MA,  Nelson  PS,  Garraway  LA,  Sawyers  CL, 
Chinnaiyan  AM;  Integrative  clinical  sequencing  analysis  of  metastatic  castration  resistant 
prostate  cancer  reveals  a  high  frequency  of  clinical  actionability;  Cell.  161 :  2015;  1215-1228; 
published;  acknowledgement  of  federal  support  (yes). 

See  Appendix  3 

Hempelmann  JA,  Scroggins  SM,  Pritchard  CC,  Salipante  SJ;  MSIplus:  integrated  colorectal 
cancer  molecular  testing  by  next-generation  sequencing;  Journal  of  Molecular  Diagnostics. 
2015;  in  press;  acknowledgement  of  federal  support  (yes). 

See  Appendix  4 

6.12  Books  or  other  non-periodical,  one-time  publications. 

Nothing  to  report. 

6.13  Other  publications,  conference  papers,  and  presentations. 

‘Pritchard  CC,  Morrissey  C,  Kumar  A,  Zhang  X,  Smith  C,  Coleman  I,  Salipante  S,  Grady  WM, 
Tait  JF,  Vessella  R,  Walsh  T,  Shendure  J,  and  Nelson  PS.  Mechanisms  of  Microsatellite 
Instability  in  Hypermutated  Advanced  Prostate  Cancer.  (2014)  Prostate  Cancer  Foundation 
Annual  Scientific  Retreat. 

‘Salipante  S,  Scroggins  S,  Hampel  HL,  Turner  EH,  and  Pritchard  CC.  Microsatellite  Instability 
Detection  By  Next-Generation  Sequencing.  (2014)  Academy  of  Clinical  Laboratory  and 
Physician  Scientists  Annual  Meeting. 

6.2  Website(s)  or  other  Internet  site(s) 

Nothing  to  report. 

6.3  Technologies  or  techniques 

We  have  developed  the  mSINGS  method  for  detection  of  microsatellite  instability  from  targeted 
next-generation  sequencing  data.  The  source  code  for  this  bioinformatics  method  is  freely 
available  for  academic  users  and  can  be  found  at:  https://bitbucket.org/uwlabmed/msinqs.  This 
source  code  has  already  been  shared  with  multiple  national  and  international  researchers. 

6.4  Inventions,  patent  applications,  and/or  licenses 
Nothing  to  report. 

6.5  Other  Products 


23 


Colin  Pritchard  MD,  PhD  2014-2015  Year  1  Progress  Report  (AWARD:W81XWH-14-1-0448) 

Characterizing  the  Hypermutated  Subtype  of  Advanced  Prostate  Cancer  as  a  Predictive  Biomarker  for 
Precision  Medicine 

The  genomic  sequencing  dataset  we  generated  as  part  of  the  research  on  aim  1  to  identify 
mechanisms  of  hypermutation  is  publically  available  at  GenBank/EMBL/DDBJ  under  the 
accession  code  SRP044943. 

7.  PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 

7.1  What  individuals  have  worked  on  the  project? 

Name:  Colin  C.  Pritchard 
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Researcher  Identifier:  ORCID  ID:  0000-0002-2461-1557 
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written  manuscripts  related  to  this  work 
Funding  Support:  (this  award) 

Name:  Robert  Livingston 
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Researcher  Identifier:  ORCID  ID:  0000-0002-7424-2956 

Nearest  person  month  worked:  1 

Contribution  to  Project:  Dr.  Livingston  has  coordinated  coordinating  genomic  sequencing  library 
preparation  and  assisted  in  data  analysis. 

Funding  Support:  (this  award) 

Name:  Christina  Smith 

Project  Role:  Senior  Clinical  Technologist 

Nearest  person  month  worked:  3 

Contribution  to  Project:  Ms.  Smith  has  helped  perform  genomic  sequencing  assays  including 
BROCA  and  UW-OncoPlex. 

Funding  Support:  (this  award) 

Name:  Mallory  Beightol 
Project  Role:  Research  Technician 
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Contribution  to  Project:  Ms.  Beightol  has  helped  perform  genomic  sequencing  assays  including 
BROCA  and  UW-OncoPlex. 
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7.2  Has  there  been  a  change  in  the  active  other  support  of  the  PD/PI(s)  or  senior/kev 
personnel  since  the  last  reporting  period? 

Colin  Pritchard:  Changes  in  current  support  since  last  reporting  period: 

RAF-CA1 4-004  (R33) 

Title:  Advanced  Development  Advanced  Development  and  validation  of  targeted  molecular 
counting  methods  for  precise  and  ultrasensitive  quantitation  of  low  prevalence  somatic 
mutations  (Salipante,  PI) 

Time  Commitment:  0.6  Calendar  Months  (5%  effort) 

Role:  Co-investigator 

Supporting  Agency:  NIH 

Address  of  Funding  Agency:  grants.gov 

Performance  Period  4/1/15-3/31/18 

Level  of  Funding:  $295,868 

Project  Goals:  The  ultrasensitive  detection  of  clinically  relevant  somatic  alterations  in  cancer 
genomes  has  great  potential  for  impacting  patient  care,  e.g.  for  early  detection,  establishing 
diagnoses,  refining  prognoses,  guiding  treatment,  and  monitoring  recurrence.  However,  current 
technologies  are  poorly  suited  to  the  robust  detection  of  somatic  mutations  present  at  very  low 
frequencies  (<1%).  Massively  parallel  sequencing  represents  an  advantageous  path  forward, 
but  its  sensitivity  to  detect  very  rare  events  is  fundamentally  constrained  by  the  sequencing 
error  rate.  We  have  recently  developed  a  new  experimental  paradigm  that  overcomes  this 
limitation.  In  our  approach,  each  copy  of  a  target  sequence  that  is  present  in  a  sample  is 
molecularly  tagged  during  the  first  cycle  of  a  multiplex  capture  reaction  with  a  unique  random 
sequence.  After  amplification,  target  amplicons  and  their  corresponding  molecular  tags  are 
subjected  to  massively  parallel  sequencing.  During  analysis,  the  molecular  tags  are  used  to 
associate  sequence  reads  sharing  a  common  origin.  The  availability  of  robust,  cost-effective, 
quantitative,  and  generically  applicable  tools  for  the  ultrasensitive,  multiplex  detection  of  rare 
somatic  events  in  the  clinical  setting  will  provide  enhanced,  transformative  capabilities  in  the 
diagnosis  and  monitoring  of  cancers.  The  methodology  will  also  have  application  to  basic 
science  cancer  research. 

No  overlap 


7.3  What  other  organizations  were  involved  as  partners? 

Nothing  to  report. 


8.  SPECIAL  REPORTING  REQUIREMENTS 

Nothing  to  report. 
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9.  APPENDICES 

Appendix  1 :  Reprint  of  Pritchard  CC  et  al.;  Complex  MSH2  and  MSH6  Mutations  in 
Hypermutated  Microsatellite  Unstable  Advanced  Prostate  Cancer;  Nature  Communications',  5: 
2014;4988. 

Appendix  2:  Supplemental  Figures  from  Pritchard  CC  et  al.;  Complex  MSH2  and  MSH6 
Mutations  in  Hypermutated  Microsatellite  Unstable  Advanced  Prostate  Cancer;  Nature 
Communications',  5:  2014;  4988. 

Appendix  3:  Reprint  of  Robinson  D  et  al.;  Integrative  clinical  sequencing  analysis  of  metastatic 
castration  resistant  prostate  cancer  reveals  a  high  frequency  of  clinical  actionability;  Cell.  161: 
2015;  1215-1228. 

Appendix  4:  Draft  of  Hempelmann  JA  et  al.;  MSIplus:  integrated  colorectal  cancer  molecular 
testing  by  next-generation  sequencing;  Journal  of  Molecular  Diagnostics.  2015;  in  press. 
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Recently  exome  sequencing  of  metastatic  prostate  cancers 
revealed  that  a  subset  of  patients  harboured  tumors 
with  markedly  elevated  single- nucleotide  mutation  rates, 
defining  a  new  hypermutated  subtype1.  This  phenotype 
was  subsequently  observed  in  primary  prostate  cancer  in  a 
tumour  that  harboured  an  MSH6  mutation2.  However, 
mechanisms  that  lead  to  hypermutation  and  the  prevalence  of 
this  distinct  subtype  have  not  been  completely  defined. 
Comprehensive  cancer  genomics  efforts  recently  published  by 
The  Cancer  Genome  Atlas  Research  Network  (TCGA)  reported 
that  16%  of  colon  cancers  and  up  to  35%  of  endometrial  cancers 
exhibit  hypermutation3,4.  For  both  colon  and  endometrial 
cancers,  about  three  quarters  of  hypermutated  tumors  were 
associated  with  phenotypic  microsatellite  instability  (MSI)  and 
loss-of-function  DNA  mismatch  repair  genes  via  mutation  or 
epigenetic  silencing.  Therefore,  we  hypothesized  that 
hypermutated  prostate  cancer  may  also  be  associated  with  DNA 
mismatch  repair  (MMR)  gene  defects  and  MSI. 

In  this  study  we  identified  hypermutation  in  7  of  60  patients 
with  advanced  prostate  cancer.  Using  a  targeted  deep  sequencing 
approach  we  find  that  all  hypermutated  tumors  have  somatic 
mutations  in  MMR  genes  and  associated  MSI.  In  four  of 
seven  hypermutated  cases  MMR  mutations  were  complex 
structural  rearrangements  in  MSH2  and  MSH6.  We  conclude 
that  somatic  rearrangements  in  MSH2  and  MSH6  are  an 
important  mechanism  leading  to  hypermutation  and  MSI  in 
advanced  prostate  cancer. 

Results 

Prevalence  of  hypermutation.  We  identified  hypermutated  cases 
in  exome  sequencing  data  sets  of  advanced  prostate  cancer 
samples  from  two  sources:  a  panel  of  patient-derived  xenografts 
(PDX)  and  metastatic  specimens  obtained  through  a  rapid 
autopsy  programme  (Supplementary  Table  1).  Exome  data  for 
PDX  tumors  was  from  Kumar  et  al.1 ,  where  hypermutation  was 
previously  characterized.  In  the  autopsy  samples  where 
hypermutation  status  had  not  been  previously  established,  we 


defined  hypermutation  as  >  300  somatic  protein  altering 
mutations  based  on  the  distribution  of  total  mutation  burden  in 
metastatic  tumors,  which  had  matched  normal  tissue  available 
(Supplementary  Fig.  1;  Supplementary  Table  1).  We  identified 
hypermutation  in  3  of  15  PDX  tumors  (Table  1),  and  in 
metastatic  tumors  from  5  of  50  autopsy  patients  (Table  2).  There 
was  partial  overlap  between  the  two  patient  groups:  five  of  the 
PDX  tumors  were  derived  from  autopsy  patients,  including  one 
with  a  hypermutated  genome  (LuCaP  147).  Therefore,  there  were 
a  total  of  7/60  unique  patients  with  hypermutated  tumors,  for  an 
overall  prevalence  of  11.6%.  Hypermutation  status  was  100% 
concordant  at  different  metastatic  sites,  and  was  also  concordant 
between  primary  tumour  and  metastasis  in  two  patients  where 
primary  prostate  tumors  were  available  (Table  2). 


Identification  of  MSH2  and  MSH6  rearrangements.  Because 
exome  sequencing  has  limitations  in  detecting  structural  rear¬ 
rangements  and  larger  insertion/ deletion  (indel)  mutations,  we 
investigated  alterations  in  DNA  MMR  pathway  genes  in  hyper¬ 
mutated  and  non-hypermutated  cases  using  a  targeted  deep 
sequencing  approach  (BROCA  assay)  that  included  capture  of 
intronic  and  flanking  DNA  sequences  (Supplementary 
Table  2)5,6.  We  developed  a  bioinformatics  pipeline  to 
accurately  detect  structural  variation,  copy  number  variation 
and  indel  mutations  of  all  sizes7. 

All  three  PDX  hypermutated  tumors  had  complex  structural 
rearrangements  in  MSH2 ,  MSH6  or  both  genes  (Table  1;  Fig.  la; 
Supplementary  Figs  2-4),  while  only  1  of  20  non-hypermutated 
xenografts  had  mutations  in  these  genes  (FuCaP  145, 
derived  from  a  patient  with  neuroendocrine  prostate  cancer, 
Supplementary  Fig.  5).  A  second  loss-of-function  mutation  in 
MSH2  or  MSH6  was  detected  in  the  three  hypermutated  PDX 
tumors,  but  not  in  FuCaP  145,  supporting  a  requirement  for  bi- 
allelic  gene  inactivation  underlying  the  hypermutated  genome. 

We  detected  mutations  with  predicted  loss-of-function  in 
MSH2 ,  MSH6  or  both  genes  in  four  of  five  rapid  autopsy  patients 


Table  1  |  MMR  gene  mutations  in  prostate  cancer  PDX. 


PDX  tumour* 

Patient-derived  from 

Hypermutated?^ 

MSI 

MMR  gene  mutation(s)^ 

LuCaP  58 

Yes 

Yes 

(1)  MSH6  del  exon  8  to  3'UTR 

(2)  MSH6  frameshift  (c.3799_3800del) 

LuCaP  73 

Yes 

Yes 

(1)  MSH2  and  M5H6  copy  loss  (del  3  Mb) 

(2)  MSH2-FBX011  inversion 

LuCaP  147,  147CR 

05-165 

Yes 

Yes 

(1)  MSH2-C2orf61  343  kb  inversion 

(2)  MSH2-KCNK12  74  kb  inversion 

(3)  MSH2-KCNK12  40  kb  inversion 

LuCaP  23.1,  23. ICR 

No 

No 

None 

LuCaP  35,  35CR 

No 

No 

None 

LuCaP  70,  70CR 

No 

No 

None 

LuCaP  77,  77CR 

No 

No 

None 

LuCaP  78 

98-328 

No 

No 

None 

LuCaP  81 

98-362 

No 

No 

Chr2  copy  losses 

LuCaP  86.2,  86.2CR 

No 

No 

None 

LuCaP  92 

99-069 

No 

No 

None 

LuCaP  96,  96CR 

No 

No 

None 

LuCaP  105,  105CR 

No 

No 

None 

LuCaP  141 

No 

No 

None 

LuCaP  145.1,  145.2 

05-144 

No 

No 

(1)  MSH2  exon  8-16  del 

(2)  MSH6-TESC  t(2;12) 

MMR,  mismatch  repair;  MSI,  microsatellite  instability;  PDX,  patient-derived  xenografts. 

^Matched  pairs  of  androgen-sensitive  and  castration-resistant  sublines  (for  example,  LuCaP  35  and  LuCaP  35CR)  and  tumour  lines  derived  from  the  same  patient  are  listed  numerically  and  grouped  in  the 
same  row. 

tHypermutation  status  was  previously  determined  in  these  samples  in  Kumar  et  a\) 

^Mosaic  MSH6  frameshift  mutations  observed  in  a  poly  G  tract  in  exon  5  (c.3261dup/del)  and  poly  A  tract  in  exon  7  (c.3573del)  were  detected  in  several  hypermutated  samples  and  are  not  included  in 
the  table  because  they  are  presumed  to  be  due  to  MSI. 
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Table  2  |  MMR  gene  mutations  in  rapid  autopsy  patients. 


Autopsy 

patient* 

Tumour  site(s)  tested  by 

BROCA  targeted  sequencing 

Mutation  burden^ 
(exome) 

Hypermutated? 

MSI 

MMR  gene  mutation(s)^ 

05-165* 

Bone,  adrenal,  liver  and  lymph 
node 

855 

Yes 

Yes 

(1)  MSH2-C2orf61  343  kb  inversion 

(2)  MSH2-KCNK12  74  kb  inversion 

(3)  MSH2-KCNK12  40  kb  inversion 

03-130 

Lymph  node 

647 

Yes 

Yes 

(1)  MSH2  translocation  splits  the  gene  t(2;18) 

(2)  MSH2  copy  loss 

(3)  MSH6  frameshift  (c.2690del) 

(4)  MSH6  copy  loss 

06-134 

Kidney  and  lymph  node 

314 

Yes 

Yes 

MLH1  homozygous  copy  loss 

00-010 

Prostate  and  liver 

673 

Yes 

Yes 

MSH2  frameshift  (c.2364_2365insTACA) 

05-123 

Prostate  and  lymph  node 

807 

Yes 

Yes 

(1)  MSH2  frameshift  (c.1124_1125insG) 

(2)  MSH2  frameshift  (c.1082del) 

(3)  MLH1  frameshift  (c.131  Odel),  lymph  node 
only 

01-095 

Liver  and  lymph  node 

149 

No 

No 

None 

05-144* 

Bone,  adrenal,  liver  and  lymph 
node 

57 

No 

No 

(1)  M5H2  exon  8-16  del 

(2)  MSH6-TESC  t(2;12) 

05-214 

Bone,  liver  and  lymph 
node  (two  sites) 

46 

No 

No 

None 

05-116 

Bone,  adrenal,  liver  and  lung 

47 

No 

No 

None 

00-029 

Liver 

37 

No 

No 

None 

00-090 

Lymph  node 

69 

No 

No 

None 

MMR,  mismatch  repair;  MSI,  microsatellite  instability. 

*Fifty  total  unique  autopsy  patients  were  assessed  by  exome  sequencing  (see  Supplementary  Table  1).  Listed  are  a  subset  of  cases  that  were  followed  up  by  targeted  deep  sequencing  for  MMR  genes. 
Clinical  data  for  this  patient  subset  is  provided  in  Supplementary  Table  6.  Patient-matched  non-cancer  tissue  was  tested  in  every  case  and  did  not  exhibit  MSI  or  MMR  mutations.  LuCaP  147  and  147CR 
are  derived  from  autopsy  patient  05-165.  LuCaP  145.1  and  145.2  are  derived  from  autopsy  patient  05-144. 

tNumber  of  protein  altering  somatic  mutations  by  exome  sequencing  with  removing  of  germline  variants  from  matched-non-tumour  samples. 

^Mutations  were  detected  at  every  tumour  site  unless  otherwise  indicated.  Mosaic  MSH6  frameshift  mutations  observed  in  a  poly  G  tract  in  exon  5  (c.3261dup/del)  and  poly  A  tract  in  exon  7  (c.3573del) 
were  detected  in  several  hypermutated  samples  and  are  not  included  in  the  table  because  they  are  presumed  to  be  due  to  MSI. 


with  hypermutated  tumors.  Mutations  included  complex  struc¬ 
tural  rearrangements,  copy  losses  and  frameshift  mutations 
(Table  2;  Supplementary  Figs  4  and  6-9).  Two  hypermutated 
patients  had  mutations  in  the  MMR  gene  MLH1.  We  interrogated 
a  subset  of  six  non-hypermutated  patients  by  deep  sequencing 
and  did  not  detect  MMR  gene  mutations  except  in  patient  05-144 
from  which  the  PDX  LuCaP  145  was  derived  (Table  2).  Like 
hypermutation  status,  MMR  mutations  were  concordant  at 
different  metastatic  sites  in  the  same  patient.  MMR  mutations 
were  also  concordant  between  primary  tumour  and  metastasis 
except  for  a  single  MLH1  frameshift  mutation  in  patient  05-123 
not  found  in  the  primary  tumour  (Table  2;  Supplementary  Fig.  9). 
Patient-matched  non-tumour  tissues  were  tested  for  the  autopsy 
patients  (Supplementary  Table  1  and  Supplementary  Data  1).  No 
MMR  mutations  were  detected  in  patient-matched  non-tumour 
tissue,  indicating  that  none  of  the  MMR  mutations  were  inherited 
in  the  germline.  Mutations  in  additional  DNA  repair  genes  are 
given  in  Supplementary  Table  3. 

Hypermutated  tumors  have  phenotypic  MSI.  MSH2  and  MSH6 
are  mismatch  DNA  repair  genes  that  act  together  as  a  hetero¬ 
dimer,  and  bi- allelic  inactivating  mutations  of  either  gene  are 
predicted  to  result  in  MSI.  PCR  of  microsatellite  loci  revealed  MSI 
in  all  hypermutated  tumors,  from  both  PDX  and  autopsy  patients 
(Fig.  lb;  Supplementary  Data  1).  Phenotypic  MSI  was  also 
detected  directly  from  targeted  next-generation  data  for  all 
hypermutated  tumors,  and  not  detected  in  any  non-hypermu¬ 
tated  tumors  (Supplementary  Data  1;  Supplementary  Fig.  10). 
Immunohistochemistry  (IHC)  for  DNA  MMR  proteins  in 
hypermutated  tumors  demonstrated  complete  loss  of  MSH2  and / 
or  MSH6  in  a  pattern  consistent  with  the  inactivating  mutations 
detected  by  sequencing  (Fig.  lc;  Supplementary  Fig.  11).  Non- 


hypermutated  tumors  were  microsatellite  stable  (Tables  1  and  2; 
Supplementary  Data  1)  and  had  intact  MSH2  and  MSH6  pro¬ 
teins,  except  LuCaP  145,  which  exhibited  heterogeneous  loss  of 
MSH6  protein  (Fig.  lc).  MLH1  methylation  was  not  detected  in 
any  of  the  MSI  positive  tumors  (Supplementary  Fig.  12),  and 
MLH1  protein  expression  was  intact  by  IHC  in  MSI-positive 
tumors  except  in  06-134  that  had  homozygous  MLH1  gene 
deletion  (Supplementary  Fig.  13),  arguing  that  MLH1  epigenetic 
silencing  was  not  responsible  for  MSI  in  any  of  the  tumors  in 
our  series. 

Discussion 

Our  findings  support  the  conclusion  that  the  hypermutated 
subtype  of  prostate  cancer  is  chiefly  due  to  loss-of-function 
mutations  in  MSH2  and  MSH6  that  result  in  MSI.  Mutations 
were  predicted  to  be  bi-allelic  in  all  cases  except  00-010,  which 
may  harbour  a  second  undetected  mutation.  Most  interestingly, 
four  of  seven  hypermutated  cases  had  complex  structural 
rearrangements  in  MSH2  and  MSH6  that  were  not  detected  by 
exome  sequencing  in  the  same  samples,  and  would  also  not  be 
expected  to  be  detected  by  traditional  exon-based  Sanger 
sequencing  methods.  Several  previous  studies  have  reported 
MMR  protein  loss  and  MSI  in  both  primary  and  advanced 
prostate  cancers,  but  very  few  MMR  mutations  have  been 
identified8-15.  We  speculate  that  technical  limitations  have  led  to 
an  underestimation  of  MMR  gene  mutations  in  prostate  cancer. 

Our  finding  of  predominantly  MSH2  and  MSH6  mutations  is 
in  contrast  to  colon  and  endometrial  cancer,  where  MSI  is  most 
often  due  to  MLH1  epigenetic  silencing3,4.  This  supports  an 
alternate  mechanism  by  which  MSI  is  acquired  in  prostate  cancer. 
A  recent  study  demonstrated  that  DNA  translocations  and 
deletions  in  advanced  prostate  cancer  occur  in  a  highly 
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Figure  1  |  MSH2  and  MSH6  rearrangements  are  associated  with  loss  of  protein  expression  and  MSI.  (a)  Four  of  seven  hypermutated  cases  had  complex 
rearrangements  in  MSH2  and  MSH6  or  both  genes.  Shown  is  a  representative  complex  MSH2  rearrangement  present  in  hypermutated  cases  LuCaP  147 
and  05-165  (LuCaP  147  was  derived  from  autopsy  patient  05-165).  Breakpoints  were  confirmed  by  Sanger  sequencing.  Genomic  coordinates  are 
hg19.  Detail  on  additional  structural  rearrangements  and  other  mismatch  repair  gene  mutations  is  provided  in  Tables  1  and  2  and  Supplementary  Figs  2-9. 
(b)  Hypermutated  tumors  exhibited  microsatellite  instability  by  PCR.  Shown  is  representative  data  for  LuCaP  58,  which  is  positive  for  MSI  in  3/5 
mononucleotide  marker  systems  (MONO-27,  BAT-25  and  NR-24,  arrows).  All  hypermutated  tumors  tested  were  MSI-PCR  positive  in  at  least  2/5  loci 
(Supplementary  Data  1).  (c)  Hypermutated  tumors  LuCaP  58,  73  and  147  have  loss  of  MSH2  and  MSH6  proteins  by  IHC.  Similar  results  were  observed 
in  hypermutated  tumors  from  rapid  autopsy  patients  (Supplementary  Fig.  11).  A  representative  non-hypermutated  tumour  (LuCaP  23.1)  has  intact 
expression.  LuCaP  145  had  mono-allelic  mutations  in  MSH2  and  MSH6  but  was  not  hypermutated.  IHC  shows  loss  of  MSH6  protein  expression  in 
some  tumour  cells.  Scale  bars,  0.1  mm. 


interdependent  manner,  a  process  termed  ‘chromoplexy’16.  This 
process  may  play  a  role  in  the  genesis  of  MSH2  and  MSH6 
structural  rearrangements  and  deserves  future  study.  Androgen 
receptor  (AR)  function  may  also  play  a  role  in  the  formation  of 
MSH2  and  MSH6  structural  alterations.  AR  has  recently  been 
implicated  in  the  genesis  of  gene  rearrangements  in  prostate 
cancer  by  facilitating  double-strand  DNA  breaks  and  inducing 
non-homologous  end-joining  (reviewed  in  refs  17,18). 

In  summary,  we  have  shown  that  complex  structural  rearrange¬ 
ments  in  mismatch  DNA  repair  genes  MSH2  and  MSH6  are  a 
major  mechanism  underlying  hypermutation  in  advanced  prostate 
cancer.  Future  studies  should  focus  on  determining  if  patients  with 
MMR  gene  defects  exhibit  a  distinct  clinical  course  and  are 
differentially  responsive  to  genotoxic  therapy. 

Methods 

Patients  and  specimens.  The  LuCaP  series  of  prostate  cancer  xenografts  were 
obtained  from  the  University  of  Washington  Prostate  Cancer  Biorepository. 

Human  primary  and  metastatic  prostate  cancer  tissues  were  obtained  as  part  of 
the  University  of  Washington  Prostate  Cancer  Donor  Rapid  Autopsy  Programme. 
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A  haematoxylin  and  eosin  slide  was  reviewed  and  scrolls  from  tissue  blocks  with 
>  50%  estimated  tumour  purity  were  used.  The  Institutional  Review  Board  of  the 
University  of  Washington  approved  all  procedures  involving  human  subjects,  and 
all  subjects  signed  written  informed  consent.  The  sample  size  was  chosen  based  on 
the  number  of  cases  with  suitable  tissues  for  exome  sequencing. 

Genomic  DNA  was  prepared  from  either  formalin-fixed  paraffin-embedded 
tissue  or  from  fresh-frozen  tissue  (for  bone  metastases)  with  the  Gentra  Puregene 
DNA  Isolation  Kit  (Qiagen,  Catalogue  #158489). 

Immunohistochemistry.  Expression  of  MMR  proteins  was  determined  by  IHC 
using  a  tissue  microarray  (UWTMA55),  that  consisted  of  155  metastatic  prostate 
cancer  sites  from  50  patients,  including  77  soft  tissue  metastases  and  83  bone 
metastases),  UWTMA52  consisting  of  primary  prostate  cancer  obtained  at  the  time 
of  radical  prostatectomy  from  127  patients,  and  UWTMA  63  that  consisted  of 
prostate  cancer  tissue  from  32  different  LuCaP  xenograft  lines.  All  the  tissue  cores 
were  duplicated. 

Formalin-fixed  paraffin- embedded  tissue  sections  (5  pm)  were  deparaffinized 
and  rehydrated  with  three  changes  of  xylene  and  graded  ethanol.  Antigen  retrieval 
was  performed  with  heat-induced  epitope  retrieval  for  20  min.  Endogenous 
peroxide  and  avidin/biotin  was  blocked  and  sections  were  then  blocked  with  5% 
normal  goat-horse-chicken  serum  at  room  temperature  for  1  h,  and  incubated  with 
primary  antibody  (listed  in  table  below)  at  4  °C  overnight.  After  washing  three 
times  with  1  x  PBS,  slides  were  incubated  with  biotinylated  secondary  antibody 
(Vector  Laboratories  Inc.),  followed  by  ABC  reagent  (Vector  Laboratories  Inc.) 


NATURE  COMMUNICATIONS  |  5:4988  |  DOI:  10.1038/ncomms5988  |  www.nature.com/naturecommunications 

©  2014  Macmillan  Publishers  Limited.  All  rights  reserved. 


NATURE  COMMUNICATIONS  |  DPI:  10.1038/ncomms5988 


ARTICLE 


and  stable  diaminobenzidine  (Invitrogen  Corp.).  All  sections  were  lightly 
counterstained  with  haematoxylin  and  mounted  with  Cytoseal  XYL  (Richard  Allan 
Scientific).  Mouse  or  rabbit  immunoglobulin-G  was  used  at  the  same  concentration 
as  the  primary  antibody  for  negative  controls.  Antibodies  and  dilutions  used  for 
IHC  are  given  in  Supplementary  Table  4. 

Immunostaining  was  assessed  using  a  quasi-continuous  score  system,  created 
by  multiplying  each  intensity  level  (‘O’  for  no  brown  colour,  T  for  faint  and  fine 
brown  chromogen  deposition  and  ‘2’  for  clear  and  coarse  granular  chromogen 
clumps)  with  the  corresponding  percentage  of  cells  expressing  the  particular 
intensity,  and  then  summing  all  values  to  get  a  final  score  for  each  sample  (scores 
ranging  from  0  to  200).  Only  nuclear  staining  was  evaluated.  Samples  with 
damaged  tissue  core,  missing  tissue  core  or  poor  quality  of  tissue  were  excluded 
from  finial  analysis. 

Microsatellite  instability  PCR.  MSI-PCR  testing  was  performed  by  the  University 
of  Washington  (UW)  clinical  genetics  and  solid  tumors  laboratory  using  the 
Promega  MSI  analysis  kit  (Promega,  Madison,  WI,  USA)  following  the  manu¬ 
facturer’s  instructions.  Specimens  demonstrating  instability  within  two  or  more  of 
the  five  mononucleotide  markers  included  in  this  panel  were  considered  ‘MSI 
positive’,  others  were  considered  ‘MSI  negative’.  The  microsatellite  loci  tested  in  the 
Promega  MSI  analysis  kit  were  NR-21,  BAT-26,  BAT-25,  NR-24  and  MONO-27 
(Genbank  Accession  #  XM_033393,  U41210,  L04143,  X60152,  AC007684, 
respectively). 

MLH1  methylation  analysis.  Two  to  four  hundred  nanograms  of  DNA  from  each 
sample  was  bisulfite  converted  using  the  EZ  DNA  Methylation  Kit  (Zymo 
Research,  Irvine,  CA,  USA)  and  eluted  in  20  |il  volume,  according  to  manu¬ 
facturer’s  protocol. 

SYBR  Green  qPCR  to  detect  methylated  and  unmethylated  MLH1  was 
performed  using  a  CFX  96  Touch  Real-Time  PCR  Detection  System  (Bio-Rad, 
Hercules,  CA,  USA)  with  a  final  reaction  volume  of  20  pi,  consisting  of  500  nM 
each  primer,  9ng  of  bisulfite-converted  genomic  DNA  and  iTaq  Universal  SYBR 
Green  Supermix  at  the  following  conditions:  95  °C  for  3.5  min,  followed  by  40 
cycles  at  95  °C  for  5  s  and  60  °C  for  30  s.  The  unique  primer  sequences  for 
methylated  MLH1  were  5LCGGATAGCGATTTTTAACGC-3'  (forward)  and 
y-CCTAAAACGACTACTACCCG-T  (reverse),  and  for  unmethylated  MLH1 
were  5'  - AATGAATTAATAGGAAGAGTGGATAGT-3'  (forward)  and  5'-TCTCT 
TCATCCCTCCCTAAAACA-3'  (reverse)  (ref.  19).  The  four  primers  each  also 
included  a  20  bp  GC-rich  tail  ( 5'  -  GCGGTCCC  A  A  A  AGGGT  C  AGT  -  3' )  at  their 
5'  end.  Repetitive  Alu  sequence  (‘AluC4’)  was  used  to  normalize  for  the  amount  of 
input  DNA2.  The  absolute  quantitation  of  methylated  and  unmethylated  MLH1  in 
each  sample  was  determined  by  using  the  Epitect  human  methylated  and 
unmethylated  DNA  (Qiagen,  Germantown,  MD,  USA)  to  create  a  standard  curve. 
The  SYBR  Green  assay  results  are  expressed  as  ratios  between  methyl-MLHJ  or 
unmethyl-MLHI  values  and  the  ALUC4  control  values.  The  error  bars  represent 
the  s.e.m. 

Exome  sequencing.  Exome  sequencing  for  autopsy  samples  was  performed  using 
the  Nimblegen  EZ  SeqCap  kit  (Roche)1,20.  Shotgun  libraries  were  constructed  by 
shearing  DNA  and  ligating  sequencing  adaptors.  Libraries  were  hybridized  to  either 
the  EZSeqCap  VI  or  V2  solution-based  probe,  amplified  and  sequenced  on  either 
the  Illumina  GAIIx  or  HiSeq  platform.  For  all  metastases,  somatic  mutations  were 
called  using  Mutect  using  default  parameters  with  matched  normal  (non-tumour) 
samples.  To  remove  common  polymorphisms  and  other  artifacts,  we  imposed  a 
number  of  additional  requirements,  including  requiring  variants  to  be  observed 
with  a  variant  allele  fraction  of  at  least  10%  within  a  tumour,  removing  variants 
present  within  dbSNP  vl37  that  had  first  been  stripped  of  all  disease-associated 
variants  and  removing  variants  that  were  present  at  an  allele  balance  of  40%  or 
more  in  any  germline  sample.  All  exome  sequencing  was  performed  on  fresh- 
frozen  tissue  samples. 

Exome  data  for  PDX  samples  was  from  Kumar  et  al.1,  where  hypermutation 
status  was  previously  characterized  based  on  the  distribution  of  mutations  across 
samples.  For  the  xenografts,  because  corresponding  normal  germline  DNA  was  not 
available,  tumour  sequences  were  compared  against  a  database  of  common 
germline  variants.  The  variants  remaining  were  termed  novel  single-nucleotide 
variants  SNVs  (‘novSNV’)  and  the  estimated  the  contribution  of  germline  variants 
was  ~  200  and  sometimes  more  per  individual.  novSNV  counts  from  Kumar  et  al. 1 
are  provided  in  Supplementary  Table  1. 

Targeted  deep  sequencing  by  BROCA.  Targeted  deep  sequencing  of  DNA  repair 
pathway  genes  was  performed  using  the  BROCA  assay  in  the  UW  clinical  genetics 
and  solid  tumors  laboratory5.  Three  micrograms  of  DNA  was  sonicated  to  a  peak 
of  200  bp  on  a  Covaris  S2  instrument  (Covaris,  Woburn,  MA,  USA).  Following 
sonication,  DNA  was  purified  with  AMPure  XP  beads  (Beckman  Coulter,  Brea  CA, 
USA)  and  subjected  to  three  enzymatic  steps:  end  repair,  A-tailing  and  ligation  to 
Illumina  paired- end  adaptors  as  described  in  the  SureSelectXT  Target  Enrichment 
for  Illumina  multiplexed  sequencing,  which  is  available  for  free  download. 
Adapter-ligated  library  was  PCR  amplified  for  five  cycles  with  Illumina  primers  1.0 
and  2.0  and  individual  paired-end  libraries  (500  ng)  were  hybridized  to  a  custom 


design  of  complementary  RNA  biotinylated  oligonucleotides  targeting  53  genes  in 
52  genomic  regions  (Supplementary  Table  2).  The  120-mer  oligonucleotide  baits 
were  designed  in  Agilent’s  e Array  web  portal  with  the  following  parameters: 
centred  tiling,  3  x  bait  overlap  and  a  maximum  overlap  of  20  bp  into  repetitive 
regions.  The  custom  design  targets  a  total  of  1.4  Mb  of  DNA.  Following  capture, 
each  library  was  PCR  amplified  for  13  cycles  with  primers  containing  a  unique  6  bp 
index.  Equimolar  concentrations  of  96  libraries  were  pooled  to  a  final 
concentration  of  10  pM,  denatured  with  3N  NaOH,  and  cluster  amplified  with  a 
cBot  instrument  on  a  single  lane  of  an  Illumina  v3  flowcell.  Sequencing  was 
performed  with  2  x  101  bp  paired-end  reads  and  a  7  bp  index  read  using  SBS  v3 
chemistry  on  a  HiSeq2500  (Illumina,  SanDiego,  CA,  USA). 

We  used  our  targeted  tumour  sequencing  bioinformatics  pipeline  for  data 
analysis21.  Reads  were  mapped  to  human  reference  genome  (hgl9/GRCh37)  and 
alignment  performed  using  BWA  v0.6.1-rl0419  and  SAMtools  vO. 1.1820.  SNV  and 
indel  calling  was  performed  through  the  GATK  Universal  Genotyper  using  default 
parameters  and  using  VarScan  v2.3.2  and  PINDEL  version  0.2.42.  Structural 
variants  were  identified  using  CREST  vl.O  and  BreakDancer  vl.l.  For  copy  number 
variant  (CNV)  analysis,  copy  number  states  for  individual  probes  were  initially 
called  using  CONTRA  v2.0.32  with  reference  to  a  CNV  control  comprised  of  reads 
from  two  independent  rounds  of  library  preparation  and  sequencing  of  HapMap 
individual  NA12878.  CNV  calls  were  made  at  the  resolution  of  individual  exons 
using  custom  Perl  scripts.  CNV  plots  were  visualized  using  the  R  package  ggplot2. 

Phenotypic  MSI  was  assessed  directly  from  BROCA  next-generation  sequencing 
data  using  mSINGS  (MSI  by  NGS)22.  This  method  evaluated  up  to  146 
mononucleotide  microsatellite  loci  that  are  captured  by  BROCA  in  both  matched 
normal  non-tumour  and  tumour  samples.  For  each  specimen,  microsatellite  loci 
covered  by  a  read  depth  of  <  30  x  were  excluded  as  not  passing  quality  filter.  For 
each  microsatellite  locus  passing  quality  filter,  the  distribution  of  size  lengths  were 
compared  with  a  population  of  normal  controls.  Loci  were  considered  unstable  if 
the  number  of  repeats  is  statistically  greater  than  in  the  control  population.  A 
fraction  of  >  0.20  (20%  unstable  loci)  was  considered  MSI-high  by  mSINGS  based 
on  validation  with  324  tumour  specimens,  in  which  108  cases  had  MSI-PCR  data 
available  as  a  gold  standard22. 


Confirmation  of  MSH2  and  MSH6  structural  rearrangements.  To  validate 
structural  rearrangement  calls,  we  designed  primers  against  regions  flanking 
putative  breakpoints  using  either  PrimerBlast  (http://www.ncbi.nlm.nih.gov/tools/ 
primer-blast/)  or  Primer3  (http://bioinfo.ut.ee/primer3-0.4.0/primer3/input.htm). 
We  used  the  iProof  High-Fidelity  PCR  kit  (Bio-Rad)  to  perform  PCR  under  the 
following  conditions:  98  °C  for  35  s  followed  by  30-40  cycles  of  55-69  °C  for  30  s, 
72  °C  for  30  s  and  72  °C  for  10  min.  Primers  are  listed  in  Supplementary  Table  5. 
We  submitted  resulting  PCR  products  to  Genewiz  for  Sanger  sequencing  and 
aligned  fragments  to  the  human  genome  reference  sequence  (hgl9)  using  BLAT 
from  the  UCSC  Genome  Browser  (http://genome.ucsc.edu/cgi-bin/hgGateway). 

Copy  number  changes  were  confirmed  by  genomic  microarray.  One  microgram 
of  high  molecular  weight  genomic  DNA  from  each  sample  was  labelled  by  random 
priming  using  the  Agilent  Genomic  DNA  Enzymatic  Labelling  Kit  (Cy3-dUTP.)  A 
pool  of  reference  normal  DNA  (Promega)  was  labelled  with  Cy5-dUTP.  Cy3  and 
Cy5  probes  were  combined  and  hybridized  to  Agilent  2  x  400K  SurePrint  G3  CGH 
Microarrays  and  washed  following  the  manufacturer’s  specifications.  Fluorescent 
array  images  were  collected  using  the  Agilent  DNA  microarray  scanner  G2505C 
and  Agilent  Feature  Extraction  software.  Data  analysis  was  performed  with 
Biodiscovery  Nexus  Copy  Number  6.0  software.  The  FASST2  segmentation 
algorithm  and  default  Agilent  settings  for  significance,  gain  and  loss  thresholds, 
with  at  least  six  probes  per  segment  were  used  to  identify  regions  of  CNV  for  each 
sample.  Results  of  copy  number  analysis  by  genomic  microarray  are  given  in 
Supplementary  Fig.  14. 
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Supplementary  Figure  1 :  Somatic  Mutation  Burden  in  Autopsy  Cases  by  Exome 
Sequencing.  Total  number  of  somatic  nonsynonymous  mutations  by  exome  sequencing 
for  rapid  autopsy  cases.  The  threshold  of  300  mutations  used  to  determine  hypermutation 
status  is  shown  with  the  dashed  line.  Median  and  average  mutation  burden  is  given  for  both 
groups.  Cases  that  had  microsatellite  instability  testing  are  shown  with  yellow  stars  (positive) 
and  black  circles  (negative). 


Pritchard  et  al.  2014 


Supplementary  Figure  2 


LuCaP  58 
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Supplementary  Figure  2:  Detail  on  Mismatch  Repair  Gene  Mutations  in  LuCaP  58.  Two 

inactivating  mutations  were  detected  in  MSH6.  The  first  (1)  is  a  deletion  of  exon  8  through  the 
3’UTR  (top).  Copy  number  was  calculated  from  normalized  depth  of  coverage  of  BROCA 
sequencing  data  and  confirmed  by  genomic  microarray  (data  not  shown).  The  blue  bars 
indicate  exons  1  -1 0  (from  left  to  right)  and  black  bars  are  the  standard  deviation  of  the 
measurement  of  Log2  ratio.  The  second  (2)  is  a  2bp  deletion  resulting  in  a  frameshift  and 
premature  truncation  of  the  MSH6  protein  (c.3799_3800del,  bottom).  Shown  is  a  screenshot 
from  the  integrated  genomics  viewer  of  representative  sequencing  reads.  The  black  bars 
indicated  the  deleted  bases.  The  frameshift  was  detected  in  314  out  of  a  total  of  360 
sequencing  reads,  strongly  supporting  that  there  is  bi-allelic  inactivation  of  MSH6. 
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Supplementary  Figure  3 
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Supplementary  Figure  3:  Detail  on  Mismatch  Repair  Gene  Mutations  in 
LuCaP  73.  Two  large  rearrangement  mutations  were  detected  at  the  MSH2/MSH6 
locus  on  chromosome  2  predicted  to  result  in  bi-allelic  inactivation  of  MSH2.  The 
first  (1)  is  a  3Mb  deletion  that  deletes  both  the  MSH2  and  MSH6  genes.  The 
breakpoints  where  confirmed  by  Sanger  sequencing  as  chr2:45,002,771  - 
48,029,401  in  hg19  genomic  coordinates  (top).  The  second  (2)  is  a  440kb 
inversion  mutation  between  MSH2  intron  8  and  FBXOII  intron  1  that  splits  the 
MSH2  gene  and  is  predicted  to  result  in  loss  of  function.  The  breakpoints  of  the 
inversion  were  confirmed  by  Sanger  sequencing  as  chr2:47687644-chr2:48131365 
(bottom).  There  is  a  short  inserted  sequence  between  the  two  breakpoints. 
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LuCaP  147  and  05-165  (see  also  Figure  1A) 

1)  MSH2-C2orf61  343kb  inversion 

2)  MSH2-KCNK12  74kb  inversion 

3)  MSH2-KCNK1 2  40kb  inversion 
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Supplementary  Figure  4:  Detail  on  Mismatch  Repair  Gene  Mutations  in  LuCaP 
147  and  05-165.  Three  different  inversion  mutations  were  detected  involving  the  MSH2 
gene  that  are  predicted  to  result  in  loss-of-function.  The  inversions  were  detected  in  all 
metastatic  sites  (bone,  adrenal,  liver,  and  lymph  node).  Each  inversion  was  confirmed 
by  Sanger  sequencing  with  breakpoints  given  in  hgl  9  genomic  coordinates.  LuCaP  1 47 
was  derived  from  autopsy  patient  05-165  and  the  same  mutations  were  detected  in 
both,  indicating  that  the  MSH2  structural  rearrangements  are  not  a  result  of 
xenografting. 
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Supplementary  Figure  5:  Detail  on  Mismatch  Repair  Gene  Mutations  in  LuCaP 
145  and  05-144.  Two  mutations  were  detected.  The  first  (1)  is  copy  loss  of  exons  8-16 
of  MSH2  (top).  Copy  number  was  calculated  from  normalized  depth  of  coverage  of 
BROCA  sequencing  data  and  confirmed  by  genomic  microarray  (data  not  shown).  The 
blue  bars  indicate  MSH2  exons  1  -1 6  (from  left  to  right)  and  black  bars  are  the  standard 
deviation  of  the  measurement  of  Log2  ratio.  The  second  (2)  is  a  translocation  between 
MSH6  intron  3  and  TESC  intron  4  on  chromosome  12  q24.22.  The  breakpoints  of  the 
translocation  were  confirmed  by  Sanger  sequencing.  These  tumors  had  neuroendocrine 
differentiation.  Unlike  the  other  tumors  MSH2/MSH6  rearrangements  these  tumors 
were  not  hypermutated  and  did  not  demonstrate  MSI,  most  likely  because  one  copy  of 
MSH2  and  MSH6  remain  functionally  intact.  One  hypothesis  is  that  the  cancer  was 
‘transitioning’  to  a  hypermutated  state  when  the  patient  died,  with  a  second  hit  in  the 
MSH2  or  MSH6  gene  not  yet  acquired.  LuCaP  145  was  derived  from  autopsy  patient 
05-144  and  the  same  mutations  were  detected  in  both,  indicating  that  the  structural 
rearrangements  are  not  a  result  of  xenografting. 
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Supplementary  Figure  6 
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Supplementary  Figure  6:  Detail  on  Mismatch  Repair  Gene  Mutations  in  03-130.  There  was 
evidence  of  bi-allelic  loss-of-f unction  mutations  in  MSH2  and  MSH6.  The  first  mutation  (1)  is  a 
translocation  between  MSH2  intron  8  and  chrl  8  q21 .1  (top  left),  in  which  the  breakpoints  were  confirmed 
by  Sanger  sequencing.  The  second  (2)  is  copy  loss  of  MSH2  (top  right,  homozygous  in  exons  1  -8,  likely 
as  a  result  of  the  MSH2 translocation).  The  third  (3)  is  frameshift  mutation  in  exon  4  of  MSH6  (c.2690del, 
p.N897lfsX9).  The  black  bars  indicated  the  deleted  bases.  The  frameshift  was  detected  in  366  out  of  a 
total  of  547  sequencing  reads  (67%),  despite  admixture  of  tumor  with  normal  cells  in  the  sample  tested, 
supporting  that  there  is  bi-allelic  inactivation  of  MSH6  in  tumor.  The  fourth  (4)  is  copy  loss  of  MSH6, 
probably  single  copy.  In  copy  number  plots  blue  bars  indicate  exons  and  black  bars  are  the  standard 
deviation  of  the  measurement  of  Log2  ratio.  Copy  number  was  calculated  by  normalized  depth  of 
coverage  of  BROCA  sequencing  and  confirmed  by  genomic  microarray  (data  not  shown) 
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Supplementary  Figure  7 
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Supplementary  Figure  7:  Detail  on  Mismatch  Repair  Gene  Mutations  in  06- 
134.  The  sample  tested  had  MLH1  copy  loss  which  is  very  likely  to  be  homozygous 
and  result  in  complete  loss  of  MLH1  protein  function.  Depicted  is  copy  number 
analysis  by  BROCA  targeted  deep  sequencing.  MLH1  deletion  was  confirmed  by 
genomic  microarray  (data  not  shown).  The  Log2  ratio  in  relationship  to  a  normal 
female  control  patient.  BRCC3  (far  right)  on  the  X  chromosome  can  be  used  to 
calibrate  the  Log2  ratio  expected  with  heterozygous  copy  loss  because  this  is  a 
male  patient.  PTEN  is  also  deleted  in  this  patient’s  tumor,  a  common  event  in 
metastatic  prostate  cancer. 


Pritchard  et  al.  2014 


Supplementary  Figure  8 
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Supplementary  Figure  8:  Detail  on  Mismatch  Repair  Gene  Mutations  in  00- 
010.  The  primary  prostate  and  liver  tumor  metastasis  samples  tested  had  a 
frameshift  loss-of-function  mutation  in  MSH2  (c.2364_2365insTACA, 
p.A789YfsX11).  The  purple  “I”  indicates  the  position  of  the  inserted  bases, 
visualized  in  the  integrated  genomics  viewer.  A  second  loss-of-function  mutation 
was  not  detected.  It  is  suspected  a  second  MSH2  loss-of-function  mutation  is 
present  because  the  case  was  MSI  high,  had  loss  of  MSH2  and  MSH6  protein  by 
IHC,  was  MLH1  unmethylated,  and  had  intact  MLH1  protein  IHC  (see  separate 
figures). 
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1 )  MSH2  frameshift  (c.  1 1 24_1 1 25insG) 

2)  MSH2  frameshift  (c.1082del) 

3)  MLH1  frameshift  (c.  1 31  Odel),  lymph  node  only 
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Supplementary  Figure  9:  Detail  on  Mismatch  Repair  Gene  Mutations  in  05-123.  The 

primary  prostate  and  lymph  node  metastasis  tumor  samples  tested  had  bi-allelic  frameshift 
loss-of-function  mutations  in  MSH2ex on  7  (c.1082del  and  c.1124_1125insG,  top).  The  black 
bars  indicate  sequencing  reads  with  the  c.1082del  deletion;  the  purple  “I”  indicates  the 
position  of  the  inserted  bases  in  the  c.1124_1125insertion,  visualized  in  the  integrated 
genomics  viewer.  Note  that  these  two  mutations  do  not  occur  on  sequencing  reads  from  the 
same  allele,  strongly  supporting  that  they  are  in  trans  (bi-allelic).  In  addition,  an  MLH1 
frameshift  mutation  in  exon  12  (c.1310delC)  was  detected  in  the  lymph  node  metastasis 
sample  only  (bottom). 
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Supplementary  Figure  10 
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Supplementary  Figure  10:  Microsatellite  Instability  Results  by  BROCA  Next- 
Generation  Sequencing.  We  developed  an  approach  to  measure  microsatellite  instability 
directly  from  next-generation  sequencing  data  that  we  call  mSINGS.  This  method  is 
described  in  the  methods  section  (Salipante  et  al.  2014  Clinical  Chemistry,  in  press).  The 
fraction  of  unstable  microsatellite  loci  out  of  a  maximum  of  146  mononucleotide  microsatellite 
loci  captured  by  BROCA  is  given  on  the  y-axis.  A  threshold  of  0.2  (20%)  unstable  loci  was 
established  as  a  cutoff  for  microsatellite  instability  (dashed  line).  Solid  lines  represent  the 
median  fraction  of  unstable  loci  for  hypermutated  and  non-hyermutated  cases.  The  raw  data 
used  to  generate  this  summary  figure,  including  which  loci  were  unstable  by  BROCA  in  each 
sample  is  given  in  Supplementary  Data  3. 
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Supplementary  Figure  11 :  IHC  results  for  prostate  rapid  autopsy  metastasis  samples. 

Hypermutated  MSI  positive  autopsy  cases  00-010,  03-130,  05-165,  which  harbored  somatic 
mutations  in  MSH2,  MSH6  or  both  genes  show  complete  loss  of  MSH2  and  MSH6 
expression  by  IHC  using  a  tissue  microarray  (top  panels).  Tissue  was  not  available  for  IHC 
studies  in  hypermutated  case  05-123.  Hypermutated  MSI  positive  case  06-134,  which  had 
somatic  deletion  of  MLH1,  has  focal  intact  nuclear  expression  of  MSH2  and  MSH6  protein 
(top  left,  insets).  By  contrast,  MSH2  and  MSH6  nuclear  staining  is  intact  in  MSI-negative 
autopsy  cases  00-029,  01-095,  05-214,  and  05-144  (bottom  panels,  examples  of  positive 
nuclear  staining  in  insets).  MSH6,  but  not  MSH2  protein  expression  was  detected  MSI- 
negative  case  05-116,  a  case  that  also  had  absent  MLH1  protein  (see  separate  figure).  This 
could  reflect  a  false  negative  result  due  to  poor  quality  tissue  for  this  sample  on  the  tissue 
microarray.  For  MSH2  and  MSH6,  heterogeneity  of  immunostaining  is  common  in  tumor 
tissue,  and  protein  expression  is  generally  considered  intact  if  any  cells  display  positive 
nuclear  staining.  Because  MSH2  and  MSH6  function  as  a  heterodimer,  mutations  in  one 
gene  frequently  result  in  loss  of  expression  of  both  proteins,  particularly  when  there  are 
MSH2  mutations.  All  samples  are  from  metastases  and  not  primary  tumors.  Scale  bar: 
0.1mm. 
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Supplementary  Figure  12:  Hypermutated  Prostate  Tumors  Do  Not  Exhibit  MLH1 
Methylation.  Results  of  a  MLH1  methylation-specific  SYBR  green  assay  are  expressed  as 
ratios  between  Methyl-MLHI  or  Unmethyl-MLHI  values  and  the  ALUC4  control  values. 
Genomic  DNA  samples  were  bisulfite-treated  with  EZ  DNA  Methylation  Kit  (Zymo  Research, 
Irvine,  CA)  according  to  manufacturer’s  protocol.  The  primers  used  in  the  SYBR  Green  Assay 
were  previously  described  (see  methods).  The  error  bars  represent  the  standard  error  of  the 
mean.  DNA  samples  from  2  known  MLH1  methylated  colon  cancer  tumors  (POS  CRC1  and 
POS  CRC2)  and  cancer  cell  lines  RKO  (methylated)  and  SW480  (unmethylated)  are  used  as 
controls  for  the  assay. 
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Supplementary  Figure  13:  Hypermutated  MSI  Positive  Prostate  Tumors  With  MSH2  or 
MSH6  mutations  have  intact  MLH1  protein  by  IHC.  Hypermutated  MSI  positive  cases 
LuCaP  58,  LuCaP  73,  LuCaP  147  (top)  and  autopsy  cases  00-010,  03-130,  05-165  (middle), 
which  harbored  somatic  mutations  in  MSH2,  MSH6  or  both  have  intact  MLH1  expression  by 
IHC  using  a  tissue  microarray.  This  corroborates  the  MLH1  methylation  studies  and  strongly 
argues  against  MLH1  epigenetic  silencing  as  a  mechanism  of  MSI  in  these  tumors. 
Hypermutated  MSI  positive  case  06-134  that  had  homozygous  deletion  of  MLH1  has  absent 
MLH1  protein.  Tissue  was  not  available  for  IHC  studies  in  hypermutated  case  05-123.  MLH1 
protein  expression  was  not  detected  MSI-negative  case  05-116  (bottom),  a  case  that  also  had 
absent  MSH2  protein  (see  separate  figure).  This  could  reflect  a  false  negative  result  due  to 
poor  quality  tissue  for  this  sample  on  the  tissue  microarray.  For  MLH1 ,  heterogeneity  of 
immunostaining  is  common  in  tumor  tissue,  and  protein  expression  is  generally  considered 
intact  if  any  cells  display  positive  nuclear  staining.  Scale  bar:  0.1  mm. 
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Supplementary  Figure  14:  Confirmation  of  MSH2,  MSH6,  and  MLH1  Copy  Number 
Status  by  Genomic  Microarray.  Genomic  microarray  (array  CGH)  was  performed  for 
all  cases  that  had  MSH2  and  MSH6  structural  rearrangements  and  also  for  case  06-134 
with  MLH1  homozygous  gene  deletion.  The  genomic  loci  are  given  along  the  top  of 
each  panel.  Y  axes  are  log2  ratio  compared  to  normal  control.  Red  indicates  deletion. 
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SUMMARY 

Toward  development  of  a  precision  medicine 
framework  for  metastatic,  castration-resistant  pros¬ 
tate  cancer  (mCRPC),  we  established  a  multi-institu¬ 
tional  clinical  sequencing  infrastructure  to  conduct 
prospective  whole-exome  and  transcriptome  se¬ 
quencing  of  bone  or  soft  tissue  tumor  biopsies 
from  a  cohort  of  150  mCRPC  affected  individuals. 
Aberrations  of  AR,  ETS  genes,  TP53,  and  PTEN 
were  frequent  (40%-60%  of  cases),  with  TP53  and 
AR  alterations  enriched  in  mCRPC  compared  to 
primary  prostate  cancer.  We  identified  new  genomic 
alterations  in  PIK3CA/B,  R-spondin,  BRAF/RAF1 , 
APC,  |3-catenin,  and  ZBTB 16/ PLZF.  Moreover,  aber¬ 
rations  of  BRCA2,  BRCA1,  and  ATM  were  observed 
at  substantially  higher  frequencies  (19.3%  overall) 
compared  to  those  in  primary  prostate  cancers. 
89%  of  affected  individuals  harbored  a  clinically 
actionable  aberration,  including  62.7%  with  aberra¬ 
tions  in  AR,  65%  in  other  cancer-related  genes,  and 
8%  with  actionable  pathogenic  germline  alterations. 
This  cohort  study  provides  clinically  actionable  infor¬ 
mation  that  could  impact  treatment  decisions  for 
these  affected  individuals. 

INTRODUCTION 

Prostate  cancer  is  among  the  most  common  adult  malig¬ 
nancies,  with  an  estimated  220,000  American  men  diagnosed 
yearly  (American  Cancer  Society,  2015).  Some  men  will  develop 
metastatic  prostate  cancer  and  receive  primary  androgen 
deprivation  therapy  (ADT).  However,  nearly  all  men  with  meta¬ 
static  prostate  cancer  develop  resistance  to  primary  ADT,  a 
state  known  as  metastatic  castration-resistant  prostate  cancer 
(mCRPC).  Multiple  “second  generation”  ADT  treatments,  like 
abiraterone  acetate  (de  Bono  et  al.,  2011;  Ryan  et  al.,  2013) 
and  enzalutamide  (Beer  et  al.,  2014;  Scher  et  al.,  2012),  have 
emerged  for  mCRPC  affected  individuals;  however,  nearly  all 
affected  individuals  will  also  develop  resistance  to  these  agents. 
In  the  U.S.,  an  estimated  30,000  men  die  of  prostate  cancer 
yearly. 

Multiple  studies  have  identified  recurrent  somatic  mutations, 
copy  number  alterations,  and  oncogenic  structural  DNA 
rearrangements  (chromoplexy)  in  primary  prostate  cancer 

(Baca  et  al.,  2013;  Barbieri  et  al.,  2012;  Berger  et  al.,  2011; 


Cooper  et  al.,  2015;  Pflueger  et  al.,  2011;  Taylor  et  al.,  2010; 

Tomlins  et  al.,  2007;  Wang  et  al.,  2011).  These  include  point 
mutations  in  SPOP,  FOXA1 ,  and  TP53 ;  copy  number  alterations 
involving  MYC,  RBI,  PTEN,  and  CHD1;  and  E26  transforma¬ 
tion-specific  (ETS)  fusions,  among  other  biologically  relevant 
genes.  Although  certain  primary  prostate  cancer  alterations 
or  signatures  have  prognostic  clinical  significance  (Hieronymus 
et  al.,  2014;  Lalonde  et  al.,  2014),  the  therapeutic  impact 
of  primary  prostate  cancer  genomic  events  has  not  yet  been 
realized. 

Genomic  studies  of  metastatic  prostate  cancers  demon¬ 
strated  additional  alterations  in  AR  (Taplin  et  al.,  1995)  and  in 
the  androgen  signaling  pathway  (Beltran  et  al.,  2013;  Grasso 
et  al.,  2012;  Gundem  et  al.,  2015;  Hong  et  al.,  2015),  although 
these  studies  were  performed  predominantly  using  autopsy 
samples  or  preclinical  models  with  limited  cohort  sizes.  Prospec¬ 
tive  genomic  characterization  of  fresh  biopsy  samples  from  living 
mCRPC  affected  individuals  has  been  limited  due  to  challenges 
in  obtaining  adequate  tumor  tissue,  especially  from  bone  bi¬ 
opsies  (Mehra  et  al.,  201 1 ;  Van  Allen  et  al.,  2014a),  which  is  the 
most  common  site  of  metastatic  disease.  Thus,  the  landscape 
of  genomic  alterations  in  mCRPC  disease  remains  incompletely 
characterized.  Moreover,  the  low  frequency  of  actionable 
genomic  alterations  in  primary  prostate  cancer  has  limited  the  in¬ 
clusion  of  mCRPC  among  cohorts  wherein  precision  cancer 
medicine  approaches  have  been  piloted  to  guide  treatment  or 
clinical  trial  enrollment. 

We  conducted  a  systematic  and  multi-institutional  study 
of  mCRPC  tumors  obtained  from  living  affected  individuals 
to  determine  the  landscape  of  somatic  genomic  alterations 
in  this  cohort,  dissect  genomic  differences  between  primary 
prostate  cancer  and  mCRPC,  and  discover  the  potential 
relevance  of  these  findings  from  a  biological  and  clinical 
perspective. 

RESULTS 

Clinical,  Biopsy,  and  Pathology  Parameters 

An  international  consortium  consisting  of  eight  academic  medi¬ 
cal  center  clinical  sites  was  established  to  capture  fresh  clinical 
mCRPC  affected  individual  samples  as  part  of  standard-of-care 
approaches  or  through  a  cohort  of  prospective  clinical  trials  (Fig¬ 
ures  1 A  and  IB).  Standard-of-care  approaches  for  mCRPC 
included  abiraterone  acetate  or  enzalutamide.  Clinical  trials 
included  in  this  study  focused  on  combination  strategies 
involving  abiraterone  acetate  or  enzalutamide,  inhibitors  of 
poly  ADP  ribose  polymerase  (PARP),  or  inhibitors  of  aurora  ki¬ 
nase.  Here,  we  report  the  results  of  genomic  profiling  from 
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Figure  1.  Overview  of  the  SU2C-PCF  IDT  Multi-Institutional  Clinical  Sequencing  of  the  mCRPC  Project 

(A)  Schema  of  multi-institutional  clinical  sequencing  project  work  flow. 

(B)  Clinical  trials  associated  with  the  SU2C-PCF  mCRPC  project. 

(C)  Biopsy  sites  of  the  samples  used  for  clinical  sequencing. 

(D)  Histopathology  of  the  cohort.  Representative  images  of  morphological  analysis  of  mCRPC  are  shown  along  with  prevalence  in  our  cohort. 


mCRPC  biopsy  samples  obtained  at  time  of  entry  into  the 
cohort  study.  Future  reports  will  include  longitudinal  clinical 
data  such  as  treatment  response.  The  consortium  utilized  two 
sequencing  and  analysis  centers,  one  centralized  digital 
pathology  review  center,  and  one  centralized  data  visualiza¬ 
tion  portal  (Cerami  et  al.,  2012;  Gao  et  al.,  2013;  Robinson 
et  al.,  2011;  Thorvaldsdottir  et  al.,  2013).  Cross-validation  of 
sequencing  data  from  the  two  original  sequencing  sites  demon¬ 
strated  comparable  variant  calls  for  adequately  powered  genetic 
loci  (E.M.V.A.,  D.R.,  C.  Morrissey,  C.C.P.,  S.L.  Carter,  M.  Rosen¬ 
berg,  A.  McKenna,  A.M.C.,  L.A.G.,  and  P.S.N.,  unpublished 
data). 

Here,  we  describe  150  affected  individuals  with  metastatic 
disease  with  complete  integrative  clinical  sequencing  results 
(whole-exome,  matched  germline,  and  transcriptome  data)  (Fig¬ 
ure  1C)  and  summarized  in  Table  SI.  189  affected  individuals 
were  enrolled  in  the  study,  and  1 75  cases  were  sequenced  after 
pathology  review  and  assessment  of  tumor  content.  Of  these, 
150  biopsies  had  >20%  tumor  content  as  defined  by  computa¬ 


tional  analysis,  based  on  mutant  allele  variant  fractions  and 
zygosity  shifts.  The  biopsies  sequenced  were  from  lymph  node 
(42%),  bone  (28.7%),  liver  (12.7%),  and  other  soft  tissues 
(16.7%).  Baseline  clinical  information  is  available  in  Table  S2.  A 
majority  of  cases  (96.4%)  displayed  typical  high-grade  prostate 
adenocarcinoma  features,  whereas  2.9%  of  cases  showed 
neuroendocrine  differentiation.  One  case  (0.7%)  exhibited 
small-cell  neuroendocrine  features  (Epstein  et  al.,  2014) 
(Figure  ID). 

Landscape  of  mCRPC  Alterations 

Somatic  aberrations  in  a  panel  of  38  statistically  or  clinically  sig¬ 
nificant  genes  are  illustrated  in  Figure  2.  Mean  target  coverage 
for  tumor  exomes  was  160x  and  for  matched  normal  exomes 
was  lOOx.  Although  the  average  mutation  rate  for  mCRPC 
was  4.4  mutations/Mb,  there  were  four  cases  that  exhibited  a 
mutation  rate  of  nearly  50  per  Mb,  three  of  which  are  likely  due 
to  alterations  in  the  mismatch  repair  genes  MLH1  and  MSH2, 
as  discussed  later. 
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Figure  2.  Integrative  Landscape  Analysis  of  Somatic  and  Germline  Aberrations  in  Metastatic  CRPC  Obtained  through  DNA  and  RNA 
Sequencing  of  Clinically  Obtained  Biopsies 

Columns  represent  individual  affected  individuals,  and  rows  represent  specific  genes  grouped  in  pathways.  Mutations  per  Mb  are  shown  in  the  upper  histogram, 
and  incidence  of  aberrations  in  the  cohort  is  in  the  right  histogram.  Copy  number  variations  (CNVs)  common  to  mCRPC  are  shown  in  in  the  lower  matrix,  with  pink 
representing  gain  and  light  blue  representing  loss.  Color  legend  of  the  aberrations  represented  including  amplification,  two  copy  loss,  one  copy  loss,  copy  neutral 
loss  of  heterozygosity  (LOH),  splice  site  mutation,  frameshift  mutation,  missense  mutation,  in-frame  indel,  and  gene  fusion.  Cases  with  more  aberration  in  a  gene 
are  represented  by  split  colors. 


Frequent  copy  number  gains  of  8q,  as  well  as  copy  number 
losses  of  8p,  13q,  16q,  and  18q,  were  also  observed.  The 
mean  number  of  identified  biologically  relevant  genetic  aberra¬ 
tions  per  case  was  7.8  (Figure  2).  All  mutations  identified  are  pre¬ 
sented  in  Table  S3.  The  landscape  of  copy  number  alterations 
demonstrated  expected  recurrent  amplification  peaks  (frequent 
AR,  8q  gain)  and  deletion  peaks  (( CHD1 ,  PTEN,  RBI,  TP53)  (Fig¬ 
ure  3A).  Additional  frequent  focal  amplifications  were  observed 
in  regions  encompassing  CCND1  and  PIK3CA  and  PIK3CB.  A 
new  recurrent  focal  homozygous  deletion  event  was  observed 
in  chr11q23,  encompassing  the  transcriptional  repressor 
ZBTB16. 

To  identify  gene  fusions,  analysis  of  215  transcriptome 
libraries  derived  from  the  150  tumor  RNAs  was  performed  and 
identified  4,122  chimeras  with  at  least  4  reads  spanning  the 


fusion  junction.  These  fusion  junctions  resulted  from  2,247 
gene  pairs,  an  average  of  15  gene  fusions  per  tumor  (Table 
S4).  Among  chimeric  fusion  transcripts  identified,  recurrent 
ETS  fusions  (Tomlins  et  al.,  2005)  were  observed  in  84  cases 
(56%),  of  which  the  majority  were  fused  to  ERG  and  others 
to  FLI1,  ETV4,  and  ETV5  (Figure  3B).  In  addition,  potential  clini¬ 
cally  actionable  fusions  (involving  BRAF,  RAF1 ,  PIK3CA/B,  or 
RSP02)  were  seen  in  eight  cases  (Figure  SI  and  covered 
subsequently). 

To  place  the  mCRPC  mutation  landscape  in  the  context  of  pri¬ 
mary  prostate  cancer  somatic  genomics,  we  performed  a  selec¬ 
tive  enrichment  analysis  to  compare  somatic  point  mutations 
and  short  insertion/deletions  observed  in  this  cohort  with  those 
observed  in  somatic  whole-exome  mutation  data  from  440  pri¬ 
mary  prostate  cancer  exomes  (Barbieri  et  al.,  2012;  The  Cancer 
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Figure  3.  Classes  of  Genomic  Aberrations  Seen  in  mCPRC 

(A)  Copy  number  landscape  of  the  SU2C-PCF  mCRPC  cohort.  Individual  chromosomes  are  represented  by  alternating  colors,  and  key  aberrant  genes  are 
indicated. 

(B)  The  gene  fusion  landscape  of  mCRPC.  Pie  chart  of  all  driver  fusions  identified  and  the  box  plot  represents  specific  ETS  fusions. 

(C)  Mutations  enriched  in  mCRPC  relative  to  hormone  naive  primary  prostate  cancer.  Primary  prostate  cancer  data  derived  from  published  studies  (Barbieri  et  al., 
201 2;  The  Cancer  Genome  Atlas,  201 5).  Level  of  CRPC  enrichment  is  represented  by  the  x  axis,  and  MutSig  CRPC  significance  analysis  is  provided  by  the  y  axis. 
Diameters  are  proportional  to  the  number  of  cases  with  the  specific  aberration.  Genes  of  interest  are  highlighted. 

(D)  Classes  of  driver  aberrations  identified  in  mCRPC. 

(E)  Classes  of  clinically  actionable  mutations  identified  in  mCRPC. 


Genome  Atlas,  2015)  (Figure  3C  and  Table  S5).  Focusing  on 
genes  previously  implicated  in  cancer  (n  =  550),  somatic  TP53 
mutations  were  the  most  selectively  mutated  (q  <  0.001 ;  Benja- 
mini-Flochberg),  followed  by  AR,  KMT2D,  APC,  BRCA2,  and 
GNAS  (q  <  0.1;  Benjamini-Hochberg;  Table  S6).  Both  AR  and 
GNAS  were  mutated  exclusively  in  mCRPC.  We  found  no  genes 
selectively  mutated  in  primary  prostate  cancer  compared  to 
mCRPC. 

We  identified  an  established  biological  “driver”  aberration 
in  a  cancer-related  gene  (i.e.,  known  oncogene  or  tumor  sup¬ 
pressor;  Table  S7)  in  nearly  all  the  cases  (Figure  3D).  Although 
99%  of  the  mCPRC  cases  harbored  a  potential  driver  single¬ 
nucleotide  variant  (SNV)  or  indel,  other  classes  of  driver  aberra¬ 
tions  were  also  highly  prevalent.  These  include  driver  gene 
fusions  in  60%,  driver  homozygous  deletions  in  50%  and 
driver  amplifications  in  54%.  Although  informative  mutations 
were  present  in  virtually  all  mCRPC  cases,  63%  harbored 
aberrations  in  AR,  an  expected  finding  in  castrate-resistant 


disease  but  with  higher  frequency  than  in  prior  reports  (Fig¬ 
ure  3E).  Interestingly,  even  when  AR  was  not  considered, 
65%  of  cases  harbored  a  putatively  clinically  actionable  alter¬ 
ation  (defined  as  predicting  response  or  resistance  to  a  ther¬ 
apy,  having  diagnostic  or  prognostic  utility  across  tumor  types) 
(Table  S8)  (Roychowdhury  et  al.,  2011;  Van  Allen  et  al.,  2014c). 
Non-AR  related  clinically  actionable  alterations  included  aber¬ 
rations  in  the  PI3K  pathway  (49%),  DNA  repair  pathway 
(19%),  RAF  kinases  (3%),  CDK  inhibitors  (7%),  and  the  WNT 
pathway  (5%).  In  addition  to  somatic  alterations,  clinically 
actionable  pathogenic  germline  variants  were  seen  in  8%  of 
mCRPC  affected  individuals,  potentially  emphasizing  the  need 
for  genetic  counseling  in  affected  individuals  with  prostate 
cancer. 

Genomically  Aberrant  Pathways  in  mCRPC 

Integrative  analysis  using  both  biological  and  statistical  frame¬ 
works  (Lawrence  et  al.,  2013,  2014)  of  somatic  point  mutations, 
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Figure  4.  Aberrations  in  the  AR  Pathway  Found  in  mCRPC 

(A)  Cases  with  aberrations  in  the  AR  pathway.  Case  numbering  as  in  Figure  2. 

(B)  Key  genes  found  altered  in  the  AR  pathway  of  mCRPC.  DHT,  dihydrotestosterone. 

(C)  Point  mutations  identified  in  AR.  Amino  acids  altered  are  indicated.  NTAD,  N-terminal  activation.  DBD,  DNA-binding.  LBD,  ligand  binding. 

(D)  Splicing  landscape  of  AR  in  mCRPC.  Specific  splice  variants  are  indicated  by  exon  boundaries,  and  junction  read  level  is  provided.  SU2C,  this  mCRPC  cohort. 
PRAD  tumor,  primary  prostate  cancer  from  the  TCGA.  PRAD  normal,  benign  prostate  from  the  TCGA. 

(E)  Homozygous  deletion  oiZBTB16.  Copy  number  plots  with  x  axis  representing  chromosomal  location  and  the  y  axis  referring  to  copy  number  level.  Red  outline 
indicates  region  of  ZBTB16  homozygous  loss. 


short  insertion/deletions,  copy  number  alterations,  fusion  tran¬ 
scripts,  and  focused  germline  variant  analysis  identified  discrete 
molecular  subtypes  of  mCRPC  (Figure  2).  These  subtypes  were 
classified  based  on  alteration  clustering  and  existing  biological 
pathway  knowledge  and  implicated  the  AR  signaling  pathway, 
phosphatidylinositol-4,5-bisphosphate  3-kinase  (PI3K),  WNT, 
DNA  repair,  cell  cycle,  and  chromatin  modifier  gene  sets,  among 
others.  The  most  frequently  aberrant  genes  in  mCRPC  included 
AR  (62.7%),  ETS  family  (56.7%),  TP53  (53.3%),  and  PTEN 
(40.7%)  (Figure  2). 

AR  Signaling  Pathway 

In  aggregate,  107/150  (71.3%)  of  cases  harbored  AR  pathway 
aberrations,  the  majority  of  which  were  direct  alterations 
affecting  AR  through  amplification  and  mutation  (Figure  4A). 


Figure  4B  summarizes  the  key  genes  altered  in  AR  signaling, 
including  AR  itself,  FOXA1  as  a  pioneer  transcription  factor, 
NCOR1/2  as  negative  regulators  of  AR,  SPOP  as  a  putative 
androgen  receptor  transcriptional  regulator  (Geng  et  al.,  2013), 
and  ZBTB16  as  an  AR  inducible  target  gene  that  may  also  nega¬ 
tively  regulate  AR.  Recurrent  hotspot  mutations  in  AR  were 
observed  at  residues  previously  reported  to  confer  agonism  to 
AR  antagonists  such  as  flutamide  (T878A)  and  bicalutamide 
(W742C),  as  well  as  to  glucocorticoids  (L702FI).  Some,  but  not 
all,  of  these  affected  individuals  had  documented  prior  expo¬ 
sures  that  could  explain  enrichment  for  these  mutations.  Addi¬ 
tional  clinical  data  collection  is  ongoing  (Figure  4C).  Rare  AR 
mutations  not  previously  described  were  seen  in  our  cohort, 
although  these  are  of  unclear  functional  significance.  Further¬ 
more,  one  affected  individual  (Case  89)  harbored  two  putatively 
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functional  AR  mutations  (T878A  and  Q903H),  which  may  further 
suggest  intra-tumor  heterogeneity  emerging  in  the  CRPC 
setting  (Carreira  et  al.,  2014).  Analysis  of  AR  splice  variants 
from  RNA-seq  data  demonstrated  a  distribution  of  splice  vari¬ 
ants  observed  throughout  these  mCRPC  tumor  cases  (Fig¬ 
ure  4D).  Analysis  of  the  TCGA  prostate  dataset  revealed  that 
many  of  these  variants  were  also  present  at  varying  levels  in  pri¬ 
mary  prostate  cancer  and  benign  prostate  tissue.  AR-V7,  which 
has  been  implicated  in  abiraterone  acetate  and  enzalutamide 
resistance  (Antonarakis  et  al.,  2014),  was  observed  in  a  majority 
of  pre-abiraterone/enzalutamide  cases  but  at  very  low  ratios 
relative  to  full  length  AR.  Implications  for  treatment  response 
are  unknown  at  this  time. 

In  addition  to  AR  mutations  itself,  we  observed  alterations  in 
AR  pathway  members  (Figure  4A).  These  included  known  alter¬ 
ations  in  NCOR1 ,  NCOR2,  and  FOXA1  that  have  been  previously 
reported  in  primary  prostate  cancers  and  mCRPC  (Barbieri  et  al., 
2012;  Grasso  et  al.,  2012).  In  this  cohort,  truncating  and 
missense  mutations  in  FOXA1  form  a  cluster  near  the  end  of 
the  Forkhead  DNA  binding  domain  (Figure  S2). 

Recurrent  homozygous  deletions  of  the  androgen-regulated 
gene  ZBTB16  (also  known  as  PLZF)  were  seen  in  8  (5%)  cases 
(Figure  4E)  not  previously  reported  in  clinical  mCRPC  biopsies. 
Analysis  of  the  minimally  deleted  region  seen  in  this  cohort  nar¬ 
rowed  the  candidate  genes  in  the  chr11q23  region  to  ZBTB16 
(Figure  S3).  ZBTB16  has  been  previously  implicated  in  prostate 
cancer  tumorigenesis  and  androgen  resistance  in  preclinical 
models  (Cao  et  al.,  2013;  Kikugawa  et  al.,  2006),  with  loss  of 
ZBTB16  upregulating  the  MAPK  signaling  pathway  (Hsieh 
et  al.,  2015). 

New  PI3K  Pathway  Discoveries 

The  PI3K  pathway  was  also  commonly  altered,  with  somatic  al¬ 
terations  in  73/150  (49%)  of  mCRPC  affected  individuals  (Fig¬ 
ure  5A).  This  included  biallelic  loss  of  PTEN,  as  well  as  hotspot 
mutations,  amplifications  and  activating  fusions  in  PIK3CA, 
and  p.E17K  activating  mutations  in  AKT1  (Figure  S2).  Of  note, 
PIK3CA  amplifications  resulted  in  overexpression  compared  to 
the  remaining  cohort  (Figure  S3). 

Interestingly,  mutations  in  another  member  of  the  PI3K  cata¬ 
lytic  subunit,  PIK3CB,  were  observed  in  this  cohort  for  the  first 
time,  at  equivalent  positions  to  canonical  activating  mutations 
in  PIK3CA  (Figure  5B).  PIK3CB  mutations  appeared  in  the 
context  of  PTEN-deficient  cases,  which  is  consistent  with  a  pre¬ 
vious  report  demonstrating  that  some  PTEN-deficient  cancers 
are  dependent  on  PIK3CB,  rather  than  PIK3CA  (Wee  et  al., 
2008).  Furthermore,  two  affected  individuals  harbored  fusions 
involving  PIK3CA/B,  with  these  events  resulting  in  overexpres¬ 
sion  of  the  gene  relative  to  other  tumors  in  the  cohort  (Figures 
5C  and  5D). 

New  Writ  Pathway  Discoveries 

27/150  (18%)  of  our  cases  harbored  alterations  in  the  Wnt 
signaling  pathway  (Figure  6A).  Hotspot  activating  mutations  in 
CTNNB1  were  seen  (Figure  6B),  as  previously  described  (Voel- 
ler  et  al.,  1998).  Notably,  recurrent  alterations  in  APC  were  also 
observed,  which  have  not  been  previously  described  in  clinical 
mCRPC  affected  individuals.  This  prompted  a  broader  exami¬ 


nation  of  Wnt  signaling  genes  (Figure  6B).  Through  integrative 
analysis,  we  identified  alterations  in  RNF43  and  ZNRF3,  which 
were  recently  described  in  colorectal,  endometrial,  and  adreno¬ 
cortical  cancers  (Assie  et  al.,  2014;  Giannakis  et  al.,  2014)  and 
were  mutually  exclusive  with  APC  alterations  (Figure  6A).  More¬ 
over,  we  also  discovered  R-spondin  fusions  involving  RSP02, 
as  previously  observed  in  colorectal  carcinoma  (Seshagiri 
et  al.,  2012)  in  association  with  RSP02  overexpression  in  these 
cases  (Figure  6C).  RSP02  is  a  key  factor  in  prostate  cancer  or¬ 
ganoid  methodology  (Gao  et  al.,  2014).  Affected  individuals  s 
with  aberrations  in  RNF43,  ZNRF3,  or  RSP02  (overall  6%  of 
affected  individuals)  are  predicted  to  respond  to  porcupine  in¬ 
hibitors  (Liu  et  al.,  2013). 

Cell-Cycle  Pathway 

We  observed  RBI  loss  in  21  %  of  cases  (Figure  S4).  Expanding 
the  scope  of  cell-cycle  genes  implicated  in  mCRPC,  we  noted 
focal  amplifications  involving  CCND1  in  9%  of  cases,  as  well 
as  less  common  (<  5%)  events  in  CDKN2A/B,  CDKN1B,  and 
CDK4  (Figure  S4).  Cell-cycle  derangement,  such  as  through 
CCND1  amplification  or  CDKN2A/B  loss,  may  result  in  enhanced 
response  to  CDK4  inhibitors  in  other  tumor  types  (Finn  et  al., 
2015),  and  preclinical  mCRPC  models  predict  similar  activity  in 
prostate  cancer  (Comstock  et  al.,  2013). 

DNA  Repair  Pathway 

Integrative  analysis  of  both  the  somatic  and  pathogenic  germline 
alterations  in  BRCA2  identified  19/150  (12.7%)  of  cases  with 
loss  of  BRCA2,  of  which  —90%  exhibited  biallelic  loss  (Fig¬ 
ure  7A).  This  was  commonly  a  result  of  somatic  point  mutation 
and  loss  of  heterozygosity,  as  well  as  homozygous  deletion. 
One  of  the  clinical  trials  in  our  consortium  is  evaluating  poly(- 
ADP-ribose)  polymerase  (PARP)  inhibition  in  unselected 
mCRPC  affected  individuals.  Importantly,  multiple  affected  indi¬ 
viduals  in  this  trial  who  experienced  clinical  benefit  harbored 
biallelic  BRCA2  loss,  providing  further  evidence  of  clinical  ac¬ 
tionability  (Mateo  et  al.,  2014).  Eight  affected  individuals 
(5.3%)  harbored  pathogenic  germline  BRCA2  mutations  (Fig¬ 
ure  7B)  with  a  subsequent  somatic  event  that  resulted  in  biallelic 
loss,  revealing  a  surprisingly  high  frequency  relative  to  primary 
prostate  cancer. 

We  therefore  expanded  the  focus  to  other  DNA  repair/recom¬ 
bination  genes  and  identified  alterations  in  at  least  34/150 
(22.7%)  of  cases.  These  include  recurrent  biallelic  loss  of 
ATM  (Figure  7B),  including  multiple  cases  with  germline  patho¬ 
genic  alterations.  ATM  mutations  were  also  observed  in 
affected  individuals  who  achieved  clinical  responses  to  PARP 
inhibition  (Mateo  et  al.,  2014).  In  addition,  we  noted  events  in 
BRCA1,  CDK12,  FANCA,  RAD51B,  and  RAD51C.  If  aberrations 
of  BRCA2,  BRCA1 ,  and  ATM  all  confer  enhanced  sensitivity  to 
PARP  inhibitors,  29/150  (19.3%)  of  mCRPC  affected  individuals 
would  be  predicted  to  benefit  from  this  therapy.  Interestingly, 
three  out  of  four  mCRPC  tumors  exhibited  hypermutation  and 
harbored  alterations  in  the  mismatch  repair  pathway  genes 
MLFI1  or  MSFI2  (Figures  2  and  7C),  corroborating  a  recent 
report  identifying  structural  alterations  in  MSFI2  and  MSFI6 
mismatch  repair  genes  in  hypermutated  prostate  cancers 
(Pritchard  et  al.,  2014). 
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Figure  5.  Aberrations  in  the  PI(3)K  Pathway  Found  in  mCRPC 

(A)  Cases  with  aberrations  in  the  PIK3  pathway.  Case  numbering  as  in  Figure  2. 

(B)  Point  mutations  identified  in  PIK3CB.  Amino  acids  altered  are  indicated.  Analogous,  recurrent  COSMIC  mutations  in  PIK3CAare  shown  as  expansion  views. 

(C)  Outlier  expression  of  PK3CAin  CRPC  case  harboring  the  TBL1XR1-PIK3CA  gene  fusion.  Structure  of  the  gene  fusion  is  inset.  UTR,  untranslated  region.  CDS, 
coding  sequence. 

(D)  As  in  (C),  except  for  PIK3CB  and  the  ACPP-PIK3CB  gene  fusion. 


DISCUSSION 

To  effectively  implement  precision  cancer  medicine,  prospective 
identification  of  predictive  biomarkers  should  be  performed  with 
information  derived  from  the  most  contemporary  tumor  assess¬ 
ments  that  reflect  the  affected  individual’s  prior  therapies  and 
treatment  opportunities.  In  mCRPC,  precision  cancer  medicine 
activities  have  been  limited  by  difficulties  obtaining  clinical  sam¬ 
ples  from  mCRPC  affected  individuals  and  a  lack  of  comprehen¬ 
sive  genomic  data  for  potentially  actionable  alterations.  By 
demonstrating  the  feasibility  of  prospective  genomics  in  mCRPC 
and  defining  the  mutational  landscape  in  a  focused  metastatic 
clinical  cohort,  this  report  may  inform  multiple  genomically  driven 


clinical  trials  and  biological  investigations  into  key  mediators  of 
mCRPC.  In  nearly  all  of  the  mCRPC  analyzed  in  this  study,  we 
identified  biologically  informative  alterations;  almost  all  harbored 
at  least  one  driver  SNV/indel,  and  approximately  half  harbored  a 
driver  gene  fusion,  amplification,  or  homozygous  deletion. 
Remarkably,  in  nearly  90%  of  mCRPC  affected  individuals,  we 
identified  a  potentially  actionable  somatic  or  germline  event. 

The  high  frequency  of  AR  pathway  alterations  in  this  cohort 
strongly  implies  that  the  vast  majority  of  mCRPC  affected 
individuals  remain  dependent  on  AR  signaling  for  viability.  The 
“second-generation”  AR-directed  therapies  (e.g.,  abiraterone 
acetate  and  enzalutamide)  may  select  for  distinct  phenotypes 
that  may  be  indifferent  to  AR  signaling,  and  prospective 
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Figure  6.  Aberrations  in  the  WNT  Pathway  Found  in  mCRPC 

(A)  Cases  with  aberrations  in  the  WNT  pathway.  Case  numbering  as  in  Figure  2. 

(B)  Aberrations  identified  in  APC  and  CTNNB1 .  Amino  acids  altered  are  indicated.  ARM,  armadillo  repeat.  Phos,  phosphorylation  domain.  TAD,  trans- activating 
domain.  EB1,  end  binding  protein-1  domain.  CC,  coiled  coil. 

(C)  Outlier  expression  of  RSP02  in  CRPC  and  the  GRHL2-RSP02  gene  fusion.  RNA-seq  expression  across  our  CRPC  cohort.  Structure  of  the  gene  fusion  is  inset. 
UTR,  untranslated  region.  CDS,  coding  sequence. 


characterization  of  such  cases  will  be  of  particular  interest.  We 
hypothesize  that  affected  individuals  with  acquired  AR  muta¬ 
tions,  including  new  AR  mutations  discovered  in  this  cohort, 
will  harbor  differential  responses  to  these  second-generation 
ADT  therapies.  As  the  number  of  affected  individuals  in  this 
cohort  with  AR  mutations  increases,  we  will  subsequently  be 
able  to  link  specific  AR  mutations  with  clinical  phenotypes  to 
determine  which  mutations  confer  selective  response  or  resis¬ 
tance  to  subsequent  AR-directed  therapy. 

Moreover,  these  data  identify  multiple  therapeutic  avenues 
warranting  clinical  investigation  in  the  CRPC  population. 
Excluding  AR  aberrations,  65%  of  mCRPC  have  a  potentially 
actionable  aberration  that  may  suggest  an  investigational  drug 
or  approved  therapy.  For  example,  focusing  on  the  PI3K 
pathway,  PIK3CB-specific  inhibitors  may  have  utility  in  affected 
individuals  with  mutation,  amplification,  and/or  fusion  of  this 


gene  (Schwartz  et  al.,  2015);  multiple  affected  individuals  who 
achieved  durable  (>1  year)  responses  to  PIK3CB-specificin  inhi¬ 
bition  harbored  activating  mutation  or  amplification  in  PIK3CB 
(J.S.  de  Bono  et  al.,  201 5, 1 06th  Annual  Meeting  of  the  American 
Association  for  Cancer  Research,  abstract).  RAF  kinase  fusions 
in  3%  of  mCPRC  affected  individuals  would  suggest  the  use  of 
pan-RAF  inhibitors  or  MEK  inhibitors  (Palanisamy  et  al.,  2010). 
In  addition,  the  emergence  of  porcupine  inhibitors  (Liu  et  al., 
2013)  and  R-spondin  antibodies  may  warrant  investigation  in 
mCRPC  tumors  harboring  Wnt  pathway  alterations  or  specif¬ 
ically  R-spondin  fusions,  respectively.  These  observations  will 
need  to  be  prospectively  assessed  in  the  clinical  trials. 

Additionally,  biallelic  inactivation  of  BRCA2,  BRCA1 ,  or  ATM 
was  observed  in  nearly  20%  of  affected  individuals.  Previous 
work  in  other  cancer  types  suggests  that  these  affected  individ¬ 
uals  may  benefit  from  PARP  inhibitors  (Fong  et  al.,  2009; 
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Figure  7.  Aberrations  in  the  DNA  Repair  Pathway  Found  in  mCRPC 

(A)  Cases  with  aberrations  in  the  DNA  repair  pathway.  Case  numbering  as  in  Figure  2. 

(B)  Aberrations  identified  in  BRCA2,  ATM,  and  BRCA1.  Amino  acids  altered  are  indicated.  HELC,  helical  domain.  OB,  oligonucleotide  binding  fold.  FAT,  FRAP- 
ATM-TRRAP  domain.  PIK3c,  PI3  kinase  domain.  CC,  coiled  coil.  BRC,  Brea  repeat. 

(C)  Microsatellite  instability  analysis  of  representative  hypermutated  CRPC  cases  and  non-hypermutated  cases. 


Kaufman  et  al.,  2015;  Weston  et  al.,  2010)  or  platinum-based 
chemotherapy,  and  prior  reports  have  implicated  the  presence 
of  germline  BRCA2  alterations  in  primary  prostate  cancer  with 
poor  survival  outcomes  (Castro  et  al.,  201 3).  Given  the  incidence 
of  pathogenic  germline  BRCA2  mutations  in  this  cohort  with  sub¬ 
sequent  somatic  events  (5%),  along  with  enrichment  for  somatic 
BRCA2  alterations  in  mCRPC  (13%),  germline  genetic  testing  in 
mCRPC  affected  individuals  warrants  clinical  consideration. 

The  ability  to  molecularly  characterize  mCRPC  biopsy  sam¬ 
ples  from  affected  individuals  actively  receiving  therapy  will 
also  enable  focused  studies  of  resistance  to  secondary  ADT 
therapies,  including  neuroendocrine-like  phenotypes.  This  will 
require  iterative  sampling  of  pre-treatment  and  resistant  tumors 
from  matching  affected  individuals  and  may  warrant  multire- 


gional  biopsies  from  affected  individuals  (if  feasible)  given  het¬ 
erogeneity  in  mCRPC  (Carreira  et  al.,  2014;  Gundem  et  al., 

2015).  Toward  that  end,  in  some  affected  individuals,  we 
observed  multiple  AR  mutations  emerging  in  the  same  biopsy, 
which  may  indicate  clonal  heterogeneity  within  these  mCRPC  tu¬ 
mor  samples.  Additional  genomic  alterations  discovered  in  this 
cohort  (e.g.,  ZBTB16)  warrant  exploration  in  prostate  cancer 
model  systems,  including  organoid  cultures  (Gao  et  al.,  2014). 

Broadly,  our  effort  demonstrates  the  utility  of  applying 
comprehensive  genomic  principles  developed  for  primary  malig¬ 
nancies  (e.g.,  TCGA)  to  a  clinically  relevant  metastatic  tumor 
cohort.  Our  effort  may  also  catalyze  multi-institutional  efforts  to 
profile  tumors  from  cohorts  of  affected  individuals  with  metasta¬ 
tic,  treated  tumors  in  other  clinical  contexts  because  our  results 
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demonstrate  multiple  discoveries  within  this  advanced  disease 
stage  that  have  not  been  observed  in  primary  tumor  profiling. 
Moreover,  this  study  sets  the  stage  for  epigenetic  and  other 
profiling  efforts  in  mCRPC  not  taken  in  this  study,  which  may 
enable  biological  discovery  and  have  immediate  therapeutic 
relevance  in  mCRPC  (Asangani  et  al.,  2014).  Overall,  our  efforts 
demonstrate  the  feasibility  of  comprehensive  and  integrative  ge¬ 
nomics  on  prospective  biopsies  from  individual  mCRPC  affected 
individuals  to  enable  precision  cancer  medicine  activities  in  this 
large  affected  individual  population. 

EXPERIMENTAL  PROCEDURES 
Affected  Individual  Enrollment 

Affected  individuals  with  clinical  evidence  of  mCRPC  who  were  being  consid¬ 
ered  for  abiraterone  acetate  or  enzalutamide  as  standard  of  care,  or  as  part  of 
a  clinical  trial,  were  considered  for  enrollment.  Affected  individuals  with  meta¬ 
static  disease  accessible  by  image-guided  biopsy  were  eligible  for  inclusion. 
All  affected  individuals  provided  written  informed  consent  to  obtain  fresh  tu¬ 
mor  biopsies  and  to  perform  comprehensive  molecular  profiling  of  tumor 
and  germline  samples. 

Biopsies  and  Pathology  Review 

Biopsies  of  soft  tissue  or  bone  metastases  were  obtained  under  radiographic 
guidance.  Digital  images  of  biopsy  slides  were  centrally  reviewed  using 
schema  established  to  distinguish  usual  adenocarcinoma  from  neuroendo¬ 
crine  prostate  cancer  (Epstein  et  al.,  2014).  All  images  were  reviewed  by  geni¬ 
tourinary  oncology  pathologists  (M.R.,  J.M.M.,  L.P.K.,  S.A.T.,  R.M.,  V.R.,  A.G., 
M.L.,  R.L.,  and  M.B.). 

Sequencing  and  Analysis 

Normal  DNAs  from  buccal  swabs,  buffy  coats,  or  whole  blood  were  isolated 
using  the  QIAGEN  DNeasy  Blood  &  Tissue  Kit.  Flash-frozen  needle  biopsies 
with  highest  tumor  content  for  each  case,  as  determined  by  pathology  review, 
were  extracted  for  nucleic  acids.  T umor  genomic  DNA  and  total  RNA  were  pu¬ 
rified  from  the  same  sample  using  the  AllPrep  DNA/RNA/miRNA  kit  (QIAGEN) 
with  disruption  on  a  Tissuelyser  II  (QIAGEN).  RNA  integrity  was  verified  on  an 
Agilent  2100  Bioanalyzer  using  RNA  Nano  reagents  (Agilent  Technologies). 

Whole-exome  capture  libraries  were  constructed  from  1 00  ng  to  1  |ig  of  DNA 
from  tumor  and  normal  tissue  after  sample  shearing,  end  repair,  and  phos¬ 
phorylation  and  ligation  to  barcoded  sequencing  adaptors.  Ligated  DNA  was 
size  selected  for  lengths  between  200  and  350  bp  and  subjected  to  hybrid  cap¬ 
ture  using  SureSelect  Exome  v4  baits  (Agilent).  Exome  sequence  data  pro¬ 
cessing  and  analysis  were  performed  using  pipelines  at  the  Broad  Institute 
and  the  University  of  Michigan.  A  BAM  file  aligned  to  the  hgl  9  human  genome 
build  was  produced  using  lllumina  sequencing  reads  for  the  tumor  and  normal 
sample  and  the  Picard  pipeline.  Somatic  mutation  analysis  was  performed  as 
described  previously  (Cibulskis  et  al.,  2013;  Van  Allen  et  al.,  2014c)  and  re¬ 
viewed  with  Integrated  Genomics  Viewer  (IGV)  (Robinson  et  al.,  2011). 

Copy  number  aberrations  were  quantified  and  reported  for  each  gene  as  the 
segmented  normalized  log2-transformed  exon  coverage  ratios  between  each 
tumor  sample  and  matched  normal  sample  (Lonigro  et  al.,  2011).  To  account 
for  observed  associations  between  coverage  ratios  and  GC  content  across 
the  genome,  lowess  normalization  was  used  to  correct  per-exon  coverage  ra¬ 
tios  prior  to  segmentation  analysis.  Mean  GC  percentage  was  computed  for 
each  targeted  region,  and  a  lowess  curve  was  fit  to  the  scatterplot  of  log2- 
coverage  ratios  versus  mean  GC  content  across  the  targeted  exome  using 
the  lowess  function  in  R  (version  2.13.1)  with  smoothing  parameter  f  =  0.05. 
The  resulting  copy  ratios  were  segmented  using  the  circular  binary  segmenta¬ 
tion  algorithm  (Olshen  et  al.,  2004). 

Statistical  analysis  of  recurrently  mutated  genes  was  performed  using  Mut- 
Sig  (Lawrence  et  al.,  2013).  Selective  enrichment  analysis  (Van  Allen  et  al., 
2014b)  of  mutations  observed  in  mCRPC  compared  to  primary  prostate  can¬ 
cer  was  performed  by  tabulating  the  frequency  of  affected-individual-normal¬ 
ized  mutations  observed  in  either  CRPC  or  primary  prostate  cancer  and 


performing  a  two-sided  Fisher’s  exact  test  using  allelic  fraction  cut  off  of  0.1 
or  greater  and  a  set  of  biologically  relevant  cancer  genes  (n  =  550  genes)  (Fu- 
treal  et  al.,  2004).  Multiple  hypothesis  test  correction  was  performed  using 
Benjamini-Hochberg  method. 

Transcriptome  libraries  were  prepared  using  200-1,000  ng  of  total  RNA. 
PolyA+  RNA  isolation,  cDNA  synthesis,  end-repair,  A-base  addition,  and  liga¬ 
tion  of  the  lllumina  indexed  adapters  were  performed  according  to  the  TruSeq 
RNA  protocol  (lllumina).  Libraries  were  size  selected  for  250-300  bp  cDNA 
fragments  on  a  3%  Nusieve  3:1  (Lonza)  gel,  recovered  using  QIAEX  II  reagents 
(QIAGEN),  and  PCR  amplified  using  Phusion  DNA  polymerase  (New  England 
Biolabs).  Total  transcriptome  libraries  were  prepared  as  above,  omitting  the 
poly  A  selection  step  and  captured  using  Agilent  SureSelect  Human  All  Exon 
V4  reagents  and  protocols.  Library  quality  was  measured  on  an  Agilent  2100 
Bioanalyzer  for  product  size  and  concentration.  Paired-end  libraries  were 
sequenced  with  the  lllumina  HiSeq  2500,  (2x100  nucleotide  read  length) 
with  sequence  coverage  to  50  M  paired  reads  and  100  M  total  reads. 

Paired-end  transcriptome  sequencing  reads  were  aligned  to  the  human 
reference  genome  (GRCh37/hg19)  using  a  RNA-seq  spliced  read  mapper 
Tophat2  (Kim  and  Salzberg,  2011)  (Tophat  2.0.4),  with  “-fusion-search”  option 
turned  on  to  detect  potential  gene  fusion  transcripts.  Potential  false-positive 
fusion  candidates  were  filtered  out  using  “Tophat-Post-Fusion”  module. 
Further,  the  fusion  candidates  were  manually  examined  for  annotation  and 
ligation  artifacts.  Gene  expression,  as  fragments  per  kilobase  of  exon  per 
million  fragments  mapped  (FPKM;  normalized  measure  of  gene  expression), 
was  calculated  using  Cufflinks  (Trapnell  et  al.,  2012). 
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Molecular  analysis  of  colon  cancers  currently  requires  multiphasic  testing  that  uses  various  assays  with 
different  performance  characteristics,  adding  cost  and  time  to  patient  care.  We  have  developed  a  single, 
next-generation  sequencing  assay  to  simultaneously  evaluate  colorectal  cancers  for  mutations  in  relevant 
cancer  genes  [KRAS,  NRAS,  and  BRAF)  and  for  tumor  micro  sate  Hite  instability  (MSI).  In  a  sample  set  of  61 
cases,  the  assay  demonstrated  overall  sensitivity  of  100%  and  specificity  of  100%  for  identifying  cancer- 
associated  mutations,  with  a  practical  limit  of  detection  at  2%  mutant  allele  fraction.  MSIplus  was  97% 
sensitive  (34  of  35  MSI-positive  cases)  and  100%  specific  (42  of  42  MSI-negative  cases)  for  ascertaining 
MSI  phenotype  in  a  cohort  of  78  tumor  specimens.  These  performance  characteristics  were  slightly  better 
than  for  conventional  multiplex  PCR  MSI  testing  (97%  sensitivity  and  95%  specificity),  which  is  based  on 
comparison  of  microsatellite  loci  amplified  from  tumor  and  matched  normal  material,  applied  to  the  same 
specimen  cohort.  Because  the  assay  uses  an  amplicon  sequencing  approach,  it  is  rapid  and  appropriate  for 
specimens  with  limited  available  material  or  fragmented  DNA.  This  integrated  testing  strategy  offers 
several  advantages  over  existing  methods,  including  a  lack  of  need  for  matched  normal  material,  sensitive 
and  unbiased  detection  of  variants  in  target  genes,  and  an  automated  analysis  pipeline  enabling  prin¬ 
cipled  and  reproducible  identification  of  cancer-associated  mutations  and  MSI  status  simultaneously. 
(J  Mot  Diagn  2015 ,  17:  1-10;  http://dx.doi.org/10.1016/jfmotdx.2015.05.008) 


After  initial  diagnosis,  the  molecular  characterization  of 
colorectal  cancers  may  require  several  separate  clinical  tests. 
Evaluation  of  microsatellite  instability  (MSI)  status  is  recom¬ 
mended  testing  on  all  primary  colon  cancers  from  patients  50 
years  or  younger  (and  in  older  patients  if  specific  pathological 
features  are  present1)  to  serve  as  a  screening  test  for  Lynch 
syndrome,  a  disease  of  hereditary  cancer  predisposition. 
Moreover,  MSI  status  provides  diagnostic  information  about 
disease  prognosis  and  predicted  treatment  response  to  fluo- 
rouracil,  and  can,  therefore,  directly  inform  patient  care.  In 
the  case  of  metastatic  disease,  additional  molecular  testing 
beyond  MSI  status  is  indicated.  Recently  updated  guidelines 
from  the  National  Comprehensive  Cancer  Network  recom¬ 
mend  that  extended  NRAS  and  KRAS  testing  is  performed  on  all 
stage  IV  cancers  to  identify  mutations  conferring  resistance  to 
epidermal  growth  factor  receptor  inhibitors.  If  RAS  gene 


mutational  status  is  negative,  testing  for  BRAF  mutations  is 
then  also  advised,  given  the  poor  response  of  BRAF- mutated 
tumors  to  cetuximab  plus  irinotecan,  fluorouracil,  and  leuco- 
vorin  combination  therapy,11  and  also  to  epidermal  growth 
factor  receptor  inhibitors  used  beyond  first-line  treatment. 

Independent  molecular  assays  are  currently  used  to  test 
tumor  specimens  for  MSI  status  and  gene  mutations.  Molec¬ 
ular  diagnosis  of  MSI  is  implemented  using  multiplexed  PCR- 
based  MSI  testing  (MSI-PCR),  wherein  a  limited  number  of 
informative  micro  satellite  markers1  are  PCR  amplified 
from  tumor  and  matched  normal  material,  products  are 
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resolved  using  capillary  gel  electrophoresis,  and  the  presence 
of  additional  alleles,  which  are  the  hallmark  of  the  MSI 
phenotype,  is  qualitatively  ascertained.  In  some  cases, 
immunohistochemical  (IHC)  detection  of  mismatch  repair¬ 
pathway  protein  loss  may  additionally  or  alternatively  be 
performed,  although  some  studies  suggest  that  IHC  does 
not  sensitively  identify  all  MSI-positive  tumors.  In 
contrast,  clinical  testing  for  NRAS ,  KRAS ,  and  BRAF  gene 
mutations  is  commonly  achieved  using  single-gene  assays 
using  melting  curve  analysis^  real-time  PCR,  or  con¬ 
ventional  Sanger  sequencing.  Peptide  nucleic  acid  clamp¬ 
ing25  or  selective  amplification  of  mutant  alleles  is  sometimes 
used  to  improve  sensitivity.26-28 

Next-generation  sequencing  (NGS)  is  becoming  increas¬ 
ingly  used  by  clinical  laboratories  as  a  cost-effective  and 
scalable  method  to  interrogate  multiple  genetic  targets  in 
parallel,  and  could  be  adapted  for  the  focused  purpose  of 
characterizing  colorectal  cancers  for  molecular  workup. 
Although  integrating  colorectal  cancer  testing  into  an  NGS 
diagnostic  would  offer  practical  advantages  in  eliminating  the 
need  for  multiple,  separate  diagnostic  tests,  the  unique  analytic 
properties  of  NGS  could  also  translate  to  performance  ad¬ 
vantages  over  existing  testing  methods.  Such  benefits  include 
the  following:  i)  The  ability  of  the  technology  to  detect  low- 
prevalence,  cancer-associated  mutations  is  greater  than  that 
of  conventional  methods,  because  each  DNA  molecule  is 
examined  independently  by  NGS,  potentially  providing 
increased  sensitivity  for  detecting  relevant  cancer-associated 
mutations,  ii)  Our  group  and  others36-  have  developed 
methods  to  computationally  infer  MSI  status  from  NGS  data 
on  the  basis  of  the  quantification  and  distribution  of  observed 
allele  lengths  at  microsatellite  loci.  This  offers  a  standardized, 
statistical  approach  for  interpreting  MSI  testing  results,  in 
contrast  to  the  current  practice  of  subjective  interpretation  of 
MSI-PCR  electropherogram  traces,  iii)  Owing  to  the  digital 
nature  of  NGS  data,  primary  analysis  can  be  readily  auto¬ 
mated,  ensuring  consistency  in  test  interpretation. 

Herein,  we  describe  the  clinical  validation  of  a  novel  NGS 
assay,  MSIplus,  for  simultaneously  evaluating  tumor  MSI 
status  and  mutational  hotspots  in  KRAS ,  NRAS ,  and  BRAF 
genes,  and  provide  interpretive  guidelines  for  its  diagnostic 
use.  In  contrast  to  earlier  research  work,  which  has  used  NGS 
to  evaluate  MSI  status  from  exome  and  targeted  gene  capture 
designs,  our  assay  uses  an  amplicon  sequencing  approach, 
which  enables  effective  targeting  and  high  depth  of  coverage 
for  the  loci  of  interest.  The  assay  is  suitable  for  molecular 
characterization  of  colorectal  cancers,  does  not  require 
matched  normal  material  for  inference  of  MSI  status,  and  is 
rapid,  cost-effective,  sensitive,  and  specific. 

Materials  and  Methods 

Selection  of  Target  Sequences  and  Primer  Design 

For  the  purpose  of  inferring  MSI  status,  we  selected  a  panel 
of  11  microsatellite  markers  from  our  earlier  analysis  of 


colorectal  cancer  exome  data  that  were  empirically  found  to 
be  both  most  discriminatory  for  MSI  and  most  frequently 
unstable  in  MSI-positive  tumors  (Table  1).  We  also  incor-  [Tl] 
porated  the  mononucleotide  A/T  tract  of  HSPH1,  the 
instability  of  which  predicts  sensitivity  to  particular  anti¬ 
cancer  agents,  and  the  five  microsatellite  markers  (MONO- 
27,  BAT-25,  BAT-26,  NR-21,  and  NR-24)  that  compose  a 
performance-enhanced  derivative  of  the  Bethesda  panel 
used  in  current  clinical  MSI-PCR  assays.  Separate  primers 
were  designed  to  span  exons  containing  relevant  mutational 
hotspots  in  KRAS  (exons  2,  3,  and  4;  codons  12,  13,  61,  1 17, 
and  147),  ARABS' (exons  2, 3,  and  4;  codons  12, 13,61, 117,  and 
147),  and  BRAF  (exon  15;  codons  599,  600,  and  601). 

Multiplexed  primer  design  was  performed  using  MPprimer  Q6 
version  1.441  ( https://code.google.eom/p/mpprimer ,  last 

accessed  October  15,  2013)  with  some  manual  curation. 
Genomic  coordinates  for  each  locus  (human  genome  hgl9/ 
GRCh37)  and  PCR  primer  sequences  are  provided  in  Table  1 . 
Primers  were  concatenated  at  the  5'  end  to  partial  Illumina 
sequencing  adaptors,  which  were  extended  in  downstream 
steps.  All  oligonucleotides  were  synthesized  by  Integrated 
DNA  Technologies  (Coralville,  IA). 

Tumor  Specimens  and  Clinical  Testing 

DNA  from  tumor  specimens,  which  were  predominantly 
colorectal  cancers,  but  included  a  small  subset  of  endo¬ 
metrial  cancers,  lung  cancers,  ovarian  cancers,  melanoma, 
and  additional  tumor  types,  was  extracted  from  fresh- 
frozen  tissue  or  formalin-fixed,  paraffin-embedded  tissue 
blocks.  All  tumors  had  >10%  neoplastic  cellularity,  as 
estimated  by  review  of  hematoxylin  and  eosin— stained 
slides.  Clinical  specimens  were  obtained  in  accordance 
with  the  Declaration  of  Helsinki  and  the  ethics  guidelines 
of  the  human  subjects  division  of  the  University  of 
Washington  (Seattle,  WA). 

Our  study  design  is  summarized  in  Figure  1.  Clinical  [FI] 
testing  for  mutations  in  BRAF ,  KRAS ,  and  NRAS  was  per¬ 
formed  in  61  specimens  by  the  University  of  Washington 
Clinical  Molecular  Genetics  Laboratory.  Variants  were 
identified  using  either  the  UW-OncoPlex  NGS  oncology  Q7 
assay3 1  {http : //tests,  labmed.  Washington.  edu/UW-OncoPlex , 
last  accessed  June  4, 2014)  or  PCR-amplification  and  melting 
curve  analysis  assays  for  each  individual  gene  target. 

Clinical  MSI-PCR  testing  of  81  colorectal  tumor  speci¬ 
mens  was  performed  by  the  University  of  Washington 
Clinical  Molecular  Genetics  Laboratory  using  the  MSI 
analysis  kit  (Promega,  Fitchburg,  WI).  Samples  demon¬ 
strating  instability  at  two  or  more  of  the  five  mono¬ 
nucleotide  markers  included  in  this  panel  were  considered 
MSI  positive  [MSI  high  (MSI-H);  diagnosis,  44  specimens]. 

All  other  specimens  analyzed  in  this  study  did  not  demon¬ 
strate  any  unstable  loci  by  MSI-PCR,  and  were  considered 
MSI  negative  (microsatellite  stable;  diagnosis,  37  speci¬ 
mens).  IHC  staining  for  MLH1,  MSH2,  MSH6,  and  PMS2 
was  performed  using  standard  diagnostic  techniques. 
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Table  1  Loci  and  Primer  Sequences 
Assay  Target  Primer  Repeat 

stage  Locus  coordinates  type  Forward  primer  sequence  Reverse  primer  sequence 


Stage  1 

Bat-25 

Chr4:  55598177- 

(A)22 

PCR 

55598271 

Stage  1 

Bat-26 

Chr2:  47641524- 

(T)  19 

PCR 

47641622 

Stage  1 

MONO- 

Chr2:  39564859- 

(T)  28 

PCR 

27 

39564957 

Stage  1 

NR-21 

Chrl4:  23652311- 

(A)22 

PCR 

23652403 

Stage  1 

NR-24 

Chr2:  95849327- 

(T)24 

PCR 

95849419 

Stage  1 

MSI-01 

Chrl:  201754376- 

(T)17 

PCR 

201754446 

Stage  1 

MSI-03 

Chr2:  62063059- 

(A)17 

PCR 

62063129 

Stage  1 

MSI-04 

Chr2:  108479588- 

(T)18 

PCR 

108479658 

Stage  1 

MSI-06 

Chr5:  172421726- 

(T)15 

PCR 

172421796 

Stage  1 

MSI-07 

Chr6:  142691916- 

(T)17 

PCR 

142691986 

Stage  1 

MSI-08 

Chr7:  1787485- 

(A)17 

PCR 

1787555 

Stage  1 

MSI-09 

Chr7:  74608706- 

(7)13 

PCR 

74608776 

Stage  1 

MSI-11 

Chrll:  106695477- 

(7)12 

PCR 

106695550 

Stage  1 

MSI-12 

Chrl5:  45897737- 

(T)14 

PCR 

45897807 

Stage  1 

MSI-13 

Chrl6:  18882625- 

(A)15 

PCR 

18882695 

Stage  1 

MSI-14 

Chrl7:  19314883- 

(T)18 

PCR 

19314953 

Stage  1 

HSPH1  - 

Chrl3:  31722570- 

(A)17 

PCR 

T17 

31722746 

Stage  1 

BRAF 

Chr7:  140453095- 

NA 

PCR 

exon  15 

140453431 

Stage  1 

KRAS 

Chrl2:  25378395- 

NA 

PCR 

exon  4 

25378686 

Stage  1 

KRAS 

Chrl2:  25380233- 

NA 

PCR 

exon  3 

25380491 

Stage  1 

KRAS 

Chrl2:  25398245- 

NA 

PCR 

exon  2 

25398504 

Stage  1 

NRAS 

Chrl:  115252168- 

NA 

PCR 

exon  4 

115252401 

Stage  1 

NRAS 

Chrl:  115258629- 

NA 

PCR 

exon  3 

115258838 

Stage  1 

NRAS 

Chrl:  115256475- 

NA 

PCR 

exon  2 

115256731 

Stage  2 

NA 

NA 

NA 

PCR 

Index 

NA 

NA 

NA 

read 

5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATC  TGG AGG ATG AC  G AGTTGGC  C  C  T AGAC - 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG 
ATC  T AGTGG AGTGG AGG AGGGG AG AG AAA - 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTCTACTGTCCTACTGTGCCTGGCTCC - 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTCTGGTGCACAGAGCAGAACCATCCT-3/ 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTCTGTAGTCCC AGCTATTCGGGAGGC - 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTTTGATGTCCTGCGTCTAGGGTCTGC - 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTGCCACTGCTATTTGAAAGAGTTGCTC - 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTTCCAAGATTCCTTCCCTGGCCACTC - 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTAGCAGCAAACTGAACAGGTCACCAAC - 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTGCTGAAAGCAACCTAAGCTGTGGTGA- 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTCCAGCCCCCATGTACACTGTAGTCG-3/ 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTGTCTCGGCTACTTGGGAGGCTTAGG-3' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTAGCATGTTTGCAGCCTTCTTCTGGA- 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTGCTGAGGCTAAAC ACTATC ATGCCA- 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTAC ATCTTCAGGTC AGGAAAACAGCTCG- 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTCATTTCAACTGACCTGCCTGGCCTC - 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTTGGAAAAGGAACTGCATCTGTGACGG- 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTAC AACTGTTC AAACTGATGGGACC - 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTTTTCAGTGTTACTTACCTGTCTTGTC - 3 ' 
S'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTCCCAGTCCTCATGTACTGGTCCCT- 3' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTTGAATTAGCTGTATCGTCAAGGCACTC - 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTAATGCTGAAAGCTGTACCATACC - 3 ' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTTGTGGCTCGCCAATTAACCCTG-3' 
5'-ACACTCTTTCCCTACACGACGCTCTTCCG- 
ATCTAGGAAGCCTTCGCCTGTCCTCA- 3 ' 

5 '  - AATG AT AC GGC G AC C AC C GAGATC T AC - 
ACTCTTTCCCTACACGACGC - 3 ' 

5'  - AGATCGGAAGAGCGGTTCAGCAGGA- 
ATGCCGCGCCCG- 3 ' 


S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCCCAAAGAGACAGCAGTTGGAACATGA- 3 ' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCTTGCAGTTTCATCACTGTCTGCGGT-3/ 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCCAGCCTGGGCAAGATAATGAGACCC - 3 ' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCGCAACCTCAAAAGCTGCCTCCCTTT-3/ 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCAAATGACCCCTTCCTGCCCATCACT-3/ 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCGACTGGAGCCTTGGACAGGTTGAGA- 3 ' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCGCC ACTGCTATTTGAAAGAGTTGCTC - 3 ' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCACTGTCTGTAGTCCTGGCTTCGTGG-3/ 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCAGC AGCAAACTGAACAGGTCACCAAC - 3 ' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCTGCTATAAGAGCTGAGCAGACGAC A- 3 ' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCCCCACCCCAAGGCCAAAATCAGTAA-3/ 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCCCGACTAAAGAGGTCATTCACTTGT-3' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCAGC ATGTTTGCAGCCTTCTTCTGGA- 3 ' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCCAGAGGTTGCAGTGAGCCGAGATTG-3' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCTAATGACTTGGGCTTTGGAAGCAGC - 3 ' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCCTTGGCAAACGGGCAAGTCTTCAGT-3' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCTTTTCCTAATCCCCTGTGAAACCTGT- 3 ' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCCCTCATCCTAACACATTTC AAGCCCCA- 3 ' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCTGACAAAAGTTGTGGACAGGT- 3 ' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCCCGTCATCTTTGGAGCAGGAACA- 3' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCACACGTCTGCAGTCAACTGGAATTT-3/ 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCCCCAGCCTAATCTTGTTTTTCTT- 3' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCTGAGAGACAGGATCAGGTCAGCGG- 3' 
S'-CGGTCTCGGCATTCCTGCTGAACCGCTCTTGTG- 
TTTCACC AGATAGGCAGAAATGGGCTTGA- 3 ' 

5 '  - C AAGC AG AAG AC GGC AT AC GAGATXXXXXXXXC G - 
GTCTCGGCATTCCTGCTGAACCG- 3 '  * 

NA 


Q20 


*X  indicates  the  presence  of  an  8-bp  sample-specific  index  sequence. 
NA,  not  applicable. 


Determination  of  mSINGS  Baseline  Reference  Values 

Determining  MSI  status  by  mSINGS  analysis  of  NGS  data 
entails  comparing  experimental  results  against  baseline 
reference  values  at  each  microsatellite  locus  to  assess  its 
instability.  Because  amplifying  microsatellite  loci  by  PCR 


generates  a  distribution  of  alternate  fragments  (stutter  arti¬ 
fact)  that  results  from  template  slippage  during  cycling,42,43 
it  is  necessary  to  establish  assay-specific  baseline  values  for 
each  locus.  To  establish  baseline  reference  values,  we 
extracted  DNA  from  42  peripheral  blood  specimens  and 
analyzed  them  using  the  MSIplus  assay.  Q8 
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Figure  1  Study  design  and  summary  of  results.  Diagnostic  testing  algorithm  and  results  are  depicted.  The  number  of  specimens  at  each  stage  is  indicated 
numerically  at  corresponding  nodes.  Fifteen  specimens  were  tested  for  both  microsatellite  instability  (MSI)  status  and  mutational  hotspots;  thus,  inclusion  in 
these  categories  is  not  mutually  exclusive.  IHC,  immunohistochemistry. 


Library  Preparation  and  Sequencing 

Sequencing  libraries  were  generated  by  PCR  amplification 
in  two  separate  stages.  The  purpose  of  the  first  stage  was  to 
simultaneously  amplify  the  loci  of  interest  from  the  genome 
and  to  incorporate  partial  Illumina  sequencing  adaptors  into 
the  amplification  product.  The  second  stage  of  PCR  fully 
extended  the  sequencing  adaptors  and  incorporated  unique 
8-bp,  sample-specific  index  sequences  (Table  1),  which 
enabled  the  multiplexing  of  multiple  specimens  together 
onto  the  same  sequencing  run.  In  addition  to  the  tumor 
samples,  control  DNA  from  HapMap  individual  NA12878 
(Coriell  Institute,  Camden,  NJ)  was  included  with  each  li¬ 
brary  preparation  and  sequencing  run,  and  served  as  a 
negative  control  for  both  MSI  status  and  RAS  and  BRAF 
mutations,  as  was  a  nontemplate  amplification  control. 

The  first  stage  of  PCR  was  performed  in  two  separate 
reactions,  one  using  an  equimolar  pool  of  the  microsatellite 
primers,  and  the  other  using  primers  targeting  mutational 
hotspots.  In  the  latter  primer  pool,  primers  were  combined 
in  equimolar  amounts,  except  for  the  primer  pairs  ampli¬ 
fying  NRAS  exon  2  (included  at  0.5  x  concentration)  and 
NRAS  exon  4  (included  at  1.5  x  concentration).  Both  first- 
stage  PCRs  were  performed  in  a  25 -pL  volume  using  the 
Qiagen  Multiplex  PCR  Kit  (Qiagen,  Valencia,  CA),  and 
incorporating  50-ng  template  DNA,  0.25  pmol/L  of  the 
appropriate  primer  pool,  and  1  x  Qiagen  Q-solution.  PCR 
cycling  for  both  primer  pools  was  as  follows:  5  minutes’ 
incubation  at  95°C;  30  cycles  of  94°C  for  30  seconds,  60°C 
for  90  seconds,  and  72°C  for  60  seconds;  and  a  final 
extension  at  72°C  for  10  minutes.  Before  the  second  stage  of 
PCR,  amplification  products  were  purified  using  a 
0.8  x  volume  of  Agencourt  AMPure  XP  magnetic  beads 
(Beckman-Coulter,  Indianapolis,  IN),  according  to  the 
manufacturer’s  instructions. 


The  second  stage  of  PCR  was  performed  using  5  ng  of 
amplification  products  from  the  first  stage  of  PCR  as  template. 
PCR  was  performed  in  a  25 -pL  volume  using  KAPA  HiFi 
HotStart  ReadyMix  PCR  Kit  (KAPABiosy stems,  Wilming¬ 
ton,  MA),  and  0.3  pmol/L  of  each  of  the  second-stage  primers 
(Table  1).  This  phase  incorporated  a  reverse  primer  carrying  a 
sample-specific  index  sequence;  however,  the  same  reverse 
primer  was  used  in  separately  amplifying  the  microsatellite 
and  mutational  hotspot  amplicons  derived  from  the  same 
specimen.  PCR  cycling  conditions  were  as  follows:  3  minutes’ 
incubation  at  95°C;  five  cycles  of  98°C  for  20  seconds,  65°C 
for  15  seconds,  and  72° C  for  60  seconds;  and  a  final  extension 
at  72°C  for  5  minutes.  PCRs  were  purified  using  a 
1.8  x  volume  of  Agencourt  AMPure  XP  magnetic  beads. 

Amplicons  derived  from  the  same  specimen  were  pooled 
in  an  8:1  volumetric  ratio  of  microsatellite  PCR  product/ 
mutational  hotspot  PCR  product  before  sequencing. 
Sequencing  was  performed  on  an  Illumina  MiSeq  (San 
Diego,  CA)  using  200-bp,  single-ended  reads  and  an  8-bp 
index  read,  with  the  addition  of  a  custom  index  read  primer 
(Supplemental  Table  SI).  A  5%  concentration  of  PhiX 
Control  version  3  (Illumina)  was  included  in  each  Q10 
sequencing  run.  Sequencing  used  a  Micro  or  Nano  MiSeq  Q11 
Reagent  version  2  300-cycle  kit  (Illumina),  depending  on 
the  number  of  samples  pooled  for  sequencing  (up  to  32 
samples  and  up  to  9  samples,  respectively). 

Data  Analysis 

Single-ended  sequence  reads  were  initially  aligned  to  the 
human  genome  (hgl9/GRCh37)  using  bwa  version  0.6.1-  Q12 
rl0444  ( http://sourceforge.net/projects/bio-bwa ,  last  accessed 
October  28,  2013)  and  SAMtools  version  0.1.18  (http:// 
sourceforge.net/projects/samtools ,  last  accessed  October  28, 
2013).  Sample-level,  fully  local  indel  realignment  was  then 
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Q'3  performed  using  Genome  Analysis  Toolkit  version  3.2  {https:// 
www.broadinstitute.org/gatk ,  last  accessed  May  28,  2014), 39,46 
followed  by  quality  score  recalibration,  to  generate  a  final, 
realigned,  and  recalibrated  alignment,  which  was  used  for 
subsequent  analyses. 

Identification  of  single-nucleotide  variants,  insertions, 
and  deletions  in  KRAS ,  NRAS ,  and  BRAF  was  performed 
Q14  using  VarScan  version  2.3.7  ( http ://var scan. sourcef or ge . 
net ,  last  accessed  May  16,  20 14), 47  with  the  minimum 
variant  frequency  set  to  0.005  reads,  the  minimum  number 
of  variant  reads  set  to  2,  and  strand  filtering  disabled.  In 
addition  to  primary  variant  calls,  we  also  tabulated  the  ab¬ 
solute  number  of  sequence  reads  matching  specific  variants 
of  clinical  actionability  (Supplemental  Table  S2)  using  the 
VarScan  readcounts  function. 

Q15  MSI  status  was  determined  using  the  mSINGS  package 
{https://bitbucket.org/uwlabmed/msings,  last  accessed  April 
27,  2015)  with  multiplier  set  to  1.75,  msi_min_threshold 
set  to  0.27,  and  msi_max_threshold  set  to  0.54. 

Results 

Sensitivity  and  Specificity  for  KRAS,  NRAS,  and  BRAF 
Mutations 

We  evaluated  the  ability  of  MSIplus  to  detect  mutations  in 
KRAS,  NRAS,  and  BRAF  for  a  panel  of  61  formalin-fixed, 
paraffin-embedded  tumor  specimens  known  to  be  positive 
for  mutations  in  these  genes  on  the  basis  of  prior  clinical 
[T2]  testing  (Table  2).  We  first  estimated  the  frequency  of  false¬ 
positive  sequence  reads  at  each  clinically  significant  site  in 
KRAS,  NRAS,  and  BRAF  by  tallying  mutant  reads 
(Supplemental  Table  SI)  at  known  wild-type  sites,  on  the 
basis  of  prior  clinical  testing.  Of  a  total  of  517,062  sequence 
reads  examined  from  a  subset  of  27  specimens,  2127 
sequence  reads  (0.4%)  carried  a  false-positive  mutation. 
Regardless,  multiple  reads  must  carry  the  same  artifact 
mutation  for  a  variant  to  be  called;  thus,  this  analysis 
overestimates  false-positivity  rate.  We,  therefore,  considered 
each  of  the  mutant  calls  independently,  which  yielded  an 
average  allele  fraction  of  0.07%  (SD,  0.17%)  for  false¬ 
positive  calls  of  any  particular  mutation.  We  set  a  mini¬ 
mum  threshold  of  2%  allele  fraction  for  calling  mutations 
using  MSIplus,  a  threshold  that  should  exclude  virtually  all 
false-positive  variant  calls  (z-score  =  1.5  x  10-25)  and  that 
defines  the  practical  limit  of  detection  for  this  assay. 

We  next  evaluated  the  assay’s  ability  to  detect  clinically 
relevant  mutations  within  the  three  target  genes  for  each  of 
61  positive  control  specimens  (Figure  1  and  Table  2). 
Average  read  depth  across  the  seven  separate  amplicons 
covering  mutational  hotspots  was  1652  reads  (interquartile 
range,  225  to  2306  reads).  We  achieved  100%  sensitivity 
Q16  [61  of  61  expected  variants  recovered;  95%  Cl,  94.1%— 

100%  by  the  Clopper-Pearson  (exact  binomial)  method]. 

No  specimen  demonstrated  a  false-positive  variant  call 
occurring  at  or  above  a  2%  allele  fraction  for  any  of  the 


Table  2  Detection  of  KRAS,  NRAS,  and  BRAF  Mutations 


Mutation 

No.  of  samples 

No.  detected 

KRAS  p.G12S 

1 

1 

KRAS  p.G12D 

5 

5 

KRAS  p.G12C 

2 

2 

KRAS  p.G12V 

6 

6 

KRAS  p.G13C 

1 

1 

KRAS  p.G13D 

8 

8 

KRAS  p.Q61H 

3 

3 

KRAS  p.K117N 

1 

1 

KRAS  p.A146V 

1 

1 

NRAS  p.G13V 

1 

1 

NRAS  p.G13R 

1 

1 

NRAS  p.G12S 

1 

1 

NRAS  p.G12D 

6 

6 

NRAS  p.Q61R 

8 

8 

NRAS  p.Q61L 

2 

2 

NRAS  p.Q61K 

1 

1 

NRAS  p.Q61H 

1 

1 

BRAF  p.V600K 

1 

1 

BRAF  p.V600E 

11 

11 

Total 

61 

61 

44  possible  nucleotide  sequence  variants  encoding  a  clinically 
relevant  mutation  (Supplemental  Table  SI).  These  results 
equate  to  a  specificity  of  100%  (2684  of  2684  true-negative 
variant  calls;  95%  Cl,  >99.3%  to  100%).  In  addition,  we 
observed  a  high  degree  of  correlation  {R2  =  0.84)  between  the 
allele  fraction  of  mutations  detected  by  MSIplus  and  the  esti¬ 
mated  mutant  allele  fraction  from  previous  clinical  testing  using 
targeted  gene-capture  NGS  methods  (subset  of  55  specimens) 
(Figure  2  and  Supplemental  Table  S2).  Approximate  linearity  [F2] 
between  these  estimates  was  observed  over  a  range  of  2.0%  to 
59.2%  estimated  variant  allele  fraction.  This  finding  suggests 
that  the  amplification  bias  or  other  artifacts  potentially  affecting 
the  calculated  allele  fraction  are  not  pronounced  in  the  MSIplus 
assay.  Furthermore,  variants  in  two  specimens  with  2.0% 
mutant  allele  frequencies  were  successfully  identified,  sup¬ 
porting  the  assay’s  theoretical  limit  of  detection. 

To  evaluate  the  reproducibility  of  mutation  detection,  we 
examined  a  subset  of  eight  control  specimens  across  three  or 
more  independent  batches  of  library  preparation  and 
sequencing  (Supplemental  Table  S3).  The  expected  muta¬ 
tions  were  recovered  in  all  36  independent  technical  repli¬ 
cates.  The  CV  for  the  estimated  allele  fraction  was  0.06  Q17 
(range,  0.03  to  0.15). 

Determination  of  MSI  Status 

We  separately  assessed  the  assay’s  ability  to  detect  the 
MSI-positive  phenotype  on  the  basis  of  mSINGS  analysis 
of  targeted  microsatellite  loci.  By  using  MSIplus,  we 
typed  a  collection  of  81  specimens  (Figure  1)  previously  Q18 
subjected  to  MSI-PCR  testing  in  our  laboratory  (44  mi¬ 
crosatellite  unstable,  or  MSI-H,  results  and  37  microsat¬ 
ellite  stable  results).  We  first  evaluated  the  reproducibility 
of  mSINGS  score  determination  (corresponding  to  the  total 
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Allele  fraction  by  targeted  gene  capture 

Figure  2  Correlation  of  allele  fraction  by  MSIplus  and  targeted  gene- 
capture  next-generation  sequencing.  Allele  fractions  are  estimated  by 
either  the  MSIplus  assay  or  the  UW-OncoPlex  targeted  gene-capture 
sequencing  panel.  The  subset  of  55  specimens  for  which  data  from  both 
assays  were  available  is  shown.  Dashed  gray  line  indicates  a  theoretical 
perfect  correlation  between  the  two  estimations. 

fraction  of  typed  micro  satellite  loci  that  are  unstable)  by 
examining  at  least  two  separately  prepared  and  sequenced 
technical  replicates  each  for  a  subset  of  48  specimens 
(Supplemental  Table  S4).  Because  several  specimens  had 
mSINGS  scores  of  0,  calculating  the  CV  was  not  meani¬ 
ngful.  On  average,  the  SD  observed  among  technical  rep¬ 
licates  was  0.07  mSINGS  units,  and  the  SD  of  values 
around  this  mean  was  0.03. 

On  the  basis  of  the  qualitative  separation  of  microsatellite 
stable  from  MSI-H  groups35  (Figure  3),  we  initially  set  an 
empirical  threshold  of  a  0.40  mSINGS  score  for  differenti¬ 
ating  MSI-positive  from  MSI-negative  tumors.  However,  in 
light  of  mSINGS  score  variability,  specimens  with  scores 
falling  close  to  this  threshold  could  be  assigned  the  wrong 
MSI  status  by  chance  alone.  To  prevent  improper  MSI 
classifications  resulting  from  assay  variability,  specimens 
with  mSINGS  values  falling  between  0.27  and  0.54 
[threshold  zb  (average  mSINGS  variability  +  standard  de¬ 
viation  of  mSINGS  variability  x  2)]  were  consequently 
considered  uninterpretable.  In  contrast,  specimens  with 
mSINGS  scores  lower  or  higher  than  those  values  could  be 
confidently  classified  as  MSI  negative  or  MSI  positive, 
respectively.  Seventy-seven  specimens  had  interpretable 
mSINGS  scores,  and  we  repeated  library  preparation  for  the 
remaining  four  samples  that  did  not.  On  retyping,  mSINGS 
scores  for  one  of  the  four  specimens  became  interpretable, 
whereas  the  remaining  three  results  remained  ambiguous. 
Thus,  the  MSI  status  of  95%  of  all  specimens  was  initially 
interpretable,  and  on  repeat  typing,  this  proportion  increased 
only  modestly,  to  96%. 


Most  MSI-negative  specimens  were  readily  distinguished 
from  MSI-positive  specimens  on  the  basis  of  mSINGS 
score,  and  mSINGS  interpretations  were  fully  concordant 
with  MSI-PCR  interpretations  in  most  cases  (Figure  3  and 
Supplemental  Table  S5).  However,  four  samples  had 
discordant  results  between  the  two  methods,  which  war¬ 
ranted  further  investigation. 

One  specimen  was  classified  as  MSI  negative  by  MSI-PCR, 
but  had  a  high  mSINGS  score  (0.63)  and  was,  therefore, 
interpreted  as  MSI  positive  by  the  MSIplus  assay.  Review  of 
laboratory  records  indicated  that,  on  subsequent  workup,  the 
tumor  was  found  to  have  deficient  MHS6  expression  by  IHC 
and  that  gene  sequencing  identified  a  germline  MSH6  muta¬ 
tion  in  the  patient,  establishing  a  diagnosis  of  Lynch 
syndrome.  We  conclude  that  this  instance  represents  a  false¬ 
negative  MSI-PCR  result,  and  that  the  correct  diagnosis  was 
achieved  by  MSIplus.  The  other  three  discrepancies  corre¬ 
sponded  to  specimens  typed  as  MSI  positive  by  MSI-PCR  but 
classified  as  MSI  negative  using  MSIplus.  In  two  cases, 
clinical  IHC  testing  was  performed  and  did  not  reveal  loss  of 
expression  in  any  mismatch  repair  proteins.  Although  MSI 
positivity  has  been  observed  in  a  background  of  normal  MMR 
protein  expression  in  approximately  5  %  of  cases, 1  these  cases 
are  most  consistent  with  false-positive  MSI-PCR  results. 
Conversely,  the  remaining  discrepant  case  likely  represents  a 
false  negative  by  MSIplus  because  loss  of  MSH2  and  MSH6 
expression  was  seen  by  IHC  and  a  germline  MSH6  mutation 
was  identified  by  genetic  testing. 

Accounting  for  these  additional  clinical  data,  and 
assuming  that  all  specimens  receiving  a  concordant 


MSI-PCR  result 

Figure  3  Inferring  microsatellite  instability  (MSI)  status  by  MSIplus. 
Data  are  stratified  according  to  testing  results  by  MSI-PCR.  The  mSINGS 
score  (the  fraction  of  interrogated  microsatellite  loci  that  are  unstable)  is 
plotted  for  each  specimen.  The  dashed  lines  at  mSINGS  scores  of  0.27  and 
0.54  indicate  the  cutoffs  for  delineating  MSI-positive  and  MSI-negative 
specimens  by  MSIplus:  mSINGS  scores  falling  below  these  values  were 
interpreted  as  negative,  scores  above  those  values  were  interpreted  as 
positive,  and  scores  falling  between  the  values  could  not  be  reliably 
interpreted.  Arrowheads  indicate  specimens  misclassified  by  MSI-PCR; 
asterisk  indicates  specimen  misclassified  by  MSIplus,  as  resolved  by 
alternative  clinical  testing  results  (immunohistochemical  and/or  genetic 
testing).  For  specimens  that  were  typed  multiple  times,  one  representative 
mSINGS  score  is  displayed. 
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diagnosis  by  MSI-PCR  and  MSIplus  were  correctly 
assigned  an  MSI  status,  we  conclude  that  the  validation 
set  used  in  our  study  comprised  a  total  of  36  MSI-positive 
cases  and  42  MSI-negative  cases.  We  calculate  that 
MSIplus  has  an  overall  sensitivity  of  97.1%  (34  of  35 
MSI-positive  cases  identified;  95%  Cl,  85.1%— 99.9%) 
and  a  specificity  of  100%  (42  of  42  MSI-negative  cases 
identified;  95%  Cl,  91.6%— 100%)  for  determining  tumor 
MSI  status.  The  positive  predictive  value  and  negative 
predictive  value  for  MSIplus  were  97.1%  and  97.7%, 
respectively.  For  the  same  panel  of  test  specimens,  MSI- 
PCR  demonstrated  a  sensitivity  of  97.1%  (34  of  35  MSI- 
positive  cases  identified;  95%  Cl,  85.1%— 99.9%),  a 
specificity  of  95.2%  (40  of  42  MSI-negative  cases  iden¬ 
tified;  95%  Cl,  83.8%— 99.4%),  a  positive  predictive 
value  of  94.6%,  and  a  negative  predictive  value  of  95.6% 
for  the  same  set  of  validation  specimens. 

Correlation  of  MSIplus,  MSI-PCR,  and  IHC 

Although  our  primary  objective  was  to  compare  the  results 
of  molecular  MSI  testing  by  MSI-PCR  and  MSIplus,  we 
additionally  correlated  our  findings  with  IHC  for  40  cases 
where  this  information  was  available  (Figure  1  and 
Supplemental  Table  S5).  The  diagnosis  rendered  by  each 
of  the  three  approaches  (IHC,  MSI-PCR,  and  MSIplus) 
was  fully  concordant  in  32  of  37  cases  where  MSIplus  was 
interpretable.  One  discrepant  case,  discussed  above,  evi¬ 
denced  isolated  loss  of  MSH6  by  IHC,  and  was  negative 
by  MSI-PCR  but  positive  by  MSIplus,  suggesting  that 
MSIplus  and  IHC  obtained  the  proper  diagnosis.  Three 
additional  cases  receiving  a  MSI-positive  diagnosis  by 
both  MSI-PCR  and  MSIplus  showed  no  loss  of  MMR 
protein  expression  by  IHC,  suggesting  a  MSI-positive 
status  without  loss  of  MMR  protein  expression.  One 
specimen  was  negative  by  MSI-PCR  and  MSIplus,  but 
demonstrated  reduction  of  MLH1  and  PMS2  expression 
by  IHC:  it  is  possible  that  this  latter  case  represents  either 
a  false-negative  molecular  result  or  a  false-positive  result 
by  IHC.20 

Quantitative  Correlation  of  MSIplus  and  Conventional 
MSI-PCR 

Last,  we  examined  whether  there  was  a  correlation  between 
the  fraction  of  unstable  markers  characterized  by  MSIplus 
(ie,  the  mSINGS  score)  and  the  fraction  of  unstable  mi- 
[F4]  crosatellite  loci  detected  by  MSI-PCR  (Figure  4).  We 
excluded  from  analysis  specimens  with  discordant  results 
between  the  two  assays  and  those  three  having  non¬ 
interpretable  mSINGS  scores. 

Overall,  the  two  measurements  of  MSI  demonstrated  a 
strong,  positive  correlation  ( R 2  =  0.89).  Of  42  specimens 
with  no  unstable  loci  detectable  by  MSI-PCR,  20  demon¬ 
strated  nonzero  mSINGS  scores;  nevertheless,  a  trend  of 
mSINGS  overestimating  the  fraction  of  unstable  loci  was 


Fraction  loci  unstable  by  MSI-PCR 

Figure  4  Correlation  of  locus  instability  measured  by  microsatellite 
instability  (MSI)-PCR  and  MSIplus.  The  mSINGS  score  (the  fraction  of 
interrogated  microsatellite  loci  that  are  unstable)  is  plotted  against  the 
fraction  of  unstable  loci  determined  by  MSI-PCR.  Overall  mSINGS  score  for 
the  MSIplus  assay  is  shown  separately  from  mSINGS  score  calculated  for  the 
five  loci  in  common  with  MSI-PCR.  Dashed  gray  line  indicates  a  theoretical 
perfect  correlation  between  the  two  measures  of  instability. 

not  consistent  across  all  specimens  deemed  MSI  positive  by 
MSI-PCR.  We  noted  that  if  MSI-negative  MSI-PCR  results 
were  removed  from  consideration,  the  correlation  of 
mSINGS  score  and  the  fraction  of  unstable  loci  observed  for 
MSI-H  specimens  was  not  statistically  meaningful 
C R 2  =  9.1  x  10  3).  This  finding  suggests  that,  above  the 
threshold  for  delineating  an  MSI-positive  phenotype,  the 
fraction  of  unstable  markers  identified  by  MSI-PCR  cannot 
be  generalized  to  infer  the  degree  of  overall  genomic 
instability. 

We  separately  evaluated  correlation  between  the  fraction 
of  unstable  loci  identified  by  either  MSI-PCR  or  MSIplus 
for  the  subset  of  five  loci  that  are  represented  in  both  assays. 
Again,  a  strong  positive  correlation  was  observed  between 
these  two  metrics  ( R 2  =  0.80),  although  the  lack  of  identity 
between  them  indicates  subtle  differences  in  individual  loci 
being  scored  as  stable  or  unstable  between  the  two  assays. 

Discussion 

Because  the  number  of  single-gene  molecular  tests  needed 
to  adequately  characterize  tumors  continues  to  increase,  the 
practical  consequences  of  increased  health  care  costs, 
increased  test  turnaround  times,  and  the  potential  to  deplete 
available  tissue  material  during  the  course  of  testing  become 
an  increasing  concern.  The  use  of  highly  scalable  NGS 
technologies  has  proved  a  means  to  overcome  this  chal¬ 
lenge,  in  many  cases  improving  the  quality  and  capabilities 


The  Journal  of  Molecular  Diagnostics  ■  jmd.amjpathol.org  7 

FLA  5.4.0  DTD  ■  JMDI435_proof  ■  12  September  2015  ■  8:25  pm  ■  EO:  JMD15_0081 


869 

870 

871 

872 

873 

874 

875 

876 

877 

878 

879 

880 

881 

882 

883 

884 

885 

886 

887 

888 

889 

890 

891 

892 

893 

894 

895 

896 

897 

898 

899 

900 

901 

902 

903 

904 

905 

906 

907 

908 

909 

910 

911 

912 

913 

914 

915 

916 

917 

918 

919 

920 

921 

922 

923 

924 

925 

926 

927 

928 

929 

930 


931 

932 

933 

934 

935 

936 

937 

938 

939 

940 

941 

942 

943 

944 

945 

946 

947 

948 

949 

950 

951 

952 

953 

954 

955 

956 

957 

958 

959 

960 

961 

962 

963 

964 

965 

966 

967 

968 

969 

970 

971 

972 

973 

974 

975 

976 

977 

978 

979 

980 

981 

982 

983 

984 

985 

986 

987 

988 

989 

990 

991 

992 


ARTICLE  IN  PRESS 


Hempelmann  et  al 


of  molecular  testing.  In  light  of  these  considerations, 
we  developed  MSIplus,  an  NGS  assay  for  characterizing 
colorectal  tumor  specimens  that  integrates  extended  RAS 
and  BRAF  gene  testing  with  MSI  analysis  and,  therefore, 
encompasses  recommended  molecular  testing  guidelines  in 
a  single  assay. 

MSIplus  had  an  overall  sensitivity  of  100%  and  a 
specificity  of  100%  for  identifying  cancer-associated 
mutations  in  BRAF ,  KRAS,  and  NRAS  genes,  as  calcu¬ 
lated  from  a  panel  of  61  specimens  carrying  a  spectrum  of 
known  mutations  (Table  2).  The  ability  of  a  molecular 
assay  to  detect  mutations  present  in  a  tumor  specimen  is 
dependent,  in  part,  on  the  proportion  of  neoplastic  cells  in 
the  sample  and  the  fraction  of  total  tumor  cells  that  carry 
the  mutation  of  interest.  1  The  limit  of  detection  for 
MSIplus  is  a  2%  mutant  allele  fraction,  lower  than  other 
prevalent  clinical  diagnostic  methods.  Because  the 
limit  of  detection  in  this  assay  partially  reflects  the  rate  of 
false-positive  sequence  reads,  the  sensitivity  of  MSIplus 
for  low-prevalence  mutations  could  potentially  be 
improved  in  future  iterations  through  practices  such  as 
incorporating  molecular  tagging-mediated  sequencing 
error  correction.  However,  such  methods  present 

technical  challenges  and  could  negatively  affect  other 
aspects  of  the  assay’s  performance. 

Because  it  is  a  sequencing-based  approach,  MSIplus 
should  be  able  to  identify  most  single-nucleotide  poly¬ 
morphisms  or  small  insertions  or  deletions  occurring 
within  the  selected  target  regions.  Compared  with  real¬ 
time  PCR  or  melt-curve  assays,  which  are  designed  to 
detect  only  a  specific  subset  of  actionable  mutations, 
MSIplus  is  capable  of  identifying  variants  without  prior 
knowledge  or  expectation  of  the  underlying  genetic 
change.  This  property  greatly  improves  the  ability  of 
MSIplus  to  identify  rare  or  unusual  clinically  significant 
nucleotide  alterations. 

MSIplus  had  97.1%  sensitivity  and  100%  specificity  for 
characterizing  the  MSI  phenotype  in  a  test  set  of  78  tumors, 
performance  characteristics  that  were  similar  to  those  re¬ 
ported  for  MSI-PCR  in  other  studies5  and  slightly  better 
than  the  performance  of  MSI-PCR  for  the  same  set  of 
specimens.  Although  MSIplus  does  not  substantially 
improve  sensitivity  or  specificity  compared  with  MSI-PCR, 
it  has  important  advantages  over  existing  assays.  Evaluation 
of  MSI  by  mSINGS  analysis  eliminates  the  need  for 
matched  normal  patient  material,  which  is  typically 
required  for  conventional  MSI-PCR,  thereby  expanding  the 
scope  of  available  patient  specimens  that  can  be  successfully 
typed  using  MSIplus.  The  interpretation  of  MSIplus  uses  an 
automated  analysis  pipeline  that  is  based  on  quantitative, 
descriptive  statistics.  This  feature  may  improve  the  consis¬ 
tency  of  MSI  diagnosis  and  reduce  interlaboratory  and 
intralaboratory  variation. 

We  were  unable  to  confidently  call  MSI  status  using 
MSIplus  for  a  small  subset  of  specimens  (95%  of  specimens 
were  interpretable  with  the  first  round  of  testing,  whereas 


only  5%  were  not).  Three  of  the  four  specimens  with 
indeterminate  mSINGS  scores  again  yielded  uninterpretable 
results  after  repeated  testing,  suggesting  that  the  assay 
cannot  confidently  type  MSI  status  for  a  small  fraction  of 
samples.  Although  the  biological  significance  of  these 
persistently  indeterminate  mSINGS  scores,  if  any,  is  un¬ 
clear,  1  studies  examining  larger  numbers  of  microsatel- 
lite  markers  have  suggested  that  such  cases  do  not  represent 
a  distinct  disease  category  of  subtype  of  MSI.  These 
indeterminate  specimens  could  potentially  be  resolved  by 
increasing  the  number  of  microsatellite  loci  examined  in  the 
assay  to  enable  more  accurate  assessment  of  the  mSINGS 

35 

score. 

Our  study  identified  several  cases  where  IHC  and  mo¬ 
lecular  testing  were  discordant,  as  expected.  It  is  known 
that  isolated  MSH6  deficiencies  may  result  in  false-negative 
MSI-PCR  results,5  and  our  sample  cohort  contained  one 
such  case.  Unlike  MSI-PCR,  MSIplus  identified  this  MSH6- 
dehcient  specimen  as  MSI  positive,  rendering  a  proper 
molecular  diagnosis.  Moreover,  ah  indeterminate  MSIplus 
results  occurred  in  cases  where  IHC  indicated  reduced 
MMR  protein  expression,  including  two  testing  negative  by 
MSI-PCR.  Although  anecdotal,  these  findings  suggest  an 
improved  ability  of  MSIplus  to  identify  MSH6-mutated 
specimens  and  improved  sensitivity  for  detecting  MMR 
pathway  deficiencies  compared  with  MSI-PCR.  With 
further  refinement,  the  MSIplus  assay  may  offer  significant 
performance  advantages  over  MSI-PCR  for  such  cases. 

We  anticipate  that  MSIplus  will  prove  useful  in  charac¬ 
terizing  colorectal  cancers  while  potentially  reducing  oper¬ 
ating  costs  and  standardizing  the  interpretation  of  testing 
results.  The  assay  is  rapid,  and  compatible  with  a  2-  or  3-day 
turnaround  time:  library  preparation,  approximately  8  hours; 
sequencer  setup,  approximately  1  hour;  sequencing, 
approximately  8  hours;  data  analysis,  approximately  2  hours 
per  specimen  on  a  four-processor  machine,  but  scalable  on 
larger  computing  systems  to  process  multiple  specimens 
simultaneously.  Because  the  assay  uses  PCR-mediated  li¬ 
brary  preparation,  it  requires  minimal  input  DNA  (50  ng) 
and  should  function  even  for  partially  degraded  specimens. 
The  modular  nature  of  the  multiplexed  primer  design  will 
enable  relatively  straightforward  expansion  of  the  assay  to 
additional  molecular  targets,  as  necessary,  in  response  to 
future  diagnostic  requirements.  The  approach  of  focused, 
integrated  NGS  testing,  tailored  to  a  specific  tumor  type  or 
diagnostic  workflow,  is  a  powerful  paradigm  that  can  be 
adapted  to  other  clinical  scenarios,  and  will  become  more 
common  as  NGS  technologies  are  increasingly  integrated 
into  clinical  laboratories. 
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